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Résumé 


L’érosion du sol est un phénomène naturel qui rend compte de certains processus 
(détachement des particules, transport et dépôt) induits par différents agents érosifs. Dans 
les zones semi-arides, les précipitations sont souvent intenses sur les sols secs à végétation 
pauvre, ce qui peut entraîner une forte susceptibilité à l'érosion des sols. Par conséquent, 
une évaluation appropriée de l'érosivité des précipitations revêt une importance particulière 
en raison des effets négatifs causés par l'épuisement de la couche superficielle du sol et la 
charge excessive de sédiments dans les eaux réceptrices des réservoirs. 

La présente étude a été menée sur le bassin versant de l’Oued K'sob situé au Nord-Est 
de l'Algérie, dans le but d’étudier l’influence la variabilité climatique sur la production des 
sédiments. Le bassin draine une superficie de 1480 km? située entre les altitudes 585 et 
1888 m. L’Oued principal s’écoule sur une longueur 73 km et alimente le barrage K’sob 
d’une capacité initiale de 29.5 Mm3, mis en service en 1940. La région d’étude présente un 
climat semi-aride a tendance continentale avec un hiver relativement pluvieux et un été sec 
et chaud, avec une pluie moyenne interannuelle est de 340 mm. La construction des 
courbes IDF a permis de déterminer un exposant climatique de la région, b=0.75. 

Pour l’érosivité des précipitations, le calcul d'un tel indice est basé sur les pluies dépassant 
un seuil spécifique et nécessite des données pluviométriques avec une fine résolution 
temporelle, qui, souvent, sont rares ou difficiles à acquérir. L'examen des pluies 
quotidiennes survenant avant une inondation portant une charge sédimentaire a montré une 
variabilité spatiale des seuils d'érosivité des précipitations. Les valeurs saisonnières 
des seuils sont faibles et se situent entre 2 mm en été et 6 mm en hiver, ce qui met 
en évidence un processus d'érosivité typiquement élevé dans les régions semi-arides. 

Des relations empiriques, établies à l'échelle saisonnière, ont été proposées comme solution 
alternative au calcul de l'indice R dérivé de l'équation révisée des pertes en sol. 
Les modèles déterminés ont permis de simuler l'érosivité des événements pluvieux 
en fonction de la pluie quotidienne. Entre 68 et 78 % de la variance de l'érosivité 
des précipitations s'explique par la pluie quotidienne donnant lieu à un événement pluvieux 
érosif. Ensuite, la moyenne spatiale de l'indice d'érosivité annuel a fluctué entre 228 et 386 
MJ mm ha-1 h-1 an-1 avec une moyenne interannuelle de 302 MJ mm ha-1 h-1 an-1, 
ce qui a sous-estimé de 6% l'indice d'érosivité quantifié selon l'équation universelle révisée 
de perte de sol. 

L'érosivité des précipitations est le facteur déterminant du rendement en sédiments avec un 
degré de fixation différent au cours de l'année. En automne, 68 % de la variance de la 
production de sédiments s'explique par l'érosivité des pluies, contre seulement 42 % au 
printemps en raison des changements des conditions du sol, notamment la présence d'une 
couverture végétale qui protège le sol contre l'érosivité des pluies. 


Mots-clés : Précipitations, indice d'érosivité, transport de sédiments, semi-aride, bassin 
versant du Wadi K'sob, Nord-Est de l'Algérie. 


Abstract 


Soil erosion is a natural phenomenon that accounts for certain processes (particle 
detachment, transport and deposition) induced by different erosive agents. In semiarid 
areas, rainfall is often intense on dry soils with poor vegetation which might lead to high 
soil erosion susceptibility. Therefore, an appropriate assessment of the rainfall erosivity is 
of particular importance due to negative effects caused by top soil depletion and excessive 
sediment loading to receiving waters of reservoirs. 

The present study was carried out on Wadi K'sob watershed located in northeastern 
Algeria, in order to study the influence of climate variability on sediment yield.. The basin 
drains an area of 1480 km? located between altitudes 585 and 1888 m. The main wadi 
flows over a length of 73 km and supplies the K'sob dam with an initial capacity of 29.5 
Mm*, commissioned in 1940. The study area presents a semi-arid climate with a 
continental tendency with a relatively rainy winter and a dry and hot summer, with an 
interannual average rainfall of 340 mm. The construction of the IDF curves determined a 
regional climatic exponent b=0.75. 

with regard to rainfall erosivity, the calculation of such an index is based on rains 
exceeding a specific threshold and requires rainfall data with a fine temporal resolution, 
which, often, are rare or difficult to acquire. The examination of daily rains occurring 
before a flood event carrying sediment load showed a spatial variability of the thresholds 
of rainfall erosivity. The seasonal values of the thresholds are low and lying between 2 mm 
in summer and 6 mm in winter highlighting an erosivity process characteristically high in 
semi-arid regions. 

Empirical relationships, established at seasonal scale, were proposed as an alternative 
solution to the R-index calculation derived from the Revised Soil Loss Equation. 
The determined models allowed us to simulate the erosivity of rainfall events as a function 
of daily rain. Between 68 and 78% of the variance of rainfall erosivity is explained by the 
daily rainfall giving rise to an erosive rainfall event. Then, the spatial mean of the annual 
erosivity index fluctuated between 228 and 386 MJ mm ha’! h! y‘! with an interannual 
average of 302 MJ mm ha! h?! yt, which underestimated by 6% the erosivity index 
quantified according to the Revised Universal Soil Loss Equation. Rainfall erosivity is the 
determining factor in sediment yield with a different degree of binding during the year. 
In autumn, 68% of the variance in sediment production is explained by rainfall erosivity, 
compared to only 42% in spring due to changes in soil conditions, including the presence 
of a vegetation cover that protects the soil against rainfall erosivity. 


Keywords: Rainfall, erosivity index, sediment transport, semiarid, Wadi K’sob watershed, 
Northeast Algeria. 
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INTRODUCTION GENERALE 


Dans les zones semi-arides, généralement couvertes d'une végétation clairsemée pendant 
la majeure partie de l'année, les précipitations sont souvent intenses et responsables 
d'une forte érosion et d'un apport important de sédiments (Langbein et Schumm 1958 ; 
Scott et al. 2006). Par conséquent, une évaluation appropriée de l'érosivité 
des précipitations revêt une importance particulière en raison des effets négatifs causés 
par l'épuisement des nutriments et de la matière organique (Vörösmarty et al. 2003), 
et de la charge excessive de sédiments dans les eaux réceptrices des réservoirs (Dutta 
2016). En Algérie, où les ressources en eau sont rares et irrégulières, la production 
de sédiments est l'une des plus élevées du monde (Probst et Amiotte-Suchet 1992 ; 
Megnounif et al. 2003 ; Achite et Ouillon 2016). Selon Meddi et al. (2016), les fortes 
pluies menacent 14 millions d'hectares de terres arables, avec des conséquences qui 
affectent le stockage de l'eau. En moyenne, 45 millions de m° d’envasement sont déposés 
chaque année au fond des barrages, ce qui entraîne une réduction de la capacité de 


stockage estimée à 0,7 % par an (Remini et al. 2008). 


De nombreuses études examinent l'indice d'érosité, R, présenté dans différentes versions 
de l'équation universelle de perte de sol, USLE (Wischmeier et Smith 1958) pour évaluer 
l'agressivité des précipitations comme agent causal de l'érosion des sols (Fournier 1969). 
Wischmeier et Smith (1978) ont décrit l'érosivité d'un événement pluvieux, R, comme 
le résultat d'une quantité d'énergie, Æ, libérée lorsque les gouttes de pluie touchent le sol 
et de l'intensité maximale, /30, d'un événement pluvieux pendant une période de 30 
minutes. Par conséquent, un événement pluvieux est considéré comme érosif lorsque 
la pluie est efficace en termes de rendement sédimentaire et que les caractéristiques 
suivantes sont vérifiées : (i) le cumul des pluies dépasse le seuil T=/2,7 mm, sans que 
le cumul des pluies soit inférieur à 7/2=6,35 mm pendant 6 heures ; ou (ii) le cumul des 
pluies dépasse 7/2=6,35 mm en 15 minutes. Dans ces conditions, l'évaluation de l'érosivité 
des précipitations nécessite des données pluviométriques avec une résolution temporelle 


fine, à des intervalles d'une demi-heure ou moins. 


Néanmoins, dans le monde entier, ces séries sont souvent courtes, difficiles à acquérir 


et longues à traiter car elles nécessitent l'analyse de documents provenant de pluviographes 
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(Bertoni et Lombardi 1985). En réponse à ces contraintes, des solutions alternatives sont 
proposées lorsque l'indice R est expliqué par des données pluviométriques facilement 
accessibles couvrant de longues périodes telles que les précipitations quotidiennes 
(Richardson et al. 1983 ; Choukri et al. 2016 ; Xie et al. 2016 ; Begueria et al. 2018 ; 
Bouderbala et al. 2019) ou mensuelles (Renard et Freimund 1994 ; Terranova et Gariano 
2015 ; Ballabio et al. 2017) ou annuelles (Heusch 1970 ; Bonilla et Vidal 2011 ; Yin et al. 
2015 ; Toubal et al. 2018). Oduro-Afriyie (1996), Silva (2004), et Fernandez et al. (2019) 
ont associé l'érosivité des précipitations à d'autres caractéristiques des précipitations telles 
que l'indice de Fournier, ZF (Fournier 1961), ou l'indice de Fournier modifié, IMF 
(Amoldus 1980), tandis que Choukri et al. (2016) et Farhan et Alnawaiseh (2018) ont 
souligné la variabilité spatio-temporelle de l'érosivité des précipitations en utilisant les 


précipitations annuelles, saisonnières et mensuelles. 


Les études de dégradation des sols montrent que la charge sédimentaire transportée par le 
débit naturel des cours d'eau est étroitement liée à l'érosivité des précipitations (Hicks et al. 
2000 ; Lal 2001) et que le seuil de 12,7 mm présenté par (Wischmeier et Smith 1978) est 
souvent excessif (Bollinne 1978 ; Yin et al. 2017). Denis et al. (2013) ont adopté un seuil 
de précipitations érosives de 2,5 mm pour une région semi-aride de l'Inde. Mannaerts et 
Gabriels (2000) ont utilisé un seuil de 9 mm pour caractériser les précipitations érosives 
au Cap-Vert et Xie et al. (2016) ont estimé le seuil à 9,7 mm pour la Chine. Selon 
McGregor et al. (1995) et Xie et Liu (2002), une diminution du seuil de précipitation 
érosive génère une augmentation non significative de la quantification de la production de 


sédiments. 


Le nord de l'Algérie est une région montagneuse caractérisée par de petits et moyens 
bassins de moins de 10 000 km? avec une forte variabilité spatio-temporelle de l'intensité 
de l'érosion en réponse à une végétation clairsemée et un climat irrégulier et agressif 
(Ghenim et Megnounif 2016). Outre l'irrégularité du climat, la disparité observée entre les 
différentes estimations de l'érosivité des précipitations citées dans la littérature est 
principalement due à la grande panoplie d'approches et de modèles utilisés pour estimer ce 
paramètre. En outre, la plupart des modèles sont adaptés à des études locales, dans des sites 
éloignés et souvent en négligeant les conditions pédologiques et climatiques (Meddi et al. 


2016 ; Hasbaia et al. 2017 ; Toubal et al. 2018 ; Bouderbala et al. 2019). 
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Cette thèse vise à examiner la relation entre l'érosivité des précipitations et la production de 
sédiments dans le bassin de l'oued K'sob (1480 km?), qui alimente le barrage de K'sob 
d'une capacité initiale de 29,5 Mm*, mis en service en 1940. L'étude bathymétrique, 
réalisée en 2010 par l'ANBT (Agence nationale des barrages et transferts), a révélé un taux 
élevé d'alluvions de barrage estimé à 67,6%. L'étude s’est focalisée sur l’effet de la 
variabilité du milieu environnant sur l’ampleur de l’érosion. A cet effet, l'importance des 
précipitations dans le processus érosif est analysée à l’échelle saisonnière en examinant le 
seuil effectif de précipitations par rapport à la production de sédiments. L’aboutissement 
de cette thèse est de développer un modèle permettant d’estimer l’érosivité des pluies dans 
la région d’étude à partir de données de précipitations journalières et en tenant compte de 


la variabilité saisonnière. 
Organisation de la thèse 
Cette thèse est rédigée sous la forme de cinq chapitres : 


Le premier est consacré à une synthèse bibliographique. Il donne un aperçu général et des 


notions de base sur l’érosion et le transport des sédiments dans les bassins versants. 


Le deuxième chapitre est une présentation de la zone d’étude Oued K’sob avec une 


description de ses caractéristiques physiques. 


Le troisième chapitre traite les aspects hydrologiques et climatiques de la région d’étude. 
Une analyse statistique des différents paramètres météorologiques (Pluie, Température, 


Evapotranspiration, Humidité et Vent) est effectuée. 


Le quatrième chapitre est dédié à l’étude du transport solide dans le bassin Oued K’sob, en 
se basant sur les données hydrométriques mesurées au niveau de la station de jaugeage 


Medjaz 


Le dernier chapitre, aborde l’étude de l’érosivité des précipitations et son rôle dans la 


genèse des sédiments. 
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NOTA : 
La réalisation de cette thèse est basée sur : 


o Les données hydrométriques, pluviométriques, enquêtes menées sur le terrain par 
les services de l’Agence Nationale des Ressources Hydrauliques (ANRH) et 
l'Office National de la Météorologie ONM ; 

o Cartes d’état-major & photos satellitaires (ASTER GDEM, 30m) ; 

o Une bibliographie récente et rapport interne de différents organismes liés à la 


gestion de la ressource en eau, les terres agricoles et l’aménagement du territoire. 


Les données de mesure couvrant une longue période, permettent d’examiner la sensibilité 
des bassins versants quant à la variabilité climatique à différentes échelles temporelles. 
L’étude s’appuiera sur l’outil statistique et emploi d’un système d’information 


géographique. 


Chapitre I 


Notions et concepts sur l’érosion 
et le transport des sédiments 
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Chapitre I : Notions et concepts sur 


l'érosion et le transport des sédiments 


I.1 Notions préliminaires 

L.1.1 Le Bassin versant 
Un bassin versant ou un bassin hydrographique est une zone de terre drainé par un cours 
d'eau et ses affluents (CE, 2000), (Figure I.1). L'ensemble des eaux qui tombent sur cet 


espace convergent vers un même point de sortie appelé exutoire. 


Sources 


Tributaires 


Sous-bassin —__ 
versant Na 
Ligne de partage 
des eaux (crêtes) 


x Exutoire 
A Confluent 


Figure I.1 Schéma d’un bassin versant (Source : Office national de 
l'eau et des milieux aquatiques ONEMA , France) 

Dans un bassin versant, l'eau se fraye des chemins sur et dans le sol et prend en charge des 
particules solides comme des fractions minérales, de la matière organique ainsi que de la 
matière dissoute. 

Les actions menées en amont d’un bassin versant, comme les aménagements inadaptés et les 
mauvaises pratiques agricoles peuvent fortement accentuer ou diminuer la qualité des eaux 
du ruissellement. Ainsi, les pratiques à risques se répercutent, généralement, en aval et la 
multiplication de petites perturbations peuvent modifier considérablement la réponse 


hydrologique en aval en termes de volume d’eau et matières associées. 
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1.1.2 Le cycle hydrologique 

Le cycle de l'eau, appelé aussi cycle hydrologique, est un phénomène naturel qui 
représente le parcours entre les grands réservoirs d'eau liquide, solide ou de vapeur d'eau sur 
Terre à partir des océans, de l'atmosphère, des lacs, des cours d'eau et nappes d'eaux 
souterraines ou encore des glaciers (Figure I.2). Le « moteur » du cycle est l'énergie solaire 
qui favorise l'évaporation de l'eau et entraîne l’ensemble des cheminements que peut suivre 
une particule d'eau (Laborde 2009). Ces mouvements, accompagnés de changements d'état, 
peuvent s'effectuer dans l'atmosphère, à la surface du sol et dans le sous-sol. Chaque 
particule n'effectue qu'une partie de ce cycle et avec des durées variables. Une goutte de 
pluie peut retourner à l'océan en quelques jours alors que sous forme de neige, en montagne, 


elle pourra mettre des dizaines d'années (Laborde, 2009). 


Le cycle hydrologique 


Condensation 


= 
Evaporation 


Ecoulement souterrain 


Figure I.2 Le cycle hydrologique de l’eau 
(Source : https://www150.statcan.gc.ca/n1/pub/16-201-x/2017000/sec-1/fig/fig1.1-fra.htm) 


1.1.3 Les précipitations 

Sous l’impulsion des vents, les nuages se déplacent dans l’atmosphère. Lors d’un 
changement climatique et par effet de gravité, les nuages s’alourdissent et retombent sur le 
sol sous forme d’eaux pluviales, de grêle ou de neige, 79 % des précipitations tombent sur 
les océans et les 21 % restants sur la terre et alimentent les nappes phréatiques, soit par 
infiltration, soit par ruissellement. Un peu moins de la moitié des précipitations va servir à 


recharger les nappes souterraines, le reste repart en évaporation (Musy et al, 2014). L’eau 
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qui ne parvient pas à s’infiltrer directement dans le sol, ruisselle le long des pentes pour se 
jeter ensuite dans les lacs et les rivières. Cette eau suivra leurs cours pour rejoindre les mers 
et les océans. Toutes les eaux issues du ruissellement (ruisseaux, rivières, fleuves, lacs) sont 
appelées cours d’eau de drainage. L’eau dans la nappe peut stagner jusqu’à des milliers 


d’années avant de retourner dans les océans. 


1.1.4 L'évaporation 
L’évaporation est un phénomène de changement de phase lent de l’état liquide vers l’état 
gazeux d’un liquide à une température inférieure à la température d’ébullition à la pression 
ambiante (Forestier, 2011). 
Les facteurs qui conditionnent l'évaporation sont les suivants : les rayonnements solaires et 
atmosphériques, la température de l'eau et de l'air, l'humidité de l'air, la pression 
atmosphérique, le vent, la profondeur et la dimension des nappes d'eau, la qualité de l'eau et 


les caractéristiques du bassin (exposition des versants au soleil, au vent, pentes, sol...etc). 


L.1.5 L'interception et le stockage dans les dépressions 

L’interception est le processus par lequel une partie des précipitations est captée et 
retenue par la végétation puis évaporée sans avoir atteint la surface du sol. Son importance 
et souvent marginale et donc souvent négligée sous climat semi-aride où la végétation est 
éparse. Le stockage est l'eau retenue dans les creux et les dépressions du sol pendant et après 
une averse. Comme l'interception, le stockage dans les dépressions est, souvent associé aux 
pertes. 

La rétention d'eau varie considérablement. Elle peut atteindre jusqu'à 30% de la 
précipitation totale pour une forêt mixte, 25% pour les prairies et 15% pour les cultures 
(André et Christophe, 2004). L'effet respectif de l'interception et du stockage dans les 
dépressions est très variable et diminue au cours de l'averse. Il provoque en général un retard 
dans le démarrage et la réaction hydrologique qui peut être perçue à l'exutoire du bassin. 

1.1.6 L'infiltration 
En hydrologie, l'infiltration est le processus par lequel les précipitations ou l'eau pénètrent 
dans les sols souterrains et se déplacent dans les roches à travers les fissures et les pores. La 
majeure partie de l'eau de pluie et de la neige fondue finit par s'infiltrer. 

L'eau peut être absorbée par le sol et y rester longtemps jusqu'à ce qu'elle s'évapore 


progressivement. Si la végétation est abondante, l'eau infiltrée puis absorbée par les racines 
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des plantes et transpirée plus tard. L'infiltration se produit dans les couches supérieures du 
sol mais peut aussi se poursuivre plus bas dans la nappe phréatique. 

1.1.7 Les écoulements 
La pluie ou la fonte des neiges excessives peuvent produire un écoulement de surface vers 


" 


les ruisseaux et les fossés. En hydrologie de surface, le terme " écoulement " désigne 
exclusivement la circulation de l'eau dans le réseau hydrographique. Il s'agit d'un phénomène 


qui peut se quantifier par des mesures directes de débits. 


En dehors du réseau hydrographique, on parle de ruissellement, pour qualifier un écoulement 


de surface se produisant exclusivement sur les versants. 


1.1.8 Le sol 

Le sol joue un rôle prépondérant sur le cycle de l'eau et sur l’émission des éléments 
chimiques vers l’eau. Il joue un rôle sur le partage entre infiltration et ruissellement. Il joue 
le rôle de réserve en eau pour les plantes et de milieu pour les êtres vivants présents dans le 
sol. Il est au cœur des transferts d’eau entre l'atmosphère, les nappes d'eau souterraine et les 
cours d’eau. Le sol stocke et transforme de nombreux produits épandus au sol, qu'ils soient 
solides ou liquides, minéraux ou organiques, d’origine agricole, industrielle ou urbaine. Les 
sols des zones humides, saturés en eau en hiver, ont un rôle tampon d’un point de vue 
hydrologique, permettant l’extension latérale de la nappe, et hydrochimique, favorisant les 
transformations biogéochimiques en lien avec les battements de nappe. 


I.2 Le concept “ érosion ” 


L’érosion est un processus physique et chimique par lequel le sol et les roches de la 
croûte terrestre sont continuellement soumis a une abrasion et à une corrosion.. Selon le 
domaine scientifique, le mot “ érosion” prend différent sens. Dans le damaine de 
Vhydrologie, l’érosion du sol correspond au “détachement ”, et au “transport ” des particules 
de sol, de leurs emplacements d’origine à un endroit de “dépôt” à l’aval. Les mécanismes 
d’érosion sont conditionnés par la gravité, l’eau et le vent. Ils se distinguent entre érosion 


éolienne et érosion hydrique. 
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1.3 formes de l’érosion hydrique 

1.3.1 L’érosion en nappe 

L’érosion en nappe appelée aussi érosion diffuse (“sheet erosion” pour Anglophone) 
résulte principalement de l’énergie des gouttes de pluies qui s’abattent sur la surface du sol. 
Par l’effet du ‘splash’, la pluie éboule les agrégats du sol et détache les fines particules qui 
se dispersent sur de courtes distances. Ces fines particules mobilisées sont ensuite 
transportées par le ruissellement en nappe (Figure. I.3) (Le Bissonnais, Y., & Le Souder, C. 


(1995). 


Ce type d’érosion dépend de deux processus : le détachement des particules du sol causé par 
l’impact de gouttes de pluie et le ruissellement superficiel (si l’intensité dépasse la vitesse 


d infiltration). 


Dans les sols cultivés, l’érosion en nappe peut être considérable et se traduit surtout par une 
« dilution » de la matiére organique, du fait du labour qui remonte en surface une nouvelle 
tranche de limon (Le Bissonnais et al 2000). Elle s’amplifie avec l’intensité des pluies, la 
raideur des pentes, le manque de couvert végétal et la vulnérabilité des sols (Bahri et al., 
1995). Ceci fait que l’érosion en nappe est plus intense en automne. En Algérie, l’érosion en 


nappe produit entre 0,10 a 20 t/ha/an de sédiment (Roose, 1993; Chibbani et al. 1999). 


Figure. I.3 : Un paysage qui illustre l’effet de l’érosion en nappe (photo : Guesri 
2015 a Bousaada 
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1.3.2 L’érosion concentrée (linéaire) 

Elle apparaît à l’intersection de deux plans de versant, et s’étend jusqu’au fond des 
vallées à travers le réseau hydrographique. Lorsque les eaux de ruissellement se concentrent, 
elles peuvent selon la nature du sol et l’intensité du relief creuser la terre en formant des 


griffes, des rigoles ou des ravines. 


On parle de griffes, lorsque les petits canaux ont quelques centimètres de profondeur, 
de rigoles, lorsque les canaux dépassent 10 cm de profondeur et de ravines lorsque les creux 


atteignent plusieurs dizaines de cm (plus de 50 cm) (Roose, 1994). 


1.3.2.1 L’érosion en griffes et rigoles 
L’érosion en rigoles commence lorsque l’écoulement de surface se concentre dans 
les points bas de la terre, l’énergie issue du débit concentré est capable de détacher les 
particules du sol, ce qui permet avec le temps de former des griffes et des rigoles (Figure 


14). 


Formation des rigoles 


Figure. I.4 : Le début de formation des rigoles a une terre située a la région de 


Maadid, M’sila (Photo Guesri 2017) 


1.3.2.2 L’érosion en ravines 


Elle est causée par les débits générés des ruissellements concentrés sur un versant (Figure 
I.5). Elles constituent des réseaux pour rejoindre le réseau hydrographique. Le ravinement 
est favorisé lorsque les sols sont sur les terrains imperméables soumis a des averses courtes 


et intenses. 
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En Algérie, l’érosion en ravines déplace ne moyenne entre 90 à 300 t/ha/an ce qui correspond 


à une vitesse d’altération des versants de 8 à 17 mm/an (Chebbani et al, 1999). 


Figure. LS : L’érosion en ravines est plus spectaculaire à l’Oued Rithem, Sud de 


M’sila. La profondeur de la ravine est d’une dizaine de mètres en comparaison avec la 


longueur des poteaux électrique (Photo Guesri 2017) 


1.4 Les mécanismes d'érosion, transport et dépôt de sédiments 

L’érosion hydrique des sols résulte de divers processus : le détachement des 
particules du sol puis leur transport par l’eau. Lorsque l’énergie nécessaire au transport des 
particules par l’eau n’est plus suffisante, un troisième processus intervient. Il s’agit du dépôt 
des particules, encore appelé sédimentation. Chacun de ces processus est causé par un 


ensemble de mécanismes. 


1.4.1 Les mécanismes de l'érosion 

Les mécanismes de l’érosion se différencient selon la nature des interactions en jeu 
et de l’énergie nécessaire à leur rupture, des conditions physico-chimiques nécessaires pour 
leur mise en œuvre, de la cinétique des processus ainsi que de la nature et la taille des 
particules résultant de la désagrégation. Plusieurs mécanismes ont été identifiés (Emerson, 
1967 ; Boiffin, 1984 ; Emerson et Greenland, 1990 ; Le Bissonnais, 1996): L’effet 
« Splash », l’éclatement, la désagrégation par gonflement différentiel, la dispersion physico- 


chimique et le ruissellement (Figure I.6) 
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par l'impact des gouttes d’eau 


W- Désagrégation mécanique 


Gonflement différentiel 
des argiles 


Eclatement par (sols argileux) 


compression de 
lair piégé 
(sols limoneux) 
= —— 
Dispersion physico-chimique 
des argiles 
(sols saturés Na, K) 


Figure. I.6 : Mécanismes de désagrégation à l’eau (Chenu et Cosentino, 2011) 


1.4.1.1 L’impact des gouttes de pluie « l’effet Splash » 

Ce type de mécanismes est dû au choc de la goutte de pluie sur le sol. Chaque goutte 
possède une énergie cinétique fonction de son poids et de sa vitesse de chute (Leguedois, 
2003 ; Kinnell, 2005). L’impact des gouttes de pluie peut fragmenter les agrégats et surtout 
détacher les particules de leur surface (Figure I.7). Les grains de sables, les fragments 
d'agrégats et les argiles dispersés sont projetés autour du lieu de l'impact. L'intensité de ce 


processus et la distance de projection dépendent de l'énergie cinétique des pluies, de la 


texture du matériau et de la pente (Le Bissonnais et al, 1996). 


Figure. I.7 : L’effet « splash » sur le sol sous l'impact d'une goutte de pluie 


(Photo: Pierre-Olivier Cochard Référence !) 
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L’énergie cinétique d’une goutte de pluie lors de son contact au sol est fonction de la masse 


et de la vitesse de chute de la goutte. Elle est donnée par l’équation : 


E, zipp 
2 


m : la masse de la goutte 


V : Ja vitesse de la goutte 


E, : Pénergie cinétique en joule (si ™ en Kg et V en m/s). 
Le diamètre des gouttes de pluie, leur vitesse de chute et leur masse totale déterminent la 


capacité de détachement des particules du sol. 


1.4.1.2 L’éclatement 
Il s’agit de la désagrégation par compression de l'air piégé lors de l'humectation. 
L'intensité de ce mécanisme dépend entre autres, du volume d'air piégé, donc de la teneur en 
eau initiale des agrégats et de leur porosité. Les sols argileux sont moins affectés par l'éclatement 
que les sols limoneux ou sableux (Le Bissonnais, 1990). L’humidification des sols avant une 
averse diminue considérablement le taux de détachement de particules (Roose, 1986). 
L’éclatement est qualifié comme un bon indicateur permettant de décrire la susceptibilité 


des sols à l’érosion (Barthès & Roose, 2001). 


1.4.1.3 La désagrégation par gonflement différentiel 
Ce phénomène intervient suite à l'humectation et la dessiccation des argiles, 
entraînant des fissurations dans les agrégats. L'importance de ce mécanisme dépend en 
grande partie de la teneur et de la nature de l’argile des sols. Ce mécanisme contribue à la 
réduction de la taille des agrégats à la surface du sol. La fissuration par gonflement augmente 
lorsque la teneur en argile augmente, du fait de l'augmentation du gonflement potentiel. Les 
conséquences directes de ce mécanisme sont moins importantes que les autres mécanismes 


(Le Bissonnais et Le Souder, 1995). 


1.4.1.4 La dispersion physico-chimique 
Ce mécanisme est le résultat de la réduction des forces d'attraction entre particules 
colloïdales lors de l'humectation (Emerson, 1967). L’eau s’infiltre par les fissures, dissous 
le sel contenu dans le sol et creuse des tunnels provoquant des effondrements qui se 
manifestent, parfois par d’importantes ravines qui marquent le paysage des versants. Ce 


mécanisme de désagrégation est le plus efficace, car il concerne la structure élémentaire de 
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l’agrégat. Il décuple les effets des autres mécanismes (Le Bissonnais et al, 1995). La 
dispersion physico-chimique concerne les particules élémentaires de l’agrégat, et dépend de 
de la taille et la valence des cations (particulièrement du sodium) pouvant lier les charges 
négatives dans la matière colloïdale. Ces dernières constituent des ponts entre les particules 
chargées négativement. La dispersion dépend alors, essentiellement, de la capacité d'échange 


cationique portée par le complexe agrégat. 


Les différents mécanismes de désagrégation et ses caractéristiques sont résumés dans le 


tableau suivant : 


Tableau I.1 : Caractéristiques des mécanismes de désagrégation (Le Bissonnais, 2000) 


Mécanisme Éclatement L'effet Désagrégation Dispersion 
« splash » par gonflement | physico- 
différentiel chimique 
Nature des Pression interne | Impact des Pression interne | Réduction des 
forces de l’air piégé gouttes de par gonflement | forces 
en jeu lors de pluie de l’argile d’attraction 
l’humectation entre 
particules 
Intensité de la Totale Cumulative Partielle Totale 
désagrégation 
Propriétés du Porosité, Cohésion Gonflement Statut ionique, 
sol en jeu mouillabilité, humide potentiel, minéralogie de 
cohésion (argile, cohésion, l’argile 
interne, matière conditions, 
mouillabilité, organique, d’humectation. 
oxydes) 
Particules Petits Particules et | Macro et | Particules 
produites microagrégats micro-agrégat | microa-grégats | élémentaires 
(< 500 um) 
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L.4.2 Les mécanismes du transport des sédiments 
1.4.2.1 Le ruissellement 

Suite à la phase de détachement des particules du sol les sédiments arrachés peuvent 
être mise en mouvement. L’eau qui ruisselle va les emporter et suivre le sens de la pente. 

Selon Fournier (1969), les agrégats du sol sont en général de dimensions trop grandes 
pour être facilement entraînés par le ruissellement. Ce déplacement des particules de sol est 
favorisé par plusieurs autres facteurs moins importants, tels que la gravité, le labour et le 
vent.Le transport concerne particulièrement les fines : comme l’argile et le limon, les 
matières organiques de faible masse volumique telles que les résidus de culture et les 
déjections animales ou encore les fertilisants épandus sur les parcelles cultivées. Cependant, 
les dépôts sont proches sous forme de colluvionnement. La vitesse du ruissellement sera 
d’autant plus rapide que la pente est longue, fortement inclinée et dépourvue de végétation. 
Plus la vitesse de l’eau est rapide et plus son énergie est importante, ce qui lui permet de 
transporter des particules plus grosses et plus loin. 

1.4.3 Les mécanismes de dépôt des sédiments 

L’eau ruisselée à la surface du sol est généralement chargée de sédiments mobilisés 
par les différents mécanismes d’érosion et de transport. Quand la pente devient moins 
inclinée ou que l’eau ruisselée rencontre un obstacle, son énergie diminue. A partir de ce 


moment, les particules transportées vont se déposer ou, en d’autres termes, sédimenter. 


La sédimentation constitue le dépôt en aval de particules ou sédiments arrachés au 
sol ou à d’autres éléments (roche) en amont. Ces éléments peuvent se déposer sur le sol ou 
au fond des cours d’eau. Quand on observe les sédiments, on remarque un classement des 
particules selon leur taille. En effet, les particules les plus grossières sont les premières à 
sédimenter, les plus fines étant transportées plus loin. On parle de tri granulométrique 


(Kauark, 1990). 


1.5 Facteurs de l’érosion hydrique 
1.5.1 Climat 

Le climat est considéré comme le facteur moteur qui influence l’érosion. L’efficacité 
de la désagrégation des sols est principalement liée à l’’intensité de la pluie et du vent ainsi 


que la hauteur de précipitation. 
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1.5.1.1 Précipitation 


L’agressivité des précipitations est considérée comme l’agent causal de l’érosion 
hydrique (Fournier 1969). Cette agressivité est appréciée par l’indice de l’érosivité R, utilisé 
dans les différentes versions de l'équation universelle de la perte des sols USLE, MUSLE, 
RUSLE, RUSLE2 (Modèles empiriques établis pour évaluer la perte en sol). L’indice R, 
quantifie l’énergie de l'impact des gouttes de pluie sur le sol (Wischmeier et Smith 1958). 
L’érosivité d’un évènement de pluie se calcule par le produit de l’énergie de la pluie (E) et 
de l’intensité maximale de la pluie pour une période de 30 minutes (Wischmeier et Smith 


1958 ; Wischmeier et Smith 1978 ; Renard et al 1997). 
1.5.1.2 Influence du vent 


Lorsque la vitesse du vent dépasse un seuil critique, elle donne suffisamment 
d’énergie cinétique pour faire bouger les particules du sol. Les particules commencent par 
rouler sur la surface et sur de courte distance. Si l’intensité augmente, les particules 
rebondissent dans les airs, où le mouvement du vent est plus rapide. L'effet de la gravité fait 
retomber ces particules à la surface où elles rebondissent ou entrent en collision avec d'autres 
particules. Ce processus est connu sous le nom de " saltation ". Par contre, les très petites 
particules sont entrainées dans l'air en ‘vol libre’ et sont transportées sur de longues distances 


(Figure I.8). 


Figure L8 : Le vent provoque la saltation, ou le saut des grains, du côté au vent des 
dunes de sable (Source : https://www.wikiwand.com/eu/Prozesu_eoliko ). 
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1.5.2 La topographie 

La longueur et l’inclinaison de la pente sont parmi les facteurs qui conditionnent 
considérablement l’érosion des sols. Une forte inclinaison de la pente agit directement sur 
la vitesse du ruissellement vers le bas et augmente l’énergie cinétique des eaux concentrées 


et par conséquent, le potentiel érosif. 


Fox (2008) trouve que la longueur de pente est moins importante que l’inclinaison. Au-delà 
d’une centaine de mètres, la longueur a peu d’influence sur la vitesse d'écoulement 
maximale. De plus, il a souvent été remarqué que le ruissellement diminue le long d’un 
versant de l’amont vers l’aval, peut-être dû à des zones d’infiltration plus élevées le long de 


la trajectoire de l’écoulement de surface. 


1.5.3 L’état du couvert végétal 

La nature et l’extension de la couverture végétale jouent un rôle important dans l’interception 
et la formation du ruissellement superficiel. Le feuillage des plantes permet de dissiper 
l’énergie cinétique des gouttes de pluie et par conséquent la diminution de l’effet “splash”. 
Par ailleurs, par leurs racines, les plantes favorisent l’infiltration de l’eau vers le sous-sol et 


retardent le ruissellement superficiel. 


1.5.3.1 Le défrichement 

La perte des arbres, qui ancrent le sol avec leurs racines, provoque une érosion généralisée 
dans les zones semi-arides. Seule une minorité de zones ont de bons sols qui, après avoir été 
défrichés, sont rapidement emportés par les fortes pluies. Ainsi, les rendements des cultures 
diminuent et les gens doivent dépenser leurs revenus pour importer des engrais étrangers ou 
défricher d'autres forêts 

En générale, sur les versants dénudés (défrichement,) le ruissellement se génère rapidement 
et intensément, s’organise de manière diversifiée en fonction de la pente, de la nature des 
matériaux et de la rugosité de surface et attaque le sol, de plus en plus en profondeur, et 
transporte les sédiments jusqu’à l’affleurement de la roche en place. Le processus continue 


jusqu’à une généralisation du ravinement et la formation de badlands. 


1.5.3.2 Les incendies 

Le feu peut influer directement sur les propriétés du sol par des processus de chauffage et de 
combustion (Santin et Doerr, 2016). Plusieurs phénomènes sont mis en cause comme la 
disparition de l’humus, la formation d’une croûte liée à la cuisson du sol, le dépôt d’une 
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couche hydrophobe de cendres qui diminuent la perméabilité du sol (Gottfried et al, 2003 ; 
Neary, 2004). Les feux influencent indirectement l’érosion par des modifications de son 
couvert végétal et une meilleure redistribution du sol grâce à une érosion accélérée après 
l'incendie. 
De nombreuses études dans la littérature ont fait état d'un degré variable d'érosion du sol 
après les feux de forêt. 

1.5.4 L’état structural du sol 

La stabilité structurale correspond à la capacité d'un sol à conserver son arrangement 
entre particules solides et vides lorsqu'il est exposé à différentes contraintes, comme l’action 
désagrégeante de gouttes de pluie ou l'humectation (Le Bissonnais, 1996). On la considère 
comme un bon indicateur de la susceptibilité des sols à la battance et à l’érosion hydrique 


(Le Bossonnais et al, 2002). 


L’état de la stabilité des agrégats permet d’évaluer la résistance de ces derniers aux 
forces de désintégration. Plus la stabilité des agrégats est faible dans l'eau, plus la 
susceptibilité du sol à la détérioration de sa structure par l'impact des gouttes d'eau ou par 
l'action des outils de préparation culturale du sol, ou son tassement par le passage des 


machines agricoles, sera grand. 


La matiére organique du sol (M.O.S.) est un des facteurs importants de l'agrégation 
des particules solides dans le sol. Dans les sols non calcaires ou qui contiennent des teneurs 
faibles d'oxydes de fer et aluminium, elle est méme le facteur principal de cohésion des 


agrégats entre eux (Oades, 1993 ; Angers et Carter, 1996). 


La matiére organique influence la stabilité de la structure du sol par plusieurs 
mécanismes, son effet positif sur la stabilité structure est variable en fonction de la nature, 
de la quantité et des conditions du milieu (Aoyama et al. 1999 ; Spacinni et al., 2004). 
Certains travaux montrent une liaison claire entre la teneur en carbone total et la stabilité 


structurale (Dutartre et al., 1993 ; Le Bissonnais et al., 2002 ; Chenu et al., 2011). 


1.5.4.1 La désertification 


La désertification est une forme de dégradation des terres fertiles dans lequel une zone 
relativement sèche devient de plus en plus aride, perdant généralement ses étendues d'eau 


ainsi que la végétation. 
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Une étude cartographique réalisée par Benslimane el al (2009) a permis d'évaluer la 
sensibilité de ces terres à la désertification (Figure I.9). Ses résultats montrent que 74,81 % 
des parcours steppiques sont au seuil de la désertification. L’évolution phénoménale de ce 
fléau naturel, durant ces 20 dernières années s’accélère selon un rythme d’environ 40 000 


ha/an. 


M couvent végétal très dense 

DM couvert végétal dense 

Mn état critique : zone sensible à la dégradation 
M zone dégradée 

HMM zone très dégradée 
M pans deau 


Figure I.9 Carte synthèse de sensibilité à la désertification entre 2000-2005 (Benslimane et 
al, 2009) 


1.5.5 Le ruissellement 
Quand le sol ne peut plus absorber l'eau des précipitations, une lame d'eau se forme à la 
surface du sol, correspondant à l'eau en excès. La mise en mouvement de cette eau forme le 
ruissellement. Ce dernier peut apparaitre dès que l’intensité d’une pluie devient supérieure à 
la vitesse d’infiltration du sol. Il est nettement plus élevé sur sol nu que sous cultures. 
Les plus faibles ruissellements proviennent des pluies tombantes après une période sèche. 


(Bouanani, 2004). 


Chebbani et al (1999) trouvent que les ruissellements et les pertes en terres les plus élevées 
ont eu lieu à une époque où le sol est déjà très humide et dans le cas de pluies exceptionnelles 


et des orages d’automne et de printemps. 


D’après Laouina (1998), le ruissellement ne débute qu’aprés un cumul pluviométrique 
journalier de 18 à 20 mm sur sol sec même si les pluies sont intenses. Alors qu’il apparait 


après un cumul de 4 mm sur sol humide et compact (Chebbani et al. 1999). 
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Les pluies continues de faibles intensités permettent une certaine infiltration de l’eau et 
réduisent le ruissellement sur les versants, par contre les épisodes de fortes intensités 
entraînent un ruissellement plus important. (Laouina al, 1993) 

Mais ce sont les épisodes pluvieux qui durent dans le temps avec des événements répétitifs 
de forte intensité qui occasionnent la dégradation la plus forte, et où le ruissellement en 
nappe passe au ruissellement en griffes et rigoles. 

En général, la vitesse du ruissellement est le paramètre prépondérant de l’action érosive du 


ruissellement superficiel. (Bouanani, 2004). 


1.5.6 Action humaine 
Le phénomène de l’érosion n’est pas conditionné seulement par des paramètres et 


des processus naturels, mais aussi par les actions anthropiques comme : 


- L’intensification de l'agriculture, l'élevage des animaux, les exploitations minières, 
l'abattage du bois, la construction de routes, les aménagements urbains et le 


surpâturage ; 


Ces différentes activités humaines ont une grande influence sur la terre, car elles peuvent 
avoir un effet important à cause de la destruction de la stabilité structurale du sol, 


la diminution ou la suppression de la végétation et modifications topographiques. 


Le déboisement illégal et les feux de forêts exposent le sol nu à l’agressivité des pluies, ce 
qui permet de faciliter le détachement des particules, et de favoriser le ruissellement et par 
conséquent causer une forte érosion (Benchaabane, 1998). Le surpâturage a deux actions sur 
l’augmentation des phénomènes érosifs, le premier est la dégradation du couvert végétal 
protecteur, le second est un effet mécanique sur la structure des sols en pente par les chemins 
suivis par les troupeaux. D’après (Dajoz, 1996), Les actions et les fautes de l’homme ont 


causé une perte de 5 millions km? de terres cultivables au niveau mondial. 
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1.6 Le réseau hydrographique 

Le réseau hydrographique se développe de manière à évacuer efficacement l'eau de la terre. 
L'écoulement des cours d'eau commence sous la forme d'une mince couche d'eau en surface. 
L'eau descend la pente la plus abrupte et commence à éroder la surface en créant de petits 
canaux de ruisseau. Au fur et à mesure que les rigoles fusionnent, s'approfondissent et 


s'enfoncent dans les canaux, de plus grands canaux se forment. 


L'érosion rapide allonge la pente ascendante du canal dans un processus appelé érosion vers 
l'amont. Au fil du temps, les chenaux voisins fusionnent avec des affluents plus petits qui se 
joignent à un ruisseau principal plus important. Les chenaux reliés deviennent ce que l'on 
appelle un réseau de drainage. Avec l'érosion continue des chenaux, les réseaux de drainage 


changent au fil du temps. 
I.7 Transport solide 


Le lit des cours d’eau naturels est constamment formé et déformé par les sédiments en 
mouvement. Dans les cours d’eau alluviaux, les mouvements des sédiments sont variables 
dans le temps et dans l’espace, et s’adaptent aux changements de l’écoulement, modifiant la 
morphologie du lit. Ainsi, la compréhension de la dynamique sédimentaire des cours d’eau 
nécessite la connaissance des interactions qui lient l’écoulement, le transport des sédiments 


et le développement des formes sur le lit. 
1.7.1 Modes de transport des particules solides 
On distingue deux types de transport (Figure I.10) : 


1.7.1.1 Le transport solide en suspension 

Il concerne les particules les plus fines et s'analyse comme un équilibre « entre deux 
eaux », entre le poids déjaugé de la particule (poids réduit de la poussée d'Archimède), les 
forces s'opposant au mouvement relatif de l’eau et de la particule dues a la viscosité de l’eau 
et les composantes verticales des forces turbulentes. En l’absence de turbulence s'établit une 
vitesse de chute qui résulte de l'équilibre entre poids déjaugé et forces de friction dues à la 
viscosité (Roche et al, 2012). 

1.7.1.2 Le transport solide par charriage 
Il est lié au mouvement possible de matériaux posés au fond du lit (gravier, sables, galets et 


blocs), de diamètre supérieur à 1 mm, qui se déplacent par roulement ou glissement, au 
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voisinage immédiat du fond du cours d’eau, et parcourent des distances relativement petites 
et discontinues dans le temps et dans l’espace. 

Les mesures sur le terrain du transport solide par charriage sont très onéreuses à cause des 
perturbations de l’écoulement. De ce fait, les évaluations de ce type de transport sont 


difficiles à mener en l’absence de bases de données. 


Suspention 


Saltation 


3 AE Roulement 


Charriage 


Figure 1.10 Modes de transport des particules dans un chenal (Source : http://theses.univ- 
lvon?2.fr/documents/getpart.php?id=lyon2.2008.pintomartins d&part-154405 ) 
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Chapitre II : Présentation de la zone d’étude 


IL.1. Aperçu général sur le bassin Hodna 

D’après Savornin (1908), le nom de Hodna (en arabe : 4&3), semble à répondre, 
dans la toponymie indigène, à quelque chose d’assez précis. Autres spécialistes de la langue 
y voient un terme générique et le traduisent : « Plaine entourée de montagnes », cette 
définition convient parfaitement à la région connue sous ce nom en Algérie. Le bassin 
versant du Hodna de surface de 25840 km? et de périmètre de 1059 km, est un bassin 


endoréique où ses affluents déversent dans le lac salé Chott El Hodna. 


Légende 


‘ Grands bassins versants de l'Algérie 
120 240 Kilometers) —— 


h a Bassin versant Hodna 
Code Nom du bassin versant Code Nom du bassin versant 
1 Chéliff 9 Isser 
2 Cotiers Algérois (Est-Ouest) 10 Kébir Rhummel 
3 Cotiers constantinois (Est-Centre-Ouest) 11 Macta 
4 Cotiers oranais (Est-Centre-Ouest) 12 Medjerdah 
5 El Hodna 14 Seybouse 
6 Chott Melrhir 15 Soummam 
7 Hauts plateaux constantinois 16 Tafna 
8 Hauts plateaux oranais 17 Zahrez 


Figure IL.1 Localisation du grand bassin el Hodna en Algérie du Nord 
(Database: World Imagery Basemap, ESRI) 
Il est en position centrale dans le territoire du Nord de l’Algérie (Figure IL.1) et en 


contact de ce fait avec de nombreux ensembles hydrographiques : 


o Au Nord, par les bassins versants de l’Isser et de la Soummam, 


o À Ouest, par les bassins du Cheliff et du Zahrez, 
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o À l'Est, par les Haut Plateaux Constantinois, 


o Au Sud, par le grand bassin du Chott Melrhir. 


Légende 


=] Sous bassins versants 
——— Réseau hydrographique 


a Chott el Hodna 


60 Kilometers 
Lu er Von 1 1 |) 


Figure II.2 bassin el Hodna selon le découpage de ANRH 
Selon l’agence nationale des ressources hydrauliques « ANRH », le bassin de Hodna 
peut être divisé en 23 sous-bassins endoréiques : Oued Ham (805,8 km’); Oued Sbisseb 
(1851,98 km?); Oued Guernini (1071,31 km?); Oued Djenane (860,47 km?); Oued Targa 
(1026,94 km?); Oued El Karsa (553,12 km?*); Oued Leham (647,99 km?) ; Oued Lougmane 
(333,75 km?) ; Oued K'sob (1480 km?) ; Oued M'Sila (2157,52 km?) ; Oued Soubella 
(1768,97 km?) ; Oued Barika Amont (889,77 km?) ; Oued Barika Aval (1007,05 km?) ; 
Oued Berrich (504,4 km?) ; Oued Bitam-Barika (1403,02 km?) ; Oued Maiter Amont 
(672,74 km?) ; Oued Bou Saâda (1035,15 km?) ; Oued Maiter Aval (1232,86 km?) ; Oued 
Ain Rich (1130,88 km?) ; Oued Chair Amont (1597,78 km?) ; Oued Chair Aval (877,75 
km?) ; Oued M'Sif Amont (924,2 km?) ; Oued M'Sif Aval (792,66 km?). Le 24ème central 
étant le Chott El Hodna (Figure II.2). 
Le chott el Hodna est un lac salé endoréique du nord-est de l'Algérie (Figure I.3). 
Il est le deuxième plus grand lac d'Algérie (après le chott Melrhir) et donc l'un des plus 
grands d'Afrique du Nord. Il fait partie d'une série de chotts qui se sont développés aux 
endroits de convergence des eaux provenant de l'Atlas saharien au sud (Aurès) et de l'Atlas 


tellien au Nord (Chréa et Djurdjura). 
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Figure II.3 Chott el Hodna (Photo : Amal Qahwadji, 2011) 
La cuvette (bassin) du Hodna a une superficie de 8500 km? dont 1100 km? pour le 


chott à proprement parler. Le chott est situé à 400 mètres d'altitude. Son alimentation est 
assurée par au moins 22 cours d'eau principaux auxquels il faut ajouter des sources d'eau 
douce. Il se met en eau uniquement en hiver, sec et salé en été, des croûtes de sel couvrent 
toute son étendue. À cause d'une forte évaporation, il se transforme régulièrement en un 


désert de sel. 


*Hammam|Dhalaa 
Ras/El\Qued 


Berhoum Ras|El/Aioun 
Magra 


DfAin|Khadra M Belaiba NGaous 


23 Ouled}Derradi 


: Sefiane 
DjezzarBarika 


iBen}Srour, 
Djebel Messaad à 
Légende 
AinlEl|Melh 

e Ville 


Chott el Hodna 


FaidhjEl' Botma 


60 Kilometers 


Figure IL.4 Vue Satellitaire du bassin versant Hodna (Database: WIB, ESRI) 
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Au nord du bassin, l'écran montagneux formé par la chaîne du Hodna et des monts du Titteri 
forme une véritable barrière. 

Au pied de cet obstacle, la plaine du bassin-versant est l'unité la plus vaste. Entre les massifs 
montagneux et la plaine centrale, le piedmont est constitué par de vastes glacis d'érosion 
parsemés de chebkas (ravins) (Kebiche, 1994). 

IL.1.1 La dépression Hodnéenne 


La dépression de El Hodna est une plaine où s’est creusé un réseau hydrographique dense 


qui raine les versants du cadre montagneux environnant (Figure II.5). 


Légende 


Altitude (m) 
Valu 


e 
- High : 2160 


- Low :0 


e Ville 


[F] Chott el Hodna 


60 Kilometers 


Liilia 


Figure II.5 Les grands ensembles du relief du bassin versant Hodna (Carte basée sur les 


Modèles Numériques du Terrain « ASTER GDEM », avec une résolution fine de 30 m) 


Cette dépression est subdivisée en deux milieux bien distincts : 

*Le piedmont (djerr) : 
Le piedmont ‘Djerr’ est la zone de transition entre la montagne au nord et la plaine. Il 
s’allonge d’Ouest en Est a entre les cotes 700m et 500m. C’est un milieu de collines aux 
formes très variées, certaines sont disséquées par le réseau hydrographique, ce qui leur donne 
— vu de loin — la forme de filets d’où l’appellation locale de chebka. 
Ces formes très développées à |’ Est de l’oued K’sob, constituent des dépôts continentaux du 


mio-pliocène qui forment des petits reliefs monoclinaux. 
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En plus de ces collines, ce milieu se caractérise par une autre forme : les ’draàs’ qui signifient 
bras ; c’est une série de glacis d’érosion étagés qui convergent vers l’aval et qui sont 
découpés en lanières par des vallons. Ces formes se développent sur toute la partie 
occidentale, et entre les collines de la partie orientale. 
Sur le versant Sud, la zone de piedmont est inexistante et les versants passent directement à 
la plaine sans aucune zone de transition ; le contact étant masqué par les formations 
éoliennes. (Nemouchi, 2001) 
La plaine s’étend au Sud et à l’Ouest du piedmont. Elle est circonscrite par la courbe de 
niveau de 500m, et s’étale jusqu’aux rives de la Sebkha. C’est une plaine de remblaiement 
qui couvre environ le1/3 du bassin versant mais avec des nuances entre la plaine du Nord et 
de l’Est d’un côté et celle du Sud de l’autre. La première, façonnée par les eaux de 
ruissellement issues du cadre montagneux est une plaine limoneuse parsemée de nappes de 
cailloux et de sable résultants du dépôt de la charge solide des oueds qui perdent leur 
compétence avant d’atteindre le chott. Des deltas se sont constitués à chaque embouchure 
d’Oued et gagnent en extension sur le chott. La seconde, se caractérise par la nature 
sablonneuse de ses sols, d’où son appellation locale de r’mel qui signifie sable. Elle s’étend 
directement aux pieds des reliefs qu’elle entoure par des étendues planes ou  parsemées 
d’accumulations de sables (erg) dont les plus importantes sont localisées à l’Ouest et au Nord 
de Boussaada. 
Contrairement à la plaine du Nord et de l’Est, le ruissellement n’a aucun effet important, car 
les oueds n’atteignent qu’exceptionnellement le chott — sauf les oueds M’ cif et Boussaada — 
ils forment de larges plaines d’accumulation ou les eaux s’infiltrent et s’évaporent. Par 
contre l’action éolienne est prépondérante, le vent est le principal agent d’érosion, de 
transport et d’accumulations de sable (Nemouchi, 2001) 
IL.1.2 Le couvert végétal 

La caractéristique essentielle du bassin versant du Hodna du point de vue couverture 
végétale est la prédominance des terrains qui favorisent le ruissellement notamment les 
terrains de parcours et les terrains nus dont la proportion est > 65% de la surface totale du 


bassin. Le reste est occupé par terres de culture et les forêts et matorrals 
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4 Ain El Hadjel Æ 


0 15 30 60 Kilometers 
Loria Los nr i | 
Légende 
E] Irrigated croplands EM Mosaic Forest-Shrubland/Grassland 
CZ] Rainfed croplands EM Mosaic Grassland/Forest-Shrubland 
EI Mosaic Croplands/Vegetation EM Closed to open shrubland 
EMI mosaic Vegetation/Croplands EM closed to open grassland 
O i Closed to open broadleaved evergreen or semi-deciduous forest =] Sparse vegetation 
MB Closed broadleaved deciduous forest EM Closed to open broadleaved forest regularly flooded 
==) Open broadleaved deciduous forest EM Closed broadleaved forest permanently flooded 
= Closed needleleaved evergreen forest = Closed to open vegetation regularly flooded 
MM Open needleleaved deciduous or evergreen forest GE Artificial areas 


EM Closed to open mixed broadleaved and needleleaved forest L] Bare areas 


Figure II.6 Carte de l’occupation de sol du bassin Hodna (Carte réalisée en se basant sur la 


carte de l'occupation des sols GlobCover2009, couvrant l'Afrique et la péninsule arabique) 


Les terrains boisés ont une extension faible, ils sont localisés généralement en altitude 


et sur les pentes raides ou les précipitations sont abondantes et l’action anthropique limitée. 
Les espèces dominantes sont le cèdre sur les sommets excédants 1500m, le pin d’Alep, le 
chêne vert, et l’olivier lentisque aux altitudes inférieures ou l’humidité s’atténue légèrement. 
Ces forêts sont dans un état de dégradation avancée, et la végétation climax n’a pu se 
maintenir qu’en quelques ilots sur les sommets de montagnes (Boutaleb, Maadid, ... etc) ou 
les conditions de sa conversation (humidité, absence de l’action anthropique) sont réunis. 
Ailleurs, cette forêt est dégradée, elle se présente sous forme de couverture arbustive formée 
par un maquis épineux. 


Les terrains mal protégés ou nus représentés par les parcours et les terrains dénudés 


occupent la moitié de la surface totale à cause du caractère steppique de la région. En effet, 
l’isohyète 400mm constitue la limite écologique de la steppe ou elle occupe aussi bien les 
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piedmonts, les vallons et les dépressions. C’est une formation basse, discontinue, privée 
d’arbres et d’arbustes. Les espèces dominantes sont le sparte, le Chih, le drinn, et l’armoise 
qui se répartissent selon les conditions édaphiques du milieu. 

La dégradation de ces terrains par la surcharge pastorale et certaines pratiques qui y sont 
liées (le feu engendre la disparition de cette végétation et l’apparition de terrains 
complètement nus qui s’étendent sur des surfaces relativement importantes. 


Les terres de culture occupent 1/3 de la surface totale, avec la prédominance des 


céréales et de la jachère qui s’étendent aussi bien sur les plaines, les plateaux, et les terrains 
en pente. L’arboriculture et les maraichages ont une faible extension, limitée le long des 
oueds et autour des sources (Nemouchi 2001). 

IL.1.3 L’aspect hydrogéologique 
Le bassin du Hodna se présente en deux grands ensembles de réservoirs d'eau : la nappe de 
remplissage (Mio-pliocene) en profondeur et l'aquifère superficiel, constitues par les 


formations détritiques du quaternaire (Figure II.7). 


Légende 
MM Cenozoic marine sedimentary 

[=] Cretaceous sedimentary - North 

[E] Recent-Quaternary unconsolidated sediments 
[ __| Triassic to Jurassic - North; sedimentary 


0 15 30 60 Kilometers 
Lote a bp ‘je | 


Figure II.7 Carte lithologique du Bassin Hodna (Database: UNEP 2010) 
Les principales nappes aquiféres du bassin du Hodna correspond au : 
- Nappe phréatique qui se trouve à une profondeur supérieure à 40 mètres dans la partie 


haute des plaines du Hodna ou son exploitation est la plus rentable (taux de salinité 
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de l'eau est moins de 2 g/l). Près de la Sebkha, elle est moins de 03 mètres de la 
surface, mais son utilisation est déconseillée pour l'irrigation (le taux de salinité est 
jusqu’à 250 g/l de résidu). 

- La nappe captive : se situe à une profondeur entre 150 et 250 m. Elle est exploitée 
depuis 1859. Plus de 80 forages ont été réalisé dans le but d'obtenir de l'eau artésienne 
sans frais de pompage. Elle est peu chargée en sels, le résidu sec varie de 0.5 à 02 g/l 


(Mimoune, 1995). 


IL.1.4 La géologie 
La géologie des terrains influe non seulement sur l’écoulement de l’eau souterraine, mais 
également sur le ruissellement des surfaces à travers l’étude de la perméabilité des 
formations lithologiques. Celle-ci intervient sur la vitesse de montée des crues, sur leur 
volume, ainsi que sur l’aggravation du phénomène d’inondation et par conséquent l’érosion 
et le transport des sédiments. Les caractères géologiques principaux à considérer sont la 


lithologie et la structure tectonique (Bebba 2018). 


Légende 

Es Cretaceous (undivided) ES Lower Cretaceous 
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Figure II.8 Carte géologique du Bassin versant Hodna 
(Database : USGS Geology of Africa) 


IL.1.4.1 Stratigraphie et aperçu lithologique global 
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Les terrains les plus anciens qui affleurent dans cette région sont attribués au crétacé le plus 


récent et au quaternaire le plus ancien. 


Le Primaire 
La région est dépourvue de formations géologiques paléozoïques. Toutefois, au Nord de 
celle-ci ont été observés des niveaux de formations allochtones qui comportent des copeaux 


provenant du socle primaire (roches sédimentaires, éruptives et métamorphiques). 


Le Secondaire 

Trias s’agit d’une formation classique constituée essentiellement d’argile gypsifères grises 
à minces niveaux bariolés, à blocs de dolomies empruntées de l’infra Lias et du Crétacé 
inférieur. Les affleurements triasiques sont soit pseudo diapiriques (Bou Sâada), soit 
jalonnant les contacts anormaux (Ouled Tebbane, Boutaleb,...). On cite les types qui se 


trouvent au niveau du Hodna : 


- Le Crétacé 


Dans le grand bassin versant du Hodna, le Crétacé affleure largement dans la région et on 


note la présence d’un Crétacé inférieur et d’un Crétacé supérieur. 


- Le Crétacé inferieur 


Il est essentiellement carbonaté terrigène mais on note le développement de formations 
gréseuses avec d’importantes variations de faciès. Il n’affleure qu’aux noyaux et sur les 


flancs des anticlinaux au Nord et au Sud des Monts du Hodna. (Nemouchi, 2001) 


- Le Crétacé supérieur 


Les formations du Crétacé supérieur correspondent à un milieu néritique, tous les étages sont 
représentés et l’épaisseur des formations est importante, il est essentiellement marneux et 
calcaro-dolomique, ainsi que des dépôts détritiques grossiers (djebel Bellzma et djebel 


Metlili) 
-  Jurassique 


Le Jurassique n’affleure à la faveur d’accidents majeurs que dans la partie centrale et 
orientale des monts du Hodna. 
- Le Sinémurien se résume à des calcaires dolomitiques saccharoïdes apparaissant 


localement au Djebel Soubella à la faveur de nombreuses failles. Le Pliensbachien Il s’agit 
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d’une quinzaine de mètres de calcaires olithiques et calcaires gris, surmontés de 40 m de 
calcaires sublithographiques. 

— Toarcien moyen Le Toarcien se limite à une mince série marneuse de couleur 
essentiellement grisâtre, par endroits rougeatres. -Dogger — callovo - oxfordien s’agit de 
bancs de calcaires à silex et d’une alternance de marnes et de calcaires à silex. Cet ensemble 
est surmonté de 25 à 30 m de calcaires gris-bleutés, siliceux, terminés par des dalles de 
calcaires gris-bleutés et des calcaires à filaments. Oxfordien supérieur se résume à une 
quinzaine de mètres de calcaires marneux, plus ou moins violacés ou verdâtres. 

Le Kimméridgien est caractérise par 300 à 350 m de marne gris-bleues à rares passées 
calcaires vers la base, ainsi que de calcaires compacts et des calcaires marneux alternant avec 
des marnes grises. Le Tithonique est formé de calcaires gris alternant avec des marnes grises 
sur une épaisseur d’environs 250 m. Les calcaires sublithographiques du Tithonique sont les 


derniers traversés par le forage. 


Le Tertiaire 
Le Tertiaire continental repose en discordance sur les formations crétacées, résultant du 
démantèlement du relief de l’Atlas saharien par le phénomène de l’érosion. 
Ces formations continentales du Tertiaire, se rencontrent à l’Est et au Nord de Boussaâda, 
ainsi qu’à l’Ouest du Chott el Hodna. Cet ensemble renferme des conglomérats, des galets, 
des sables, des grès et des passées marneuses. Le Tertiaire marin n’est présent que dans la 
partie orientale de l’Atlas saharien, et comprend : 
vV Le Paléocène: Qui est formé de calcaires argileux, de lumachelles, de marnes à gypses 
et de souffre noir (une série de 40 à 60 m de marnes noire). 
Y L’Eocène: Qui est constitué de calcaire gris avec des intercalations de marnes grises, 
de calcaire à silex et de lumachelles. 
L’éocène inférieur, d’une épaisseur varie entre 150 à 200m, il est représenté par des 
marno- calcaire, marnes, calcaires blancs riches en silices et souvent fissurés, où il 
est représenté par quelques dizaines de mètres de marnes jaunes, verdâtre légèrement 
gypseuses et des calcaires argileux 
L’éocéne moyen est surmonté par une assise rouge gréso-marneuse à la base 
conglomératique au sommet, dont l'épaisseur atteint au maximum 300m. 
(Belagoune, 2012) 
v Le Miocène: Les formations marines du Miocène sont transgressives dans le Hodna. 


Plusieurs cycles de transgressions y sont dénombrés. 
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On distingue de bas en haut plusieurs faciès de miocène allochtone. En haut, les 
marnes rouges, des conglomérats avec la formation des grés fins beiges très clairs, à 
la base des bancs lenticulaires de gypses massif alternant par la formation des grés 
fins très clairs, beiges très massifs d’une épaisseur de 300m. 

Grés grossier, roux à patine foncée, glauconieux à la basse 


Grés fins, jaunes ou beiges alternant avec des marnes, 200 m d’épaisseur 
Le Quaternaire 


vV Le Quaternaire ancien: Ce niveau comporte des conglomérats à gros blocs, des 
horizons à galets et feuillets d’encroûtement calcaires rouge orange, situé entre 
20-40 mètres au dessus du niveau actuel des oueds. 

V Le Quaternaire moyen: Il est représenté par de gros bancs de calcaires, de galets 
emballés dans une matrice de croûte calcaire gréseuse, des conglomérats et d’argiles 
à débris prismatique. 

Ÿ Le Quaternaire récent : Au niveau du Chott el Hodna, il est constitué par des limons, 
de couleur grise, avec des cristaux de gypse. Par contre, dans les dunes on trouve des 


sables fins, des argiles et du gypse. Les galets sont localisés au pied des reliefs. 


33 


Chapitre II Présentation de la zone d'étude 


Figure II.9 Stratigraphie synthetique du Hodna (Tarmoul 2018) 


Tectonique 

La tectonique procède en premier lieu d’une analyse de la géométrie des couches, elle 
identifie les formes, les objets tectonique (les déformations cassantes (les failles), les 
déformations souples (les plis), détermine leur position dans l’espace puis par des mesures 
de leur orientation tente de restituer les contraintes responsables des déformations 


(Guillemot, 1986). 


Du fait de sa situation entre l’ Atlas Tellien et l’Atlas Saharien, le bassin du Hodna présente 


une grande complexité tectonique caractérisée par une grande variation de faciès (c’est 
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l’ensemble des caractères lithologiques et paléontologiques qui définit un dépôt et donne des 
indicateurs sur les conditions de formation et de dépôt des roches (Guillemot, 1986). 

Les phases tectoniques affectant le bassin du Hodna sont: Emscherienne (sononier 
inférieur), Atlasique (éocène), Tellienne (miocène), Messénienne et Villafranchienne 


(Seddi, 2012). 
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Figure IL.10 : Coupes géologiques du bassin du Hodna (Source : ENAGEO, 
http://www.alnaft.cov.dz/IMG/pdf/4-Bassin du Hodna Engeo) 
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IL.1.5 Les caractéristiques lithologiques et la perméabilité 
IL.1.5.1 Les caractéristiques lithologiques 


Les caractéristiques hydrologiques des sous bassins sont dans une large mesure influencée 
par l’imperméabilité des roches et les sols qui en dérivent d’absorber l’eau et de la retenir de 
manière temporaire. D'une manière générale de l'étude géologique, nous ressortons les 


différentes formations lithologiques. 


IL.1.5.2 Perméabilité des sols des bassins versants 
En général, dans les sous bassins de notre étude on constate trois types de roches ou trois 
grands ensembles. 

vV Les formations imperméables : Constituées par des marnes, marno-calcaires et sols 
argileux de telles assises rencontrent dans la formation marno-calcaire turonienne et 
la formation marneuse miocène, ainsi que sous forme de divers recouvrements 
quaternaire (sols forestiers). Sur les zones de cette nature, l’infiltration doit être 
inférieure à 10% des précipitations, s’élève à 23% et 12% pour les sols à texture 
moyenne et lourde or ces derniers, plus hétérogénes que les terrains considérés et 
plus perméables qu’eux. 

v Les formations à perméabilité moyenne : Ces formations sont formées de calcaires 
du cénomanien supérieur et les formations dolomitiques, ainsi que les marnes du 
miocène. Il s’agit des calcaires bio détritiques et su lithographiques du crétacé, ou 
des calcaires d’algues du miocène, ou encore de la dolomie grenue, ces roches sont 
très compactes ; la formation dolomitique retenue cependant de rares poches de 
minéralisation ferrugineuse plus au moins caverneuse, seules susceptibles d’accroître 
localement sa perméabilité. 

vV Les formations à forte perméabilité : Les conglomérâtes et les brèches du miocène 
marin, dont les éléments généralement très grossiers sont soudés entre eux par un 
ciment calcaire bio détritiques fin, laissant subsister autour d’eux des rives 
importantes, ces formations occupent Ain Hadjel. 

Les grés calcareux très vacuolaires, qui correspondent aux marno-calcaires miocène 
et les divers recouvrements quaternaires, très peu consolidés formés de débris 
calcaires ou dolomitiques à taille variable, parfois arrondis, de limon sableux et 


caillouteux verts (Bendjamma, 2000). 
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II.2 Présentation du sous bassin K’sob 

Selon le découpage de l’agence nationale des ressources hydrauliques ANRH, le bassin 
versant d’ étude, Oued K’ sob, est le sous bassin N°09 qui fait partie du grand bassin du Hodna 
(N°05) et se situe dans sa partie septentrionale. Le bassin de Oued K’sob est limité au Nord 
et au Nord-Ouest par la chaine montagneuse des Bibans, au Sud et au Sud-Ouest par les 
monts du Hodna et à l’Est par les hautes plaines de Sétif (Figure IL.11). 

Il se situe dans une zone entre les longitudes (4°30’-5°17’ Est) et les latitudes (35°45’- 
36°9’ Nord). Le bassin Oued K’sob est constitué de deux ensembles physiques distincts, il 
s’agit des bordures montagneuses, entre lesquelles s’étend une vaste plaine, appelée hautes 
plaines de Bordj-Bou-Arreridj. Les bordures sont formées au sud par les monts de Hodna 
d’orientation Est-Ouest dont les altitudes sont relativement élevées (Djebel Maadid 1819 m, 
Djebel el Mzaita 1545 m et Djebel Kender 1640 m) avec un point culminant de l’ordre 
1875m à Djebel Chelendj ; cette bordure est découpée en deux compartiments par de l’oued 
Ksob et au nord par la chaine de Bibane (Djebel Moutène) d’orientation SW-NE. Les 
bordures montagneuses se caractérisent par une structure complexe et charriée ; très écaillée, 
formées dans des formations secondaires très variées ; comme les calcaires, les grés les 
schistes, les marnes, les marno-calcaires, les dolomies, et les formations triasiques. 

L’étude ou l'analyse morphométrique du bassin versant sont présentés dans la partie 
suivante. La méthodologie de calcul inclut les caractéristiques importantes du bassin versant, 
la superficie, le périmètre, le coefficient de forme, l'ordre des cours d'eau, la densité du 
drainage, la pente moyenne du bassin versant, la pente du lit principal, etc. 

Les différents paramètres sont mesurés à l’aide du logiciel ArcGis, utilisé dans les systèmes 
d'informations géographiques, et aussi à l’aide des modèles numériques de terrain MNT 
(avec une résolution de 30 mètres) établis par l’Institut d'études géologiques des États-Unis 
(USGS). 

11.2.1 Caractéristiques morphométriques 

L’étude morphométrique d’un bassin versant est la mesure et l'analyse mathématique qui 
portent les caractéristiques de la forme (surface, périmètre, longueur, allongement, ...) et du 


relief (altitude, pente et profils). 


L’examen quantitatif des attributs morphométriques revêt une importance considérable pour 
évaluer et organiser le bassin en vue de la préservation des sols et des eaux. De nombreuses 
méthodes sont utilisées pour extraire les paramètres morphométriques (Gravellius, 
1914 ;Horton, 1945 ;Strahler, 1957 ; Roche 1960 ;; Krishnamurthy et al, 1996). Au cours 


des dernières décennies, l’augmentation des capacités de calcul informatique a permis la 
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numérisation des images satellitaires. Depuis ce temps là, les études morphométriques sont 
de plus en plus faites par les méthode de télédétection et de Systèmes d’Informations 
Géographiques (Nag, 1998; Rai et al, 2017 ; Magesh et Chandrasekar, 2014 ; Singh & 
Kanhatya, 2015 ; Prakash et al, 2016). 


Les attributs géographiques et géomorphologiques d'un bassin versant sont utilisés dans les 
examens hydrologiques, y compris l'évaluation du potentiel des eaux souterraines (Sreedevi 
et al., 2005) et la hiérarchisation des bassins versants (Patel et al, 2013), mais aussi pour 
anticiper des futurs scénarios des réponses hydrologiques combinée avec les paramètres 


géographiques (Angillieri, 2008). 
IL.2.1.1 La forme du bassin versant 


Le recours aux techniques de numérisation et à l’aide de l’application « ArcMap », la surface 
du bassin versant du Oued K’sob a été estimé à A=1480 km’ et le périmètre à P=185 km. 
La forme d'un bassin versant influe sur la réponse hydrologique, notamment sur la forme de 
l’hydrogramme produit à l’exutoire du bassin. Sur un bassin allongé, l’eau ruissellera sur 
des chemins plus longs et donc aura besoin de plus de temps pour arriver à la sortie du bassin 
que si le bassin versant à forme ramassée. Inversement, un bassin versant en forme compact 
a un temps de concentration plus rapide, ce qui entraîne un temps de concentration des eaux 
plus court et un débit de pointe plus élevé (Musy & Higy, 2004; Sassolas-Serrayet et al, 
2018). 

Divers indices morphologiques peuvent être utilisés pour caractériser les débits et comparer 
différents bassins versants. L'indice de forme de Gravelius (1914), KG, compare le périmètre 
d'un bassin versant au périmètre d'un cercle de même surface. Il est exprimé par l'équation 


suivante : 


P P 
Ke = —— = 0.282. — 14 
E NTA VA D 


Où Kg est l'indice de forme de Gravelius, A est la superficie du bassin versant [km?], et P 


est le périmètre du bassin versant [km]. 
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Figure 11.11 Localisation du bassin versnt Oued K’sob 39 
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L'indice est proche de 1 pour un bassin versant de forme circulaire et supérieure à 1 pour un 


bassin versant de forme allongée (figure IL.12) 


C Basin area 
== Basin perimeter 


CJ Source area 


= Source area length 


— Tributary streams 
=== Master stream 


GC=1 GC = 1.5 GC = 1.8 GC = 2.5 


Figure IL.12 Variation de l’indce de Gravelius en fonction de la forme du bassin (Sassolas- 


Serrayet et al, 2018) 


Le coefficient de Gravelius estimé pour le bassin de l’Oued K’sob, indique une forme 


allongée avec un Kc=1.35. 


IL.2.1.2 Rectangle équivalent 

Le concept de rectangle équivalent, ou rectangle de Gravelius, est utilisé pour calculer 
la pente dans un bassin versant. Cette méthode, introduite par Roche (1963), permet 
de comparer les pentes de différents bassins versants pour comprendre l'effet des 
caractéristiques de la pente sur l'écoulement. 

Par une transformation géométrique, la forme réelle du bassin versant, est transformée 
en un rectangle ayant la même surface, le même périmètre et la même répartition 
hypsométrique. Les lignes de contour deviennent des lignes droites parallèles au côté cour 
du rectangle. De ce fait, la répartition des sols, le couvert végétal et la densité de drainage 
restent inchangés entre les lignes de contour. 

Si L et l et représentent la longueur et la largeur respectives du rectangle équivalent, alors on 


peut écrire la surface et le périmètre du bassin par les équations suivantes : 


P=2.(L+1) (11.2) 
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La surface serait égale à : 


P=L.l (11.3) 


Une combinaison des Eq. (IL.1-3) permet d'éliminer les variables P et 1 et donne un polynôme 
du second ordre dont les solutions représentent la longueur, L, et la largeur, 1, du rectangle 


équivalent et qui auront pour expression : 


1 
L =. Kç. VTA. 


7 1+ |1- avecKg > 


(11.4) 


CIS 
VT. Kg vr 


1 | 2 l 2 
ETE | avecKe 2 = (11.5) 


Au bassin versant de l'Oued K’sob correspond un rectangle équivalent ayant une 


longueur trois fois et demi plus grande que la largeur (tableau I). Le rapport de 
proportionnalité entre la longueur et la largeur est aussi déterminé à partir du coefficient de 
Gravelius. Le rapport 3.5 correspond aux dixième et centième du coefficient de Gravélius 
(Tableau 1). 


Tableau II.1 : Dimensions du rectangle équivalent du bassin Oued K’sob 


Surfaces partielles (Xi)(km?) | Longueurs partielles Li (km) | Largeur 1 (km) | Longueur L (km) 
64.72 3.12 

273.70 13.19 

506.50 24.41 

289.10 13.93 

141.60 6.82 20.75 71.29 
139.70 6.73 

47.78 2.30 

16.81 0.81 


IL.2.1.3 Hypsométrie 


Les caractéristiques du système altimétrique jouent un rôle important dans la description 
d'un bassin versant et la compréhension de son comportement hydrologique, car la principale 
force qui affecte l'écoulement de surface est la gravité. Et, de ce fait, plusieurs processus 
hydrologiques dépendent directement des caractéristiques altimétriques. De même, les 
conditions climatiques locales sont influencées par la variabilité altimétrique et de 


l’écoulement résultant d’un événement pluvieux (Musy & Higy, 2004). 
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Tableau II.2 : Répartition des surfaces partielles par rapport aux altitudes 


Bornes 
Altitudes (m) | supérieures Surface partielle Surface cumulée 


(m) 


1050 
750-900 | 900 | 273.70 | 18.49%] 1,415.19 | 95.63% 
600 - 750 

| 585-600 | 600 | oio | 0.01%] 1,479.91 | 100.00% 


La courbe hypsométrique du bassin versant de l’Oued K’sob présente une forme concave 


indiquant un état de maturité du bassin versant (Figure II.13). 
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Figure IL.13 Courbe hypsométrique du bassin versant de l’Oued K’ sob 


L'élévation maximale d’un bassin versant représente la plus haute altitude au-dessus du 
niveau de la mer. Tandis que l'altitude minimale considère le point le plus bas, généralement 
l'exutoire. Ces deux paramètres deviennent particulièrement importants lors de l'élaboration 


d'équations impliquant des variables climatiques telles que la température, les précipitations 
42 


Chapitre IT Présentation de la zone d'étude 


et la couverture neigeuse. Ils sont également utilisés pour déterminent l'amplitude 


altimétrique du bassin versant et le calcul de sa pente. 


L’altitude minimale et maximale du bassin Oued K’sob correspondent aux altitudes 


respectives de : 585 et 1888 m (Figure IL.11). 


IL.2.1.4 L'altitude moyenne 

L'altitude moyenne d'un bassin versant peut-être déduite directement de la courbe 
hypsographique ou de la lecture de la carte topographique. L'altitude moyenne d'un bassin 
versant est souvent utilisée dans l'évaluation de certains paramètres hydro-météorologiques 


et peut être calculée à l'aide de la relation suivante : 


Sihi 
Hinoj = a - (11.6) 


Où 

Hmoy : Altitude moyenne du bassin [m] ; 

Si : Superficie de la zone comprise entre deux courbes de niveau [km2] ; 
h; : Altitude moyenne entre deux courbes de niveau [m] ; 


A : Superficie totale du bassin versant [km2]. 


Tableau IL.3: Calcul de l’altitude moyenne 


Altitudes (m) h; (m) Si (km?) (Sihi) /S 
585-600 592.5 0.10 0.04 
600 - 750 675 64.62 29.47 
750 - 900 825 273.70 152.58 
900 - 1050 975 506.50 333.69 
1050 - 1200 1125 289.10 219.77 
1200 - 1350 1275 141.60 121.99 
1350 - 1500 1425 139.70 134.52 
1500 - 1650 1575 47.78 50.85 
1650 - 1800 1725 15.49 18.06 
1800 - 1888 1844 1.32 1.64 

Somme 1,479.91 1,062.61 


Hmoy = 1,063 m 


II.2.1.5 L'altitude médiane 
Ce paramètre correspond à l'altitude lue sur l'axe x de la courbe hypsométrique 
correspondant à 50% de la surface totale du bassin versant (Figure II.13). Sa valeur est 


presque égale à l'altitude moyenne lorsque la courbe hypsométrique du bassin versant a une 
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pente régulière : en d'autres termes, lorsque la densité de probabilité des pentes du bassin 
versant présente une distribution symétrique. L’altitude médiane dans le bassin versant de 


l’Oued K’sob est égale 4: H504= 1168.9 m 


IL.2.1.6 Dénivelée simple 
Ce paramètre est défini comme étant la différence d’latitude Hs% et Hos%. 


D = Hs% — Hosy (1.7) 


Les paramètres Hs% et Hosy, peuvent être estimés par interpolation linéaire à partir du tableau 
altimétrique Tableau (11.2) 
D = 1490 — 755.5 = 734.5 m 
11.2.2 Indices des pentes 
IL.2.2.1 La pente moyenne du bassin versant 

La pente moyenne est paramètre important car il nous offre des informations sur la 
topographie du bassin versant. Il est considéré comme une variable indépendante et donne 
une bonne indication du temps de parcours du ruissellement direct et donc du temps de 
concentration t, et a une influence directe sur le débit de pointe après une pluie. 

Plusieurs méthodes ont été développées pour estimer la pente moyenne des bassins versants, 
toutes basées sur la lecture d'une carte topographique exacte ou approximative. Dans les 
années 1960, Carlier et Leclerc (1964) ont proposé une méthode qui consistait à calculer la 
moyenne pondérée des pentes de toutes les surfaces situées entre deux altitudes données. 
Une valeur approximative de la pente moyenne peut alors être estimée à l'aide de l'équation 


suivante : 


D.L 
im = (11.8) 
Où imest la pente moyenne [m/km] ou en [%o] ; 
L est la longueur totale des courbes de niveau [km] (a ne pas confondre avec la longueur du 
rectangle équivalent) ; 
D est l'équidistance entre les courbes de niveau [m] ; 
et A est la surface du bassin versant [km2]. 


im = 163.1%0 = 16.31% 


44 


Chapitre IT Présentation de la zone d'étude 


IL.2.2.2 L'indice de pente de Roche Ip 

Les méthodes précédentes donnent de bons résultats dans le cas de reliefs modérés et lorsque 
les lignes de contour sont simples et régulièrement espacées. Cependant, lorsque les courbes 
de niveau se tordent et tournent, il est difficile de déterminer leur longueur totale L. 

Afin de surmonter ce problème et les incertitudes qui résultent du lissage des courbes de 
niveau, Roche (1963) a proposé un indice de pente basé sur le rectangle équivalent et la 
courbe hypsométrique du bassin versant. L'idée est d'appliquer la méthode du rectangle 
équivalent à chaque ligne de contour dans le bassin versant, de sorte que chaque ligne de 
contour soit transformée géométriquement en lignes droites parallèles sur le rectangle 


équivalent. L'indice de pente i, est exprimé en pourcentage sous la forme : 


ip = D Pete) (1,9) 


Où nh; est la fraction de la surface totale A entre deux courbes de niveau consécutives de 
l'altitude a; et a;_,. L est la longueur du rectangle équivalent. 


ip = 4.60% 
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Figure IL.14 Carte des pentes du bassin versant Oued K’sob 


IL.2.2.3 L’indice de pente globale Ig 
Cet indice également basé sur l'équation décrit la distribution des altitudes du bassin versant 
(c'est-à-dire la courbe hypsométrique), est l'indice de pente global, qui est exprimé en m/km 


et défini comme : 
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Heo — H o 
ig = (110) 


Où Hs% et Hos% sont les fractiles de 5% et 95% de la courbe hypsométrique du bassin 
et L est la longueur du rectangle équivalent. 
ig = 1.03% 

IL.2.2.4 Dénivelée spécifique Ds 
L'indice i, décroit pour un même bassin lorsque la surface augmente, il était donc difficile 
de comparer des bassins de tailles différentes. La dénivelée spécifique D, ne présente pas cet 
inconvénient, elle dérive de la pente globale Ig en la corrigeant de l'effet de surface admis 
étant inversement proportionnel à VA 
Ds = ig. VACII. 11) 
D, = 396.63m 


La dénivelée spécifique ne dépend donc que de l'hypsométrie (D = H5% — Hos%) et de la 
forme du bassin (1/L). Elle donne lieu à une deuxième classification de l'O.R.S.T.O.M., 
indépendante des surfaces des bassins. Le bassin versant Oued K’sob présente un relief fort 


(Tableau II.4). 


Tableau IL.4 Classification du relief (Laborde, 2009) 


RI Relief très faible D, <10m 

R2 Relief faible 10m<D,<25m 
R3 Relief assez faible 25m<D,<50m 
R4 Relief modéré 50 m < D, < 100m 
R5 Relief assez fort 100 m < D, < 250m 
R6 Relief fort 250 m< D, < 500 m 
R7 Relief très fort 500 m< D, 


IL2.3 Caractéristiques du réseau hydrographique 

L'analyse du réseau de drainage consiste à quantifier l'organisation topologique des cours 
d'eau qui sont contenus dans un bassin versant (Smart, 1978). La recherche 
géomorphologique classique, dont Horton a été le pionnier (Horton, 1945), a démontré 
l'association entre les mesures topologiques de la taille des cours d'eau ou de la position du 
réseau (par exemple, l'ordre des cours d'eau) et les processus qui se déroulent dans les 
paysages fluviaux. La topologie du réseau de drainage a été liée aux propriétés des eaux de 


ruissellement du bassin (Kirkby, 1976 ; Rodriguez-Iturbe et Valdes, 1979 ; Naden, 1992 ; 
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Rice et al, 2016) et aux caractéristiques des habitats des cours d'eau qui affectent l'abondance 
et la diversité des espèces aquatiques (Vannote et al, 1980 ; Stenger-Kovacs et al, 2014). 
D'autre part, les systèmes d'information géographique (SIG) sont couramment utilisés pour 
effectuer des tâches d'analyse des réseaux de drainage. En effet, les méthodes automatisées 
d'analyse basées sur les SIG sont efficaces et reproductibles par rapport aux méthodes 
manuelles (Dawson et al, 2002). Les données numériques des cours d'eau utilisées dans les 
SIG sont généralement obtenues soit en numérisant les cours d'eau à partir de cartes 
existantes, de photographies aériennes et d'images satellites. Toutefois, ce sont les méthodes 
d'analyse du réseau de drainage basées sur les MNT qui sont le plus souvent mises en œuvre 
dans les logiciels SIG disponibles et sont donc plus largement utilisées dans la pratique. 
IL.2.3.1 Hiérarchisation du réseau 
En hydrologie, la structure des réseaux hiérarchiques dans les systèmes physiques a 
longtemps été caractérisée en se basant sur nombreuses méthodes. La méthode la plus simple 
et la moins ambiguë pour effectuer une classification topologique d'un réseau de drainage a 
été proposée par Horton en 1945, et modifié par Strahler en 1957. Le système de Strahler 
« Stream Order », qui est toujours le plus utilisé, est basé sur les principes suivants : 
- tout cours d'eau n'ayant pas d'affluent est dit d'ordre 1 ; 
- au confluent de deux cours d'eau de même ordre n, le cours d'eau résultant est d'ordre n+1 
- un cours d'eau recevant un affluent d'ordre inférieur garde son ordre, ce qui se résume 
par: nt+n=n+letn+m=max(n,m) 


Dans le cas de Oued K’sob, le réseau hydrographique du bassin est d’ordre 4 (Figure IT) 


Légende 


Ordre du cours 
d'eau 


0 5 10 20 Kilometers 
Loos | 1 1 1 | 


Figure IL.15 Classification du réseau hydrographique d’Oued K’sob selon Strahler (1957) 
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IL.2.3.2 Profil en long du cours d'eau principal 

La pente hydrographique, ou profil en long du cours d'eau, peut-être déterminée sur 
la carte ou mesurée sur le terrain par un nivellement de précision. Cette pente exprimée 
généralement en m/km conditionne : 

- La vitesse de l'eau dans le chenal ; 

- La vitesse de l'onde de crue ; 

- Le tirant d'eau de la rivière : pour un même débit et une même largeur une rivière 
plus pentue a une vitesse d'écoulement plus grande et donc, généralement, une 
profondeur plus faible. 

La pente hydrographique varie plus ou moins irrégulièrement pour un même cours d'eau 
selon les structures géologiques traversées et diminue en général d'amont en aval (forme 
concave des profils en long). La figure suivante représente le profil en long du cours d’eau 


principal de l’Oued K’sob. Elle est créée en utilisant le logiciel Google Earth Pro. 


Graphique : Min. Moy.. Max Elévation : 603. 870, 1193 m 
Totaux des plages de valeurs : Distance : 78 km Gain/perte d'élévation : 819 m, -229 m 


78 km 


Figure IL.16 Profil en long du cours d’eau principal de l’Oued K’sob (réalisé à l’aide du 


logiciel Google Earth Pro) 


IL.2.3.3 La longueur du cours d’eau principal 

La longueur du cours d'eau principal (L) est la distance curviligne depuis l'exutoire jusqu'à 
la ligne de partage des eaux, en suivant toujours le segment d'ordre le plus élevé lorsqu'il y 
a un embranchement et par extension du dernier jusqu'à la limite topographique du bassin 
versant. Si les deux segments à l'embranchement sont de même ordre, on suit celui qui draine 
la plus grande surface. 


Le calcul donne pour Oued K’sob une longueur L = 80.83 km 
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IL.2.3.4 La densité de drainage 

D’après Horton (1945), la densité de drainage est la longueur totale de tous les cours d'eau 
et rivières d'un versant divisée par la superficie totale du bassin de drainage. C'est une mesure 
de la qualité ou de la faiblesse du drainage d'un bassin versant par les canaux des cours d'eau. 


Ce paramètre s’exprime par la formule suivante : 


L; 
Da = 211.12) 


Avec : 

Da : Densité de drainage [km/km2] ; 

L; : Longueur des cours d'eau [km] ; 

A : Surface du bassin versant [km2]. 

La densité de drainage dans le bassin versant de l’Oued K’sob est égale à 0.29 km/km2, ceci 
traduit que le bassin présente dans son ensemble une formation géologique perméable, dont 


l’écoulement est plutôt limité et centralisé ainsi que l’infiltration est augmentée. 


Tableau II.5 Calcul de la densité Drainage 


Ordre | Longueur d'ordre Ni | Densité de Drainage Dd 
1 235.53 0.16 
2 89.13 0.06 
3 79.58 0.05 
4 31.59 0.02 
Total 435.82 0.29 


IL.2.3.5 La densité hydrographique 


La densité hydrographique représente le nombre des cours d’eau par unité de surface : 


_ XN; 


F 
A 


(11.13) 
Où : 


F: Densité hydrographique [km2] ; N; : Nombre des cours d’eau ; A : Superficie du bassin 
[km2]. 
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Tableau II.6 Calcul de la densité hydrographique 


Ordre | Nombre d'ordre Ni Densité hydrographique F 
1 115 0.08 
2 36 0.02 
3 45 0.03 
4 12 0.01 
Total 208 0.14 


La densité hydrographique qui est égale à F=0.14, est une faible densité et donc le bassin 
d’étude présente une région à substratum perméable, à couvert végétal important et à relief 


peu accidenté. 


II.2.3.6 Constante de stabilité du cours d’eau 
Selon Schumm, la valeur inverse de la densité de drainage est une constante de stabilité du 


cours d’eau : 


1 
C, = — (11.14 


C, = 3.40 Km?! 
IL.2.3.7 Rapport de confluence 
Sur la base de la classification des cours d'eau, (Horton 1932; Schumm 1956) ont établi cette 


loi : 


Re = (11.15) 


Avec : 
R, : Rapport de confluence des cours d’eau ; 
N; : Nombre des cours d'eau d'ordre suivant ; 


N;41 : Longueur des cours d'eau d'ordre suivant. 
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Tableau IL.7 Calcul des rapports de confluence 


Ordre Nombre d'ordre N; Re 
1 115 3.19 
2 36 0.88 
3 41 10.25 
4 4 - 


Le rapport de confluence est un nombre sans dimension exprimant le développement du 
réseau de drainage. Il varie suivant l'ordre considéré. C'est un élément important à considérer 
pour établir la différence d'une région à une autre. Selon (Strahler 01964), le RQ varie de 3 à 


5 pour une région où la géologie n'a aucune influence. 


IL.2.3.8 Rapport des longueurs 


Sur la base de la classification des cours d'eau, (Horton 1932; Schumm 1956) ont établi 


cette loi : 


Li 


Lisa 


R, = (11.16) 


R, : Rapport des longueurs des cours d’eau ; 
i: Ordre d'un cours d'eau u varie entre 1 etn; 


Li : Longueur moyenne des cours d'eau d'ordre i ; 


Tableau II.8 Calcul des rapports de la longueur 


Ord Longueur Li R 
rp 183 264 
2 89.13 

79.58 L. 
31.59 2. 


II.2.3.9 Le coefficient de Torrentialité 


C’est le rapport entre la fréquence des cours d’eau d’ordre 1 avec la densité de drainage 


Ci = Da X F (11.17) 
C, =0.29*0.14= 0.04 km/km“ 
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Cette faible valeur traduit que le bassin ne représente pas des caractéristiques 
morphométriques adaptées à l’écoulement, ainsi qu’il est constitué par des formations 


perméables dont le taux de l’infiltration est important. 


IL.2.3.10 Temps de concentration 

Le temps de concentration est défini comme étant la durée que met une goutte d'eau tombée 
à un point élevé du bassin versant, pour atteindre l'exutoire, il est exprimé par la formule 
suivante de Giandott : 


AVA + 1.5L 
te = —— ~ (II. 18) 


0.8) Hmoy — Hmin 


A : Superficie du bassin [km?]; 

L : Longueur du thalweg principal [km]; 

Hmoy : Altitude moyenne [m]; 

Hmin : Altitude minimale [m]; 

La détermination du temps de concentration nous permettra de mieux connaitre la réponse 
hydrologique du bassin ; le temps de concentration des eaux est non seulement fonction de 
la surface, mais aussi d'autres paramètres tels que la lithologie, la pluviométrie et le relief. 
D’après les résultats obtenus, le bassin versant de l’Oued El K’sob présente un temps de 


concentration de l’ordre 9.15 h. 


II.2.3.11 Vitesse d’écoulement de l'eau 


Elle est donnée par l'expression suivante 


Vp =— (11.19) 


V, : Vitesse d'écoulement de l'eau [m /s]; 
L : longueur de talwegs principal [m] ; 


tc : temps de concentration en [s]. 


V, = 2.45m/s 
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II.3. Barrage K’sob 
Le barrage du K'sob du même nom que l’oued sur lequel il a été construit entre 1934 à 1939. 
Le barrage du type à voûtes multiples, de 46 m de haut se trouve à 15 Km au nord de la 


Wilaya M'sila (Figure II.17). 


Figure IL.17 Barrage K’sob (photo : Guesri 2016) 


D'une capacité de 29.5 millions de mètres cubes d'eau destinés à l'irrigation de 13 000 ha, le 
barrage ne permet, aujourd'hui, que l'irrigation du tiers de cette superficie en raison d'un état 
d'envasement avancé. Quelques chiffres sur le barrage K’sob sont résumés dans le tableau 


suivant : 


Tableau IL9 Le barrage en chiffre (Source : ANBT) 


Paramètre Chiffre 

Année de construction 1934 

Année de Mise en eau 1940 

Année de surélévation 1972-1976 
Capacité 29.5 hm° 
Capacité dernier levé 11,84hm° 
Apport moyen annuel 51,50 hm° 
Envasement annuel 0,70 hm?/an 
Surface du bassin versant 1 460 km? 
Hauteur 46,00 m 
Longueur 280,00 m 

Côte retenue Normale (R.N) 604,00 mNGA 
Côte Plus Hautes Eaux (P.H.E) 605,50 mNGA 
Déversoir à Seuil libre 850,0 m°/s 


L’analyse bi-date du barrage K’sob à partir de deux images Alsat-2A (Figure II.18), fait 
ressortir l’évolution du niveau d’eau de ce barrage, ayant atteint au courant du mois de février 
2015, un fort taux de remplissage et ce, en comparaison avec le taux de remplissage au mois 
d’octobre 2013. 

53 


Chapitre II 


Présentation de la zone d'étude 


Date de prise de vue : 16Octobre 2013 


RFE NX 


Barrage El K'sob 


Date de prise de vue : 10 Février 2015 


N’ = MIE ey) 


| 


Barrage El K'sob 


de vue : 10 Février 2015 
LR 


Limite de la tranche d’ eau du Barrage El K’sob en Octobre 2013 


Limite de la tranche d’eau du Barrage El K’sob en Février 2015 


Figure IL.18 a) Vue d'ensemble du barrage K'sob (M'sila) à partir d'image Alsat-2A, 
b+c)Extrait de l'image Alsat-2A du 10/02/2015, centrée sur le barrage K'sob les engins de 


dévasement sont identifiés en jaune d) caractérisation du changement du niveau d'eau entre 


2013 et 2015 


L’analyse comparative des deux prises de vue du satellite Alsat-2A, s’appuyant sur la bonne 


caractérisation spectrale des tranches d’eau, grace a la bonne discrimination spectrale fournie 


par ses images multispectrales, a permis de délimiter les surfaces d’eau du barrage. 


L’accroissement de la surface d’eau ainsi identifiée, est lié directement à une augmentation 


du volume d’eau stocké dans le barrage. Actuellement, le barrage fait l’objet d’une nouvelle 


opération de dévasement ; après celle de 2001 ; qui concernera un volume plafond de l’ordre 


de 10 millions de mètres cubes. Enfin, il est à noter, que les engins destinés à cette opération 


de dévasement sont visibles sur les 02 prises de vue Alsat-2A. 
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Chapitre HI: Variabilité climatique 


III.1 Introduction 

La connaissance fine des conditions météorologiques est importante pour les études 
d’impact de la variabilité climatique sur les ressources en eau ou sur les phénomènes 
naturels comme l’érosion hydrique. Pour cette fin, il est nécessaire de quantifier et 
d’estimer tous les paramètres qui interviennent dans le processus hydrologique. 

Les précipitations constituent le facteur essentiel intervenant par leur hauteur d’eau afin de 
déterminer l’abondance pluviale annuelle et saisonnière et, par leurs totaux journaliers, de 
distinguer les averses génératrices des crues et des sédiments. 

Ce chapitre est consacré pour la caractérisation et l’identification du climat de la région 
d’étude Oued K’sob, et plus spécialement l’étude des précipitations et la température et 
leur irrégularité. L’analyse statistique des différents paramètres météorologiques (Pluie, 
Température, Evapotranspiration, Humidité et Vent) est effectuée selon la disponibilité des 


mesures climatiques. 


HI.2 Collecte des données 

La consultation des archives de l’ Agence Nationale des Ressources Hydrauliques 
(ANRH) et l’Office National de la Météorologie (ONM) nous a permis d’acquérir les 
données recueillies à la station météorologique St5 de Bordj Bou Arreridj (BBA) (Fig. II- 
2). Cette dernière se trouve à la latitude de 36°04 N et la longitude de 04°46 E, et à une 
altitude de 922 m. Les périodes d’étude et résolutions de mesures de données sont décrites 
dans le tableau suivant : 
Tableau IIL.1. Récapitulatif des données fournies par ONM et mesurées au niveau de la 


station météorologique BBA (St5) 


Paramètre Période Résolution 
1990-2010 | 3h 
Pluie 
2007-2015 15min 
Température moyenne 1990-2010 3h 
Température minimale 1990-2010 | 3h 
Température maximale 1990-2010 | 3h 
Evaporation 1996-2005 24h 
Humidité relative minmale 1996-2005 mois 
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Humidité relative moyenne 1996-2005 mois 


Humidité relative maximale | 1996-2005 mois 
Vents 1980-2009 mois 


Légende 
@ Station météo 


è Ville 


Réseau hydrographique 


| Barrage K'sob 


0 5 10 20 Kilometers 
Lun nn | 1 1 1 J 


Figure III.1 Localisation de la station météo St5 de BBA 
I.3 Températures 
La température de l’air, couramment mesurée, est importante pour l’estimation du bilan 
hydrologique et de l’évapotranspiration. Elle est en interaction directe avec d’autres 
facteurs météorologiques (insolation, vitesse et turbulence du vent, précipitations) et 
biogéographiques comme l’évapotranspiration et le déficit d’écoulement. 
La figure IIL2 représente la variation des moyennes mensuelles des températures 
moyennes, minimales et maximales qui correspondent à une courte durée de 21 ans 
d’observation de Janvier 1990 à Décembre 2010. 


40 —4—T max (°C) —=—T min (°C) —e—T moy(°C) 
30 


20 


-10 


Sep Oct Nov Dec Jan Feb Mar Apr May Jun Jul Aug 


Figure IIL.2 Variations mensuelles des températures horaires à la station de Bordj Bou 


Arreridj (1990-2010) 
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À l’échelle horaire, le maximum des températures est atteint en Juillet avec une moyenne 
de 34.7 °C et le minimum se situe en Janvier à l’environ de -7.5 °C. 

Il s’avère aussi qu’à partir de Février et jusqu’à Aout, les températures moyennes 
mensuelles tendent à s’accentuer en présentant des valeurs allant de 7.5 °C à 27.4 °C, par 
contre la saison comprise entre Septembre et Janvier représente des températures 


moyennes mensuelles qui varient entre 26.9 °C et 6.4 °C. 


III. 4 Evapotranspiration 

L'évapotranspiration est le processus d'évaporation de l'eau des feuilles par la transpiration 
de la plante pendant la photosynthèse. Elle varie en fonction d'une multitude de facteurs 
comme le vent, la température, l'humidité et la disponibilité de l'eau. 

L'évaporation et la transpiration se produisent simultanément et il n'y a pas de moyen facile 
de distinguer entre les deux processus. Outre la disponibilité de l'eau dans la couche arable, 
l'évaporation d'un sol cultivé est principalement déterminée par la fraction du rayonnement 
solaire qui atteint la surface du sol. Cette fraction diminue au cours de la période de 
croissance à mesure que la culture se développe et que le couvert végétal ombrage de plus 
en plus la surface du sol. Lorsque la culture est petite, l'eau est principalement perdue par 
l'évaporation du sol, mais une fois que la culture est bien développée et recouvre 
complètement le sol, la transpiration devient le principal processus. Dans ce qui suit, nous 
présentons les méthodes classiques permettant d’estimer l’évapotranspiration potentielle et 


réelle. 


IIL.4.1 Méthode de Thornthwaite 
Afin d’estimer l’évapotranspiration, Thornthwaite (1948) a proposé une formule basée 


essentiellement sur les températures de l'air : 
10 T\* 
ETP(mm) = 16 (=) K (11.1) 


Où : ETP est l’évapotranspiration potentielle mensuelle en (mm) ; 
T est la température moyenne mensuelle en degrés celsius en (°C) ; 


I est un indice thermique annuel : somme des 12 indices mensuels calculés à partir des 


51 
températures moyennes mensuelles selon la formule : i = (<) (111.2) 
so 111.3 
a= 100 . (111.3) 
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K est un facteur correctif fonction de la saison et de la latitude de la station. 


Tableau II.2. Les valeurs mensuelles du coefficient correctif K 


Mois Jan Fev Ma Av Mai Jui Juil Août Sep Oct Nov Déc 


k 0.87 0.85 1.03 1.10 1.21 1.22 1.24 1.16 1.03 0.97 0.86 0.84 


Les valeurs de l’évapotranspiration correspondant aux températures moyennes mensuelles 


pour la station de BBA sont représentées dans la figure II.3. 


O Température (°C) mm ETP (mm) 
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Figure IIL.3 Répartition mensuelle de l’évapotranspiration et de la température. 
On observe que les valeurs d’évapotranspiration mensuelle sont en relation avec celles des 
températures mensuelles, c’est-à-dire quand les températures augmentent, l’ETP atteint ses 
valeurs maximales et contrario. Cette analyse montre l’existence de deux périodes : 
= Une période sèche pendant laquelle les températures sont maximales, généralement 
entre Avril et Octobre. 
= Une période où les températures sont inférieures à partir du mois de Novembre 
jusqu’à Mars. 
De façon générale, les valeurs de l’évapotranspiration potentielle deviennent plus faibles 
pendant les mois pluvieux qui débutent à partir du mois de novembre jusqu’au mois 
d’Avril. 
III.4.2 Méthode de Turc 
Cette méthode nécessite la connaissance des températures de l'air et des précipitations 
moyennes mensuelles. L’évapotranspiration réelle se déduit à partir d’une relation 
applicable à toutes les régions, donnant de bons résultats sur 254 bassins dans le monde 


(Turc, 1961). Cette formule est la suivante : 
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P 
ETR = ——— (III. 4) 


ETR : Evapotranspiration réelle annuelle (mm) ; 


P : Précipitations moyennes annuelles (mm) ; 


L = 300 + 25T + 0.05 T? (11.5) 
Avec : 
T : Température moyenne annuelle (°C). 
Tableau IIL.3. Calcul de PETR selon Turc 
P(mm) T°C) L ETR (mm) 
340 15.91 899,11 332.91 


Donc, ETR =332.91 mm, ce qui représente presque 98 % des précipitations. 
IIL.4.3 Méthode de Coutagne 
La méthode de Coutagne (1935) qui sert à estimer le déficit d’écoulement, est basée sur les 
mêmes hypothèses que la méthode de Turc pour l'estimation de la valeur annuelle 
moyenne de l'évapotranspiration réelle (Sakkas, 1985 ; Mimikou et al, 2016). La formule 
empirique donnée par Coutagne est la suivante : 
D =P — AP? (111.6) 


1 
0.8+0.14 T 


Où : A= (111.7) 
D : Déficit d'écoulement (m/an) ; 

P : Hauteur moyenne annuelle des précipitations (m/an); 

T : Température moyenne annuelle (°C); 


La formule précédente (II.6) est applicable dans la condition suivante : 
1 1 
—<P<— CHI. 8) 


Tableau II.4. Calcul du déficit d’écoulement selon Coutagne 
1/8A (mm) | 1/2A (mm) | T(°C) D(mm) | Ecoulement+ infiltration (mm) 
378 1513.7 15.91 301.81 38.18 mm soit 11.23% 
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HI.5 Evaporation 
L’évaporation observée à l’aide d’un bac enterré type COLORADO à la station 
météorologique BBA pour une période de 2lans, en allant de 1990/1991 jusqu’à 


2009/2010 a donné les moyennes mensuelles représentées dans la courbe suivante : 


y =+7.4885x +255.65 © 


Evaporation mensuelle moy (mm) 
— 
Un 
S 


ce K SR s OC ory EASE SEER 
9 9 SV MENNEN NN 
Années 


Figure IIL.4 Répartition mensuelle de l’évaporation de la station BBA (1990-2010) 
On note une réduction moyenne de l’ordre de 7,5 mm par an, ce qui représente 3.95 % de 
la moyenne annuelle estimée à 177 mm. Une valeur minimale d’évaporation mensuelle de 


100 mm est estimée en 2003/2004 et une autre maximale qui atteint 255 mm en 1993/1994. 


IIL.6 Diagramme pluvio-thermique 
Le diagramme pluvio-thermique de Bagnouls et Gaussen (1957) peut représenter une 
méthode utile pour mettre en évidence les limites climatiques qui peuvent fournir des 
informations sur le type de climat dans un site donné. A cet effet, cet outil a démontré un 
assez bon pouvoir prédictif. Également appelé "diagramme ombrothermique ", il fournit 
une vue d'ensemble des températures et des caractéristiques des précipitations pour une 
station donnée. 
Ce diagramme représente sur un même système d’axes, la variation des températures et 
celle des précipitations moyennes mensuelles. On porte les valeurs des températures 
moyennes à l’échelle double des précipitations moyennes mensuelles sur l’axe des 
ordonnées, les différents mois de l’année figureront sur l’axe des abscisses (Figure II.5) 
Le climat de la région se caractérise par deux périodes : 

- Une Période sèche relativement courte qui dure du mois de Juin jusqu'au mois de 


Septembre. 
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- Une Période humide qui commence du mois d'Octobre jusqu'au mois de Mai. 


Période 
sèche 


Pluie moyenne 2P(mm) 
Température moyenne (°C) 


sept oct nov déc janv févr mars avr mai juin juil août 


Figure IILS Diagramme pluvio-thermique 


Généralement, la période humide pour la station étudiée est marquée par des températures 
faibles par rapport à la température moyenne annuelle (15.91 °C) notamment pour le mois 
de Janvier (6.40 °C), et des hauteurs de précipitations maximales observées surtout entre 
les mois d’Avril et Mai. 

Les périodes sèches sont caractérisées par la rareté des pluies surtout pour les mois de 
Juillet et Août, au cours desquels les températures atteignent leurs valeurs maximales (26- 
27 °C). 

IIL.7 Humidité relative et Vent 

L'air sous haute pression se déplace vers les zones de basse pression. Plus la différence de 
pression est importante, plus l'air circule rapidement. Le vent est décrit avec sa direction et 
sa vitesse. La direction du vent est exprimée comme la direction d'où souffle le vent. Par 
exemple, les vents d'est soufflent d'est en ouest, tandis que les vents d'ouest soufflent 
d'ouest en est. Les vents ont différents niveaux de vitesse, tels que "brise" et "coup de 
vent", selon la vitesse à laquelle ils soufflent. La vitesse du vent est mesurée avec un 
anémomètre. 

Les tableaux suivants montrent les variations mensuelles de vent et humidité relative de la 
station de BBA pour nombreuses années. 


Tableau HI.5 Vents 1980-2009/ Source : ONM 
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Sept Oct Nov Déc Jan Fév Mars Avr Mai Juin Juillet Août 


Vinoy 
(m/s) 26. 26, DED 27h. 7 2H -. 32 3,3 3,2 3,3 3:2 3 
Vmoy 
(km/h) 10,08 9,36 9,72 9,72 9,72 9,72 11,52 11,88 11,52 11,88 11,52 10,8 


Tableau IIL6 Humidité relative mensuelle pour la période : 1996-2005/ Source : ONM 


Sept Oct Nov Déc Jan Fév Mars Avr Mai Juin Juillet Août 


Humidité 


; 56 61 71 75 75 69 60 58 54 44 39 43 
moy (%) 


Humidite 30 36 49 54 54 44 35 34 30 2 19 22 
min (%) 


Humidité 82 83 89 91 91 89 84 83 79 70 61 68 
max (%) 


II.8 Etude de l’irrégularité des précipitations et des événements extrêmes (Courbes 
IDF) 
Les valeurs extrêmes des variables hydro climatiques, présentent un intérêt tout particulier 
lorsqu’on parle en termes de risque. L'estimation de la récurrence de pluies extrêmes 
fournit des éléments indispensables pour la construction d'infrastructures telles les digues 
et les ouvrages d'assainissement, afin de protéger efficacement la population et leurs biens 
(Benkhaled, 2007). 
Les statistiques associées aux précipitations extrêmes jouent un rôle majeur dans la 
conception des structures hydrauliques et les études hydrologiques. En particulier, 
l'estimation des courbes IDF (intensité-durée-fréquence) s'est principalement concentrée 
sur l'analyse des enregistrements historiques des précipitations et sur l'estimation de 
l'intensité d'un événement pluvieux extrême associé à une durée de pluie donnée et à une 
probabilité de dépassement (ou période de retour). 
La construction des IDF suit généralement la nature de l’aménagement hydraulique à 
réaliser. Elles nécessitent la réalisation de lois de probabilité ajustées à des séries 
pluviométriques aux échelles temporelles similaires aux temps de concentration du bassin 
en amant de l’ouvrage projeté. 

L'acquisition des précipitations avec une résolution fine, comme au pas horaire ou 
moins, constitue la première difficulté car la plupart des séries de mesure de la 


pluviométrie sont issues de pluviomètres traditionnels conçus pour un enregistrement au 
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quotidien. De part le monde, les appareils de mensure automatisés à haute précision sont 
rares. Cependant le bassin d’étude Oued K’sob a été équipé d’un tel appareil en 1990. 
III.8.1 Notion d’averse et d’intensité 

Une averse est un ensemble de pluies liées a une perturbation météorologique. La durée 
d’une averse varier de quelques minutes à plusieurs jours, Elle peut affecter quelques 
Kilomètres carrés (précipitations convectives) à quelques milliers de kilomètres carrés 
(Pluies cycloniques). Une averse peut être donc défini comme un épisode de pluie continue 
ayant une intensité variable (Musy et Higy, 2010). Une averse peut être caractérisée par 
plusieurs paramètres qui peuvent avoir, au sein de la même pluie, citons notamment : 

- La hauteur totale de pluie, 

- La durée, 

- L’intensité moyenne, 

- Les intensités maximales sur des intervalles de temps quelconques, 

- La distribution d'intensité instantanée i(t). 
L'intensité est définie comme le rapport entre la quantité totale de pluie (hauteur des 
précipitations) tombant pendant une durée donnée. Elle est exprimée en unités de hauteur 
par unité de temps, généralement en mm par heure (mm/h). Elle est exprimée par la 
relation suivante : 

h 
i=- (II.9) 

Où: 
i: intensité moyenne de la pluie [mm/h, mm/min] ou ramenée à la surface [1/s.ha], 
h : hauteur de pluie de l'averse [mm], 
t: durée de l'averse [h ou min]. 
En réalité, l'intensité d’une averse change durant l’épisode pluvieux et il n'est pas possible 
de simuler un événement de pluie comme étant complètement continu. Par conséquent, une 
averse est généralement saisie pour un modèle déterministe. Par exemple, un pluviomètre à 
auget basculant mesure la hauteur de pluie toutes les 10 ou 15 minutes. Une représentation 


graphique d'un tel événement est appelée "Hyétographe de pluie" comme illustrée dans la 


figure II.6. 
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Figure IIL.6 Hyétogramme enregistré en 12/10/2014 à la station Sts BBA 
IIL.8.2 Période de retour : 
En hydrologie, une période de retour, également appelée intervalle de récurrence ou 
intervalle de répétition, est le temps moyen estimé entre deux événements hydrologiques 
tels que les pluies, les inondations ou le débit d'une rivière pour se produire (Handbook, 


1996; Peres et Cancelliere, 2016). 


La probabilité au non dépassement de la valeur, Xp, d’une variable hydrologique est 


donnée par : 
F,(X,) = P(X <X,) (11.10) 


La probabilité d’occurrence d’un tel événement est liée au temps de retour par l’équation 
suivante (Chow et al. 1988) 


1 


~ Prob(X > Xp) Ce) 


III.8.3 Courbes IDF et leurs utilisations : 
Les courbes IDF (Intensité-Durée-Fréquence) décrivent la relation entre l'intensité des 
précipitations, la durée des précipitations, et la période de retour (ou son inverse, la 


probabilité de dépassement) (Figure III.7). 


64 


Chapitre IT Variabilité climatique 


Inténsité de Pluie (pouce/h) 


Durée (Minutes) 


Figure IIL.7 Représentation schématique des courbes IDF (Pit et clark, 2002) 


Les courbes IDF sont couramment utilisées dans la caractérisation du climat et aussi dans 
la conception de systèmes hydrologiques, hydrauliques et de ressources en eau. Les 
courbes IDF sont obtenues par l'analyse de la fréquence des observations des 


précipitations. 


IIL.8.4 Construction des courbes IDF 
La construction des courbes IDF (Intensité-Durée-Fréquence) de la station St5 BBA, entre 
1990 et 2010 (la période disponible) comporte plusieurs étapes que nous décrivons ci- 


après: 


- Extraction des intensités horaires maximales annuelles aux différents pas de temps 
(1h, 3h, 6h, 9h, ...24h) ; 

- Ajustement de chacune des séries obtenues à une loi de probabilité. ; 

- Choix du meilleur ajustement à l’aide du test d’adéquation. Dans cette étude, 
l’adéquation est examinée par le test de Kolmogorov-Smirnov ; 

- Pour chaque pas de temps, établir différentes périodes de retour. 

- A chaque période de retour, T, on fait correspondre le pas de temps t et l’intensité 
de retour correspondante Ir ; 

- Pour chaque période de retour, faire un ajustement par régression simple aux points 


obtenus (Ir, 1/t). Les paramètres d’ajustement donnent les coefficients de Montana. 
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Une représentation simple et pratique sera donnée suivant la figure ci-dessous. 


Cumulative Distribution Function F(x) 


Figure II.8 Schéma de Dérivation des courbes IDF (Sun et al, 2019) 


III.8.4.1 L’ajustement des pluies horaires maximales annuelles 

L'analyse fréquentielle d'une longue série pluies maximales permet d’estimer le temps de 
retour d'une valeur particulière. Cette prédiction repose sur la définition et la mise en 
œuvre d’un modèle fréquentiel qui est une équation décrivant (modélisant) 
le comportement statistique d’un processus. Ces modèles décrivent la probabilité 
d’apparition d’un événement de valeur donnée. C’est du choix du modèle fréquentiel (et 
plus particulièrement de son type) que dépendra la validité des résultats de l’analyse 
fréquentielle (Musy et al, 2014). 

Des modèles fréquentiels très souvent utilisés pour décrire le comportement statistique des 
valeurs extrêmes sont : la distribution statistique de Gumbel, Log-Gumbel et la loi Loi 
d'extremum généralisée (GEV). 

L'outil utilisé pour l’ajustement des pluies horaires maximales annuelles est celui de 
l’Hydrolab2018. ce logiciel libre, a été développé par J.P. Laborde, professeur à 
l'université de Sofia Antipolis, Intégré sous Excel, évolué en plusieurs versions : 98.2, 
2010 et 2018 (Laborde et Mouhous., 1998 ; Laborde, 2010). il permet une utilisation simple 
d'outils hydrologiques. 
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III.8.4.2 L’ajustement des pluies par la loi de Gumbel 

Un modèle fréquentiel très souvent utilisé pour décrire le comportement statistique des 
valeurs extrémes est la distribution statistique de Gumbel (loi double exponentielle ou loi 
de Gumbel). La fonction de répartition de la loi de Gumbel F(x) s’exprime de la manière 
suivante : 


-5 


F(x) =e" (111.12) 


Avec la variable réduite suivante : 
X — 
b 


Où a et b sont les paramètres du modèle de Gumbel. 


u = (111.13) 


š à a = 7 ae š —p7u 
La distribution s’écrit alors de la manière suivante : F(x) = e~© 


etu = —In(—In(F(x))) (1.15). 


(111.14) 


L’avantage d’utiliser la variable réduite est que l’expression d’un quantile est alors linéaire 
:x = atbu, (11.16). 

En conséquence, dès lors que les points de la série à ajuster peuvent être reportés dans un 
système d’axes x - u il est possible d’ajuster une droite qui passe le mieux par ces points et 
d’en déduire les deux paramètres a et b de la loi. Il existe différentes méthodes 
d’ajustement : méthode graphique (ajustement à l’œil ou à l’aide d’une régression 


statistique), méthode des moments...etc. 


En pratique, il s’agit essentiellement d’estimer la probabilité de non dépassement F(x:;) 
qu’il convient d’attribuer à chaque valeur xi. Il existe de nombreuses formules d’estimation 
de la fonction de répartition à l’aide de la fréquence empirique. Elles reposent toutes sur un 
tri de la série par valeurs croissantes permettant d’associer à chaque valeur son rang r. 

Des simulations ont montré que pour la loi de Gumbel, il faut utiliser la fréquence 


empirique de Hazen : 
r — 0.5 
n 


(11.17) 


Où r est le rang dans la série de données classée par valeurs croissantes, n est la taille de 
l’échantillon, x[r] la valeur de rang r. 
Rappelons encore que le temps de retour T d'un événement est défini comme étant l'inverse 


de la fréquence d'apparition de l'événement. Soit : 
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T 


I= F) 


1 


(111.18) 


Aux différents pas de temps, les résultats d’ajustement de la série des intensités horaires 


maximales annuelles à la loi de Gumbel sont affichés dans la figure III. Rappelons que ces 


intensités sont mesurées entre 1990 et 2010 au niveau de la station météorologiques St5 de 


BBA 
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Figure III.9 Ajustements des pluies horaires maximales annuelles par la loi de Gumbel 


aux différents pas de temps 
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IIL.8.4.3 L’ajustement des pluies par la loi de Log-Gumbel 


Si la variable aléatoire Ln(x) suit une loi de Gumbel de paramètre a et b alors x suit un loi 


log-Gumbel ou de Fréchet. La loi log-Gumbel présente une dissymétrie encore plus 


marquée que la loi de Gumbel, elle s’adapte ainsi aux échantillons très dissymétriques. 


La figure IIL10 montre les graphes d’ajustement de la série des intensités horaires 


maximales annuelles à la loi de Log-Gumbel 
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Figure IIL.10 Ajustement des pluies horaires maximales annuelles par la loi de 
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IIL.8.4.4 L’ajustement des pluies par la loi GEV 
Jenkison (1955) a montré que les lois de distribution des extrêmes pouvaient se mettre sous 
une forme unique : 


KG a) 


F(x) = se (III. 19) 


Trois paramètres interviennent : x, le paramètre de position, s le paramètre d’échelle et k 
le paramètre de forme. Sous cette formulation unique on retrouve en fait les trois types de 


distribution selon les valeurs prise par k : 


Si k est positif, on retrouve une loi proposée par Weibull (GEV de type III) où x est défini 
sur l’intervalle | — co, x4[ Cette borne supérieure de l’intervalle de définition se déduit des 
paramètres s, k et x, par la relation: xg = = + Xo. En fait ce type de distribution se 


rencontre rarement en hydrologie et on la réserve à des études de valeurs extrêmes 


minimales (température, étiages...). 


1/k (x-x9) 
k(x-x, = P 
— rozo] >e s . On retrouve donc la loi de 


Si k est égal à zéro on constate que [1 


Gumbel (GEV de type I) vue précédemment et x est défini sur l’intervalle | — 00, +20[ 


Enfin, si k est négatif, on trouve une loi proposée par Frechet (GEV de type II) où x est 
défini sur l’intervalle |xç, +o[ Cette borne inférieure de l’intervalle de définition se déduit 


des paramètres s, k et x, par la relation : x9 = z +xo (111.20) 


Ce type de loi se rencontre souvent en hydrologie (crues, précipitations extrêmes...), mais 
la difficulté va résider dans l’estimation du paramètre de forme qui est très soumise aux 
aléas de l’échantillonnage. Il parait très hasardeux de caler la valeur de k sur un seul 
échantillon. Il est plus raisonnable de choisir a priori la valeur de k d’après une étude 
régionale. Koutsoyiannis (2004) a étudié 169 séries longues de pluie extrêmes (Etats-Unis, 
Royaume-Uni, France, Italie, Grèce) et il propose de fixer k à -0.15. Pour les pluies 
extrêmes d’ Algérie (508 séries de 43 années en moyenne), pour l’ Algérie, Laborde (2017) 


a proposé de retenir k=-0.073. 


La figure suivante représente les résultats d’ajustement de la série des intensités horaires 


maximales annuelles à loi des extrêmes généralisée GEV. 
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Figure III.11 Ajustement des pluies horaires maximales annuelles par la loi GEV aux 


différents pas de temps 


III.8.4.5 Application du test d’adéquation Kolmogorov-Smirnov 


Afin de déterminer la loi la mieux adaptée pour l’ajustement statistique des intensités 


maximales observées, nous allons baser sur le test de Kolmogorov-Smimov (KS) 


(Kolmogorov, 1933 ; Smirnov, 1939). C’est parmi les méthodes non paramétriques les plus 


utiles et les plus générales pour comparer deux échantillons afin de déterminer s'ils suivent 


la même distribution. 


Le test KS est un test sans distribution, qui est basé sur l'examen de la distance verticale 


maximale entre les fonctions de distribution empiriques de deux échantillons (Figure 


111.12). 
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Soit nı et 2 les tailles des deux échantillons, n1> n2. 
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Figure IIL.12 Illustration de la statistique de Kolmogorov-Smirnov 


La valeur de la statistique du test pour le test de Kolmogorov-Smirnov a deux faces et a 


deux échantillons est : 


T = Max,|F, (x) — F (x)| 


(111.21) 


L'idée est de calculer la distance maximale (en noir sur la figure IIL.12) entre les fonctions 


théoriques et empiriques. Si cette distance est inférieure de la valeur indiquée sur le tableau 


de Kolmogorov-Smirnov, on dira que l'échantillon peut suivre la fonction empirique, sinon 


on dira que la loi n’est pas adéquate. 


Tableau IIL.7. Résultats du test de Kolmogorov-Smirnov 


Test de Kolmogorov-Smirnov 
Ecart max toléré= 0,2872 (N=21, IC=95%) 


Ecrart d'érreur AF 


Loi de Loi de Log- 
Gumbel Gumbel Loi des GEV 
3h 0,0998 0,0705 0,0659 
6h 0,2200 0,1639 0,1168 
9h 0,1781 0,1270 0,1342 
12h 0,1468 0,1600 0,1488 
15h 0,1132 0,1289 0,1079 
18h 0,1160 0,1090 0,0824 
21h 0,1113 0,0741 0,0624 
24h 0,1172 0,0737 0,0692 
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D’après les résultats indiqués sur le tableau, on déduit que l’échantillon suit les trois lois 
(Gumbel, Log-Gumbel et GEV). La différence ou bien l’erreur minimale est indiquée sur 
le tableau en gras. La loi GEV paraît à priori la mieux adaptée pour l’ajustement statistique 
des échantillons des fortes valeurs des intensités observées durant la période étudiée et 


pour les différentes durées d'observation. 


IIL.8.4.6 Prévision des intensités pour différents temps de retour 
L’estimation des lames précipitées de différents temps de retour T est basé sur loi GEV en 


procédant les étapes suivantes : 


- Calcul de la fréquence de non-dépassement d’après la relation (III.11) : 
1 
EP) =1-7 (III. 22) 


- Calcul de la variable réduite de la loi GEV correspondante d’aprés la relation 
suivante : 
u = —In(— In(F,(P)) (111.23) 
- Calcul du quantile correspondant d’après la relation citée précédemment (Voir 


l’équation HI.20) : 
S 
X = Xo + ral —(-—In(F,(P))*)  (II.24) 


L’ensemble des résultats pour tous les intervalles de temps (t=3h a 24h) sont regroupés 
dans le tableau suivant : 


Tableau IILS8 Estimation des lames précipitées pour différents temps de retour 


T=2ans T=5ans T=10ans T=20ans T =50ans T =100 ans 


Durée de la Pluie Pluie Pluie Pluie Pluie Pluie 
pluie [h] (mm) (mm) (mm) (mm) (mm) (mm) 
3 21,4 30,3 37,4 45,3 57,3 67,9 

6 25,4 35,0 44,4 56,7 79,4 103,3 

9 28,6 39,7 50,3 63,7 87,5 111,8 

12 31,5 43,7 54,0 66,1 85,6 103,7 

15 33,1 46,5 57,2 69,2 87,6 103,9 

18 34,2 48,2 59,2 71,2 89,3 104,9 

21 34,6 48,2 59,2 71,7 91,2 108,7 

24 35,1 49,0 60,0 72,3 91,2 107,9 


La détermination des intensités pluviométriques prévues se fait en calculant le rapport 


entre la quantité totale de pluie (hauteur des précipitations) tombant pendant une période 
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donnée et la durée de cette période (Tableau III.9). Elle est exprimée en unités de hauteur 
par unité de temps, généralement en mm par heure (mm/h). Elle est exprimée par la 


relation suivante : 
i=- (111.25) 
Où : 
i: intensité moyenne de la pluie [mm/h] ou ramenée à la surface [1/s.ha], 


h : hauteur de pluie de l'averse [mm], 


t: durée de l'averse [h]. 


Tableau IIL.9. Estimation des intensités horaires pour différents temps de retour 


T=2ans T=5ans T=10ans T=20ans T =50ans T =100 ans 


Intensité Intensité Intensité Intensité Intensité Intensité 
Durée de la moyenne moyenne moyenne moyenne moyenne moyenne 
pluie [h] (mm/h) (mm/h) (mm/h) (mm/h) (mm/h) (mm/h) 
3 7,1 10,1 12,5 15,1 19,1 22,6 
6 4,2 5,8 7,4 9,5 13,2 17,2 
9 3,2 4,4 5,6 7,1 9,7 12,4 
12 2,6 3,6 4,5 5,5 7,1 8,6 
15 2,2 3,1 3,8 4,6 5,8 6,9 
18 1,9 24 3,3 4,0 5,0 5,8 
21 1,6 2,3 2,8 3,4 4,3 5,2 
24 1,5 2,0 2,5 3,0 3,8 4,5 


III.8.4.7 Estimation des paramètres de la loi de Montana et représentation des 
courbes IDF: 

De nombreuses formes d'équations sont utilisées pour décrire les relations entre l'intensité 
et la durée des précipitations. On cite ci-dessous la loi de Bernard (1932) ou également 


appelée la loi de Montana, l’équation la plus simplifiée à utiliser : 


a 


er: 


(III. 26) 


Où / est l'intensité des précipitations (mm/h), t est la durée de la tempête (min), et a et b 
sont des paramètres liés aux conditions météorologiques. 

L’estimation des paramètres a et b de Montana est encore simplifiée en prenant le 
logarithme de cette formule de manière à obtenir une relation linéaire : 


In(Ir(t)) = In(a) — b. In(t) (III. 27) 
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Pour chaque temps de retour T, la droite de régression passant par les couples In(t), 


In(Ip (t)) permet d’estimer les paramètres de Montana (Figure III.13) 


3.500 
y = -0.8237x + 4.1826 
3.000 | y=-0.8016x + 3.9359 x 
y = -0.7768x + 3.6157 
2.500 |y-—-0:7623x-+-3-3772 
y = -0.7525x + 3.13848 
g 2000 | -y=-0.7504x + 2.7995 
g 
= 1.500 
1.000 
0.500 
0.000 
0.000 0.500 1.000 1.500 2.000 2.500 3.000 3.500 


Ln (T) 
O- T= 2ans = O- T= Sans ©- T= 10ans A= T= 20ans x T= 50ans © T= 100ans 


Figure IIL.13 Droites de régression de Y en X calculées sur les couples In(t), In(Ir (t)) 


Le calcul des paramètres a et b de Montana se fait en prenant encore l’exponentielle de 


In(a) et le paramètre b (i.e. - pente précédente). Les résultats sont regroupés dans le tableau 


suivant : 
Tableau IIL.10 Paramètres a et b de Montana 
Temps de retour (ans) 
T=2ans T=Sans T=10ans T=20ans T= 50ans T= 100ans 
b 0,75 0,75 0,76 0,78 0,80 0,82 
a 16,44 23.07 29,29 37,18 51,21 65,53 


La figure suivante II.14 représente les courbes IDF avec la formule de Montana 
déterminés à la station météorologique St5 de Bordj Bou-Arreridj pour la période 1990- 


2010 
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Intensité moyenne de l'averse [mm/h] 


0 ; 
Durée (h) 


—©— Intensité (T=2ans) —+— l Intensité (T=5ans) —A- Intensité (T=10ans) 
—*- Intensité (T=20ans) —#— Intensité (T=50ans) —O— Intensité (T=100ans) 


Figure IIL.14 Courbes IDF (Intensité-Durée-Fréquence) d la Station de Bordj Bou-Arreridj 
(Période 1990- 2010) 


Il ressort de la figure que plus la durée augmente, plus l’intensité diminue. Les averses les 
plus brèves (3h dans notre cas) sont souvent les plus violentes, au contraire les longues 


pluies (24h dans notre cas), sont généralement assez faible (de faible intensité). 


Selon la loi de Montana (Estorge et al, 1980 ; Weisse & Bois, 2002), les pluies de courte 


durée de fréquence donnée P(t) et P(t’), correspondant à deux pas de temps différents, (t) 


i A , P(t) A 
et (t ), pour la même variable, sont rapportées par : PU — (=) (IT. 28); 
tP 
Prenant t' = 24 h, l'équation devient : P(t) = Poan max (=) (IIT. 29) 


P(t): Pluies de courte durée de fréquence donnée, en heure ; Po4p max : Pluie journalière 


maximale de fréquence donnée ; b : exposant climatique de la région 


Dans le cas de la station St5, on a estimé un exposant climatique b=0.75. Cette valeur est 


trés proche de la valeur moyenne trouvée par (Laribi, 2012). 
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Sinon l’exposant climatique estimé reste très élévé par rapport à d’autres études où b 
atteint 0.42 (ANRH) ou de celles publiés par d’autres auteurs (Boulghobra, 2013 ; Chenafi, 
2013: Bouaïchi et al, 2014). La différence peut être interprété par la variation spatiale de 
l'intensité, et aussi par la variation des paramètres spatiales de la station pluviométrique 
comme la longitude, la latitude et l’altitude X,Y,Z (Habibi et al, 2013 ; Guidoum et al, 
2014). 


77 


Chapitre IV 


Etude du transport 
des sédiments en suspension 


Chapitre IV Etude du transport des sédiments en suspension 


Chapitre IV: Etude du transport des sédiments en suspension 


IV.1 Introduction 
Le transport de sédiments dans les rivières est un élément clé dans l'étude des cycles globaux de 
l'eau (Megnounif et al, 2013). La connaissance de la quantité de sédiments transportés est 
importante pour de nombreuses raisons. Il s'agit notamment des impacts sur la fertilité des sols, 
le transfert, le stockage et le devenir des nutriments et des contaminants, les changements dans 
les tendances de la qualité de l'eau, les habitats aquatiques, l'envasement des canaux, des 
réservoirs et des ports, et la réduction de la longévité des équipements hydroélectriques 
(Williams, 1989 ; Ouillon, 1998 ; Horowitz, 2003 ; Achite et Ouillon, 2007). 
En outre, l'étude des rivières est nécessaire pour des prévisions fiables, mais c'est une tâche 
difficile en raison de la complexité des systèmes hydrologiques (Leahy et al, 2008 ; Khanchoul 
et al, 2012). La surveillance du transport de sédiments nécessite une bonne technique 
d'échantillonnage, qui est très longue et coûteuse. Par conséquent, il est nécessaire de 
développer des modèles qui puissent prédire avec précision les charges de sédiments en 
suspension à partir d'un ensemble de données continues sur l'eau. 
Un des modèles les plus courants, comportant moins de paramètres physiques, est la courbe 
d’estimation de la charge solide en suspension, ou appelé «sediment rating curve», qui 
représente une courbe fonctionnelle de puissance, relie-la concentration/charge de sédiments en 
suspension au débit du cours d'eau sont un type principal (Walling, 1974; Walling, 1978 ; 
Asselman, 2000). 
Les objectifs du présent chapitre peuvent étre résumés comme suit : 

- Présenter quelques données et chiffres sur la production des sédiments dans le monde, 

dans le grand Maghreb et en Algérie. 
- Passer en revue, les différentes méthodes empiriques d'estimation du transport solide en 


Algérie. 
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-  Quantifier la production des sédiments dans le bassin d’étude, Oued K’sob, en se basant 
sur les données instantanées disponibles (débits liquides, concentrations) pour une 
période de 20 ans. 

- Développer un modèle puissance reliant la charge sédimentaire au débit liquide, qui 


pourra être utilisé en cas d’absence des données. 


IV.2 Quelques chiffres sur l’érosion hydrique 

IV.2.1 La dégradation du sol dans le monde 

A l’échelle du globe, les bilans établis, surtout à partir de 1950, révèlent que le tonnage annuel 
des sédiments en suspension exportés vers les océans, à partir des continents, est très variable 
(Megnounif 2007). La figure I montre à l’échelle mondiale, la variabilité spatiale du transport 
des sédiments des rivières vers les océans. 


—150 —100 —50 O 50 100 150 


SFr (t km yr™') |» 
EM Above 1000 
EM 500 - 1000 
EM 100 - 500 
E 50- 100 
m] 19- 50 
| 1- 10 
EM Baw 


—150 —100 =50 0 50 100 150 


Figure IV.1 Bilan mondial des transferts fluviaux des sédiments vers les océans (Ludwig et 


Probst, 1998) 


Dans une étude menée sur plus de 200 rivières, Millimane et Syvistiky (1992) ont évalué une 


dégradation de sol de 225 t km? an‘! 
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En Amérique latine, une perte de 123 t km? an‘! a été estimé au Mississippi, et 203 t km? an’! 
en Amazone (Ludwig et Probst, 1998). Selon Shi et Zhang (2005), le fleuve jaune situé en Asie 
(Chine), transporte environ de 1.837x10!! tonnes de sédiments. En Europe, la Garonne en 
France, transporte à peu près 27.6 t km? an‘! de matières en suspension (Probst et Bazerbachi, 
1986). Au niveau du Danube, qui est le deuxième fleuve d'Europe par sa longueur, environ de 
107 t km? an‘! de sédiments transportés ont été estimés par Ludwig et Probst (1998). 

En afrique, le Sénégal (un fleuve d'Afrique de l'Ouest) transmet 8.75 t km? an! de MES, tandis 
que le fleuve Nile transporte 64 t km? an! (Ludwig et Probst, 1998). 


IV.2.2 La dégradation du sol au Maghreb et en Algérie 

Au cours des années soixante, les études menées sur l’érosion dans le grand Maghreb, ont 
souvent négligé l'importance de l'érosion au Maghreb. Selon Fournier (1960) et Strakhov 
(1967), le Maghreb est considéré comme l'une des régions les moins érodées après les zones 
désertiques. Milliman et Meade (1983) ont supposé que la production des sédiments dans 
les rivières maghrébines était d'environ 100 t km? ant. Walling (1984) a suggéré que toutes ces 
valeurs étaient sous-estimées et a proposé un rendement sédimentaire en suspension de 1000 à 
5000 t km”? an'!. Cette estimation était basée sur plusieurs études des fleuves Maghrébins dont 
celles de Heusch & Milliès-Lacroix (1971) au Maroc, en Algérie et en Tunisie, Demmak (1982) 
en Algérie, Ghorbel & Claude (1977) en Tunisie et quelques autres. 

Snoussi (1988) et Snoussi et al. (1990) ont étudié trois grands fleuves au Maroc et ont estimé le 
rendement annuel moyen en sédiments en suspension à environ 750 t km? an‘! pour le Maroc. 
La disparité entre les différentes valeurs citées ci-dessus reflète à la fois la difficulté de mesurer 
l'érosion mécanique et le transport fluvial des sédiments en suspension dans ces régions semi- 
arides et la variabilité spatiale de l'intensité de l'érosion en réponse au climat très contrasté et 


agressif, aux sols mal végétalisés et au fort relief (Probest et Suchet, 1992). 


IV.2.3 La dégradation du sol en Algérie 


En raison de son climat semi-aride, de la présence de relief abrupt et de la pauvreté de la 
végétation, l’Algérie est l’une des régions les plus vulnérables au monde à l’érosion des sols 


(Megnounif, 2007). Le nord de ce pays est une région montagneuse caractérisée par un petit et 
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moyen bassin de moins de 10 000 km? avec une forte variabilité spatio-temporelle de l'intensité 
de l'érosion en réponse à une végétation clairsemée et un climat irrégulier et agressif (Ghenim et 


Megnounif 2016). 


L’intensité de l’érosion hydrique varie d’une zone à l’autre. La partie Ouest du pays est la plus 
érodée, où l'érosion touche 47 % de l’ensemble des terres ; suivie du Centre (27%) et de l’Est 
(26%). L'érosion spécifique varie de 2000 et 4000 t/km?.an et le taux d’envasement est 
supérieur à 15 %. L'Algérie est, de ce fait l'un des pays les plus menacés dans le monde par 
l'érosion (Achite et al, 2006). Le tableau suivant montre l’ordre de grandeur de l’érosion 


hydrique dans les bassins versant Algériens. 


Tableau IV.1 Estimations de la production de sédiments dans les bassins versant algériens 


Auteurs Bassin versant Periode S (km?) [Ass 
(10°kg/km/an) 
Bourouba, 1998 Oued Medjerda 1975-1985 217 51.5-242.5 
Meddi, 1999 Oued Ebda - - 1875 
Terfous et al. 2001 Oued Mouilah 1977-1993 2650 {126 
Touaibia et al, 2001 Oued Mina 1973-1995 6200 + 
Amireche, 2001 Zardézas 1977-1986 345 321 
Megnounif et al. 2003 Haute Tafna 1988-1993 256 24 - 4288 
Terfous et al. 2003 Haute Tafna 1989-1998 1255.5 1120-6000 
Isser 1989-1998 1140 4-3400 
Mouillah 1986-1995 2650 22-980 
Boudjadja et al. 2003 Oued Allalah 295 2701 
Oued Damous 577 2879 
Oued Es Sebt - 112 2950 
Oued Messelmoun - 218 3029 
Oued Elhachem - 217 2905 
Achite et Meddi 2004 Oued Haddad 1973 - 1995 70 287 
Oued ElAbtal 1973 - 1995 4126 
Sidi A.E.K. Djilali 1973 - 1995 470 212 
Achite et Meddi 2005 Aïn Hamra 1973 -1995 2480 |117 
Kef Mahboula 1973 -1995 1680 - 
Takhmaret 1973 -1995 11553 {65 
Oued Mina 1973 -1995 4900 |187 
Larfi, et al. 2006 Oued Isser 1984-1997 B615 {1100 
Ghenim et al. 2007 Oued Sebdou 1985-1998 [256 1330 
Achite et Ouillon, 2007 Oued Abd 1973-1995 2480 1136 
Khanchoul et Jansson. 2007 Oued Saf Saf 22 ans B22 461 
Oued Kebir West - 1130 247 
Khanchoul et Jansson. 2008 Oued Mellah 1972 — 1997 1550 373 
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Meguenni et Remini. 2008 (Oued Harazza 1973 — 1997 1142 470 
Ghenim et al, 2008 Oued Mouilah 1977- 1995 7245 16.4-1038 
Cherif et al, 2009 Oued Mekerra 1950-2001 4102 114-296 
Hasbaia et al. 2012 Oued El Ham 1968-1989 5605 530 
Yles et Bouanani. 2012 Oued Saida 1980-2010 1543 55-29667 
Khanchoul et al. 2012 Oued Kebir 1976 —2008 (681 895 

El Mahi et al. 2012 Oued El Hammam 1993 —2005 8348 256 
Cherif et al. 2012 Oued Taria 1987 -2007 11365 236 
Elahcene et al. 2013 Oued Bellah 1974-2007 55 610 
Louamri et al. 2013 Oued Boumahdane 11969 -2010 11105 257 
Bouanani et al. 2013 Oued Sikak 1972 - 1990 218 170 
Megnounif et al. 2013 Oued Sebdou 1973 — 2004 256 343 
Boukhrissa et al., 2013 El Kebir 1975-1999 (681 1410 
IBenkadja et al .2015 Oued el Ksob 1973 — 1991 |1484 208 
Bouguerra et al. 2016 Oued Boumessoud |1988 — 2004 1118 518 
Tadrist et Debauche, 2018 Boukadir 1993-2013 + 446 
Selmi et Khanchoul, 2016 Mellegue 1970 — 2003 4575 589,23 
Cherif et al. 2017 Oued El Hammam 1986-2008 [7440 {111 
Hasbaia et al. 2017 Oued Soubella 1974-1989 183.5 1126 
Tani et al 2017 Oued Nachef 1990-2010 + 153 
Hallouz et al, 2017 Oued Mina 1971 —2002 + 211 
Tourki et al. 2017 Oued Kebir amont 1973 —2006 11068 884 
Madani Cherif et al, 2017 Oued Mekerra 1987-2007 + 38 
Benselama et al. 2018 Oued El Maleh 1981 — 1998 (932.5 [294.29 


IV.3 Méthodes empiriques d'estimation du transport solide 


La quantification de l'érosion est généralement difficile et dépend des échelles spatio- 
temporelles. Elle est mesurée directement sur le terrain ou elle est estimée suivant des analyses 
des sols ou des modèles empiriques ou physiques qui tiennent compte de l'impact de toutes les 
variables de l’érosion hydrique. 

Nombreuses formules et méthodes qui impliquent d’autres paramètres physiques ont été 


adoptées pour la quantification de l'érosion hydrique : 


= Formule de Henin (1950) : 
Elle symbolise la relation existante entre l'érosion et les facteurs susceptibles de la favoriser ou 
de la limiter. La quantification de l'érosion spécifique "Es" est donnée par la formule IV.1 
mer 


E 
AY KV, 


(IV.1) 
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Avec, I, : Intensité des précipitations ; 

I, : Pente du bassin versant ; 

Les facteurs "J," et "I," favorisent l'érosion et S : Susceptibilité des sols ; 

K : Perméabilité ; 

V : Végétation. 

Les facteurs S, K et V, limitent l'érosion et en plus, les facteurs S et V, semblent les plus 


difficiles à quantifier. 


= Formule de Gravilovic (1960) 


Cette formule lie l’érosion hydrique à la température, à la pluie et à un coefficient d’érosion, 


l’équation est donnée par l’expression IV.2: 


E, = 3.14T. Py VZ? (IV. 2) 


Ou 
E, : érosion spécifique dans le bassin versant (T/km2.an);- 
T: coefficient de température; 


r= [+01 IV.3 
~ 11007 ee) 


Avec: 

to : température moyenne annuelle en °C. 
Py : pluie annuelle en (mm). 

Z : coefficient d’érosion. 


Avec : 


Z = y. [x.a]. (8 + Joy) (IV. 4) 
Où 


y: valeur du coefficient de la résistance du sol dû à l’érosion. Il dépend de la roche mère, le type 


du sol et du climat. 
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[x. a] : Coefficient de régularisation du bassin versant, se rapportant à la protection des sols, des 
influences des phénomènes atmosphériques des forces érosives liées aux conditions naturelles. 


Il est aussi tabulé. 
6 : équivalant numérique des processus visibles et nettement prononcés dans le bassin versant. 
Imoy: indice de pente moyenne du bassin versant. 


= Formule de Tixeront (1960) 
Basé sur les données de 32 bassins d’Algérie et 9 bassins de la Tunisie, la formule de 


Tixeront (1960) s’écrit comme suivant : 


Es = [Si Tor + S2. Taz + S3- Taz +... +Sn. Tan] (V5) 


E,: Erosion spécifique (t/Km2.an) ; 

Tar Taz + Tan: Taux d’abrasion en fonction de la perméabilité du sol dans les sous Bassin 
(T/Km’) ; 

S: Superficie totale du bassin (Km?) ; 

S1, S2. - - Sn: Superficie des sous bassins (Km?) ; 

Maintenant, si on suppose que la totalité du bassin possède une même perméabilité, 

Ona E,=K.Le (1V.6) 


Le: Lame d’eau ruisselée moyenne interannuelle (mm). 


K : coefficient qui tient compte de la perméabilité des terrains. 
K = 8,5 pour les terrains d’une perméabilité élevée. 
K = 75 pour les terrains d’une perméabilité moyenne élevée. 
K = 350 pour les terrains d’une perméabilité faible à moyenne. K = 1400 pour les terrains 
d’une perméabilité faible. 


K = 3200 pour les terrains perméables. 


= Formule de Fournier (1960) 
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1 [PS H2, 0.46 
Ass ==. =|) |= IV.7 
Ass  : apport solide moyen inter annuel (T/Km?.an). 


P : la pluviométrie du mois le plus arrosé (mm). 
Hmoy :Ľľaltitude moyenne (mm). 
Pa : la pluviométrie moyenne inter annuelle (mm). 


S : la superficie du bassin versant (Km?). 


=" Formule de SOGREAH : 
SOGREAH (1983) a introduit la formule suivante, en se basant sur les données de 30 bassins 
versants algériens : 
Ass = aE%15 (IV.8) 
Avec : 
À, : Apport solide spécifique (T/km?/an) ; 
a : Coefficient dépendant de la perméabilité du sol ; 


E : Ecoulement annuel (mm) ; 


= Formule de L’A.N.R.H (1970) : 
Ass = 26.62 Ie + 5.0711, + 9.77 C; — 593.59 (1V.9) 

Où : 
Ass: Apport solide spécifique moyen annuel (T/Km2.an) ; 
Ie: Indice lithologique (%) ; 
I, : Indice des précipitations (%) ; 
Cı : Coefficient de torrentialité : C, = Dg. Fi (1V.10) 
Dd : Densité de drainage : 


— Sri Lı 
q A 


L; : longueur de tous les cours d’eau d’ordre « 1 » (Km). 


(IV.11) 
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F,= nombre de talweg/ surface fréquence des talwegs élémentaire Cette formule est applicable 
dans les conditions suivantes : 

Ie> 10% 

I> 300 


= Formule de Probst et Amiotte-Suchet (1992) 
Ils ont proposé une formule en utilisant les données de 130 bassins du Maghreb (coefficient de 


corrélation de ce modèle est égal à 0,71). 


Ln(A,,) = 4.79 + (0.054)Kzr + (0.004)R — (5.6 x 107°).A (1V.12) 
Où : 
Ass: Apport solide annuel moyen (T/Km*.an) ; 
Ker Érodibilité du sol ; 
R : Apport liquide annuel moyen ; 
À : Surface du bassin versant ; 
IV.4 Mesure et quantification du transports solide 

IV.4.1 Mesure des matières en suspension 
L'étude est basée sur les valeurs instantanées des débits liquides (Q en m3/s) et des 
concentrations en matières en suspension (C en g/l) mesurées entre Juin 1973 et Janvier 1992 à 
la station hydrométrique Medjez de cordonnées (X=35°53'25.78", Y=4°037'9.27", Z= 636 m). 
La localisation et le profil en travers de la station sont montrés dans les figures III.1.2. Les 
données hydrométriques utilisées sont mesurées et fournies par les services de l’ ANRH. 


(Agence Nationale des Ressources Hydraulique). 
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N Légende 


a Station hydrométrique 


Réseau hydrographique 


Er Barrage K'sob 


0 4 8 16 Kilometers 


Figure IV.2 Localisation de station hydrométrique Medjez dans le bassin K’sob 


Le protocole pratiqué par les services de ANRH pour mesurer le débit liquide et la charge 
solide en suspension est le même partout en Algérie. Les hauteurs d’eau sont transformées en 
débits liquidés par une courbe de tarage établie localement. Ces courbes sont actualisées et 
vérifiées périodiquement. Pour l’estimation du débit solide, un échantillon d’eau est prélevé au 
moyen d’un flacon de 1L. La fréquence de prise d’échantillon est adaptée au régime 
d’écoulement. Elle est intensifiée en période de crue, jusqu’à une prise par demi-heure. En 
dehors des crues, un échantillon est pris 1 jour sur deux. Lorsque les écoulements sont très 
faibles, un échantillon est pris tous les quinze jours. A chaque prise d’échantillon, le surveillant 
lit la hauteur d’eau pour en déduire le débit correspondant. La concentration des matières 
solides en suspension est estimée au laboratoire par pesée sur filtre suivant un protocole détaillé 


dans les études (Achite and Ouillon 2007 ; Megnounif et al. 2013). 
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Figure IV.3: Profil en travers de la station hydrometrique Medjez, code 05.09.01, date du 03/03/2007 cote h=-18cm 
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IV.4.2 Quantification du transport solide 
L’apport solide A, exporté par l’oued K’sob sur une période donnée (jour, mois, saison, année) 
est calculé en sommant les flux massiques instantanés répartis uniformément entre deux instants 


de mesures : 


N 


a i N C;Q; + Cir1Q jrs = 
Ass | C(t)Q(t)dt = > (Eee (tu — tj) x 10 (IV.13) 


= 2 

=0 j =0 
Où Q; (en m3/s) et C; (en g L/s) désignent respectivement le débit et la concentration de 
sédiments, mesurés au moment t; où l'ensemble de l'événement de crue est divisé en N périodes 


correspondant à (N + 1) d’échantillonnages ; 


tj 


+1 — tj : temps séparant deux prélèvements consécutifs en secondes. 

À la station de jaugeage qui contrôle Wadi K'sob, une mesure continue du niveau de l'eau, 
déduite d'un limnigraphe, est convertie en débit à l'aide d'un abaque local résultant de la 
disponibilité de mesures continues des débits quotidiens moyens sur la période 1973-1994. 


L’année 1980-1981 est exclue du calcul à cause de la non disponibilité de ses données 


Les résultats de la concentration des sédiments en suspension proviennent d'échantillons 
prélevés manuellement. Outre certaines lacunes observées, on peut soupçonner que des 
échantillons n'ont pas été prélevés lors d'une brève augmentation des rejets la nuit ou pendant 
les vacances. Cela justifie l'utilisation de la technique de la courbe d’estimation de la charge 
sédimentaire qui est largement appliquée pour estimer les apports sédimentaires tant pour les 
petits que pour les grands bassins versants (Walling 1977 ; Asselman 2000 ; Horowitz 2003 
Warrick 2015 ; Megnounif et al. 2016). 


Sur la période 1973-1994, les débits solides journaliers ont été évalués par le modèle (IV.14) et 
utilisés pour estimer la production de sédiments A, pour une période donnée (Olias et al. 2006 ; 


Louamri et al. 2013) à l'aide de la formule suivante : 
n n 
As=m) Qsa =m) aQ” (IV.14) 
j=1 j=1 
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Où Qa et Qa sont respectivement le débit moyen journalier de l'eau et le débit solide moyen 
journalier, m est une constante de conversion d'unité et n est le nombre de jours pour une 


période de temps. 


IV.4.3 Courbes de transport solide (SRC) 
Depuis plusieurs décennies, et en raison de l'absence de mesures réelles reliées au transport des 
sédiments, continues ou quasi-continues, les hydrologues utilisent souvent des courbes 
d’estimation pour prévoir les concentrations des matières en suspension « MES » quotidiennes 
et pour le calcul des flux. Bien qu'il existe nombreuses méthodes pour élaborer des courbes 


d’estimation de la charge solide en suspension, 


La plus courante est une fonction de puissance ou également appelée « sediment rating curve, 
SRC ». Cette courbe décrit empiriquement la relation entre la concentration de sédiments en 
suspension C (g/m?) et le débit liquide Q (m*/s) pour un endroit spécifique (Walling, 1974, 
Walling, 1978, Fenn et al., 1985, Jansson, 1996, Mossa, 1996, Asselman, 1999, Asselman, 
2000, Syvitski et Morehead, 1999, Syvitski et al., 2000, Morehead et al., 2003), elle est 


exprimée comme suivant : 
C = aQ? (IV.15) 
Où: 
Les coefficients a et b sont les paramètres d'évaluation des sédiments ; 


Cela s'exprime également sous la forme d'une équation linéaire basée sur une transformation 


logarithmique : 


log (C) = log(a) + blog(Q) (1V.16) 


Selon Gilroy et al. (1990), le calcul des apports solides en utilisant le modèle puissance (IV.15) 
et en le multipliant ensuite par le débit liquide est équivalent au calcul des apports solides en 


utilisant le modèle : 


Q, = 407" (IV.17). 
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La courbe de cotation des sédiments peut être considérée comme un modèle de "boîte noire" et 
les paramètres a et b sont estimés par une analyse de régression sans signification physique 
(Asselman, 2000). Néanmoins, autres chercheurs ont associé ces paramètres à la morphologie 
du lit des rivières (par exemple, la forme du chenal, la pente et l’érosivité du cours d'eau) ou à 
l'érodabilité du sol et des sections de rivière (Peters-Kümmerly, 1973, Rannie, 1978, Fenn et al. 
1985, Thomas, 1988, Morgan, 1995, Asselman, 1999, Asselman, 2000, Syvitski et al. 2000, 
Morehead et al. 2003, Yang et al. 2007a, Wang et al. 2008). 


IV.5 Analyse des résultats et discussions 

IV.5.1 Relations débit liquide- débit solide 

IV.5.1.1 A l’échelle interannuelle 

L’analyse porte sur les valeurs des débits liquides en m3/s et des débits solides en kg/s, mesurés 
au niveau du bassin versant Oued K’sob de Juin 1980 a Janvier 1992 (données instantanées 


pour la série complète) (figure IV.1). 


1.E+05 
1.E+04 
1.E+03 
1.E+02 


1.E+01 


1.E+00 Qs = 5.697Q!38! 


R? = 0.8673 


Débit solide(g s-1) 


0.01 0.1 1 10 100 1000 
Débit liquide (m3 s-1) 


Figure IV.4 Relation entre les débits solides et les débits liquides instantanés de 1973/1974 a 
1991/1992 (Oued K’ osb - Station hydrométrique Medjez). 


La figure IV.4, montre une bonne corrélation (R?=0.86) en puissance liant le débit des matières 


solides en suspension aux débits liquides, les points du nuage obtenus sont alignés autour de la 
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droite de régression. La relation puissance [Q, = aQ?] a été vérifiée pour la plupart des cours 
d’eau dans le monde. L’exposant b, fonction des caractéristiques physiques, climatiques et 
hydrologiques des bassins versants (Probst et Bazerbachi, (1986) et Vivian, (1980) ou des 
conditions hydrauliques de l’écoulement dans les cours d’eau (Kattan et Probst, 1987). On 
notera aussi que l’exposant b, de notre modèle est égale à 1.38, cette valeur est proche de celle 
trouvée par Achite et Meddi (2005) et Cherif et al (2009), et par bon nombre de chercheurs 


ayant travaillé dans des régions semi-arides (tableau IV.2). 


Tableau IV.2 Quelques valeurs de « b » publiés pour des bassins versant voisins et dans le 


monde. 

Bassin versant L’exposant «b» | Auteurs 
Oued Mekerra 1.31 Cherif et al, 2009 
Oued Mina 1,33 Achite et Meddi, 2005 
Oued K’sob 1,38 Guesri et al, 2020 

2 (Présente étude) 

$ Oued Haddad 1,45 Achite et Meddi, 2004 

ia Sikkak 1,47 Bouanani et al.2013 
Haute Tafna 1,57 Megnounif, 2007 
Tafna 1,70 Terfous et al., 2001 
Rivière Arbucies (Nort- | 1,63 Battala et Sala, 1992 

© o| Est) Péninsule Ibérique 

É 3 La Garonne (France) 1a2 Probst et al, 1986 

Lg a 

| —| Pratts-ville (U.S.A) 2,00 Leopold et al, 1953 


On observe aussi d’après la figure IV.4, que pour des faibles débits liquides (< 0.1 m3/s) 
peuvent étre associés de forts débits solides (prés de 0,1 Kg/s) déclenchés probablement a la 


suite d’orages violents ou de remobilisation de sédiments issus d’une crue précédente et 
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fraîchement déposés dans le lit de l’oued. De même, de fortes valeurs des débits liquides n’ont 
engendré qu’un faible transport solide dû probablement à la vidange des nappes survenue à la 


suite d’une crue importante. 


IV.5.1.2 À l’échelle saisonnière 

Afin d’étudier les réponses du bassin aux débits liquides et matières en suspension au cours de 
l’année hydrologique, nous avons jugé utile de regrouper suivant les différentes saisons les 
valeurs instantanées (prélevées sur les 20 années d’étude), et d’analyser la relation liant les 
débits liquides aux débits solides en suspension. Le tableau IV.3 résume les modèles trouvés 


pour les différentes saisons. 


Tableau IV.3 Relations débit solide-débit liquide pour chaque échelle de temps (1973/1974 à 
1991/1992). 


Echelle Nombre Coefficient de Relations retenues 
temporelle d’observations determination (R?)/% 

Série complète [1646 0,86 (86%) Q,=5.697Q"%! 
Automne 580 0,92 (92%) Q,=8.066Q'5 
Hiver 341 0,85 (85%) Q.=3.779Q"*" 
Printemps 478 0,82 (82%) Q,=3.655Q"%* 

re 247 0.92 (92%) Q,=9.485Q 1" 


Les valeurs des coefficients de déterminations trouvés indiquent la forte liaison entre les débits 
solides et les débits liquides pour l’ensemble des saisons. La représentation graphique des 


résultats (figure IV.5 (a), (b), (c), (d)) permet les commentaires suivants : 
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1.E+05 


a) Automne b) Hiver 
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Figure IV.5 Modèles saisonniers reliant les débits solides aux débits liquides instantanés 


1973/1974 à 1991/1992 (Oued K’sob - Station hydrométrique Medjez). 


En Automne (Figure IV.5.(a)), les paramètres a et b sont élevés par rapport aux autres 
saisons, avec une valeur de b qui est égale à (1,435). Ceci peut témoigner de l’importance de 
la dynamique érosive en cette saison. A l’instar des régions semi-arides en automne, la région 
de Oued K’sob est caractérisée par des orages violents qui surviennent en septembre, octobre, 
et novembre, qui arrachent de grandes quantités de matières solides, véhiculées en suspension 
par le cours d’eau. L’automne se distingue par un flux important de matières solides en 
suspension. Le débit solide maximal est de l’ordre de 57.62 kg/s, atteint en fin de Septembre 
1989, associé à un débit liquide de 109,9 m’/s. La relation entre le débit liquide et solide 
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semble être plus étroite pour les valeurs élevées. 


En Hiver (Figure IV.5.(b)), le modèle obtenu dénote une diminution de b pour atteindre la 
valeur de 1,361. Ceci est dû essentiellement aux réserves importantes en eau emmagasinées 
dans le sol, suite aux forts orages de l’automne. La charge en suspension n’augmente pas 
proportionnellement aux apports liquides si bien que le flux saisonnier reste inférieur à celui 
mesuré en automne. On enregistre durant cette saison de forts débits liquides associés à de 
faibles débits solides. Le débit solide maximal enregistré est de l’ordre de 5.66 kg/s résultant 
d’un débit liquide de 10,53 m’/s. En effet, durant cette saison, l’état humide du sol et son 
occupation (période végétative) sont responsables de la faiblesse des apports solides. 
L’évaporation est faible, l’eau interceptée par la végétation est en grande partie restituée au 
sol par égouttage. Cette partie des précipitations peut donc atteindre le sol après un certain 
temps et la quantité d’eau perdue est négligeable. Dans ce cas, le ruissellement est retardé, ce 
qui donne un pouvoir érosif limité, en outre la réhumectation du sol par les pluies d’automne 


réduit considérablement la vulnérabilité des particules à l’érosion hydrique. 


Au Printemps (Figure IV.5.(c)), le coefficient a= 3.655 et l’exposant b =1,364, présente des 
valeurs faibles par rapport a celles de l’automne et celles de l’été. Le modèle du printemps 
pourrait bien être un modèle représentatif moyen. Au printemps, l’état du sol et la couverture 


végétale font que l’action érosive se trouve diminuée par rapport aux écoulements abondants. 


En Eté (Figure IV.5.(d)), l’exposant b reste élevé (1,373) après celui de l’automne. 
En effet, l’oued K’sob a connu ces dernières années quelques crues estivales très importantes, 
notamment la crue du 13 juin 1974 avec un débit moyen journalier de 18.7 m°/s et un apport 
solide de As = 2332 tonnes. Le coefficient b est lié directement à la superficie de drainage en 
période de crue, mais il exprime aussi l’importance du taux des transports solides. 
Le paramètre a exprime l’état du sol et son degré de saturation, la valeur pour notre modèle 
est de 9.485. Pendant l’été la pluviométrie est presque rare par rapport aux autres saisons, 


mais elle apporte des forts débits liquides et solides. 
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IV.5.2 Analyses des apports saisonniers. 
Les modèles saisonniers retenues, ont été utilisés pour calculer les apports solides moyens 
saisonniers en suspension. Les résultats trouvés sont donnés dans le tableau IV.4. et la figure 


IV.6 (a,b) 


Tableau IV.4 Répartition saisonnière des apports liquides Al (en Hm3), des apports solides As 
(en tonnes) et des dégradations spécifiques Ass (en tonnes par Km2 et par saison) à l’oued 


K’sob pour la période 1973/1974 à 1991/1992) (Station hydrométrique Medjez) . 


Saisons Automne! Hiver | Printemps Eté Total 
Al (Hm?) 122.02 110.09 117.52 49.04 398.66 
As (MT) 3.42 0.83 1.14 1.42 6.81 
Ass 115.46 27.96 38.65 47.96 230.03 
(t/km?/an) 
a) Répartition saisonnière b) Répartition saisonnière 
des apports liquides des apports solides 


E Automne MBÉté CHiver m Printemps E Automne MEté CHiver m Printemps 
Figure IV.6. Répartition saisonnière des apports liquides et des apports solides à l’oued K’sob 
(Juin 1973 - Janvier 1992) (Station hydrométrique Medjez). 
L’analyse de ces résultats, permet de révéler l’importance des apports solides durant la saison 


d’automne entre 1973-1992, avec un apport solide de 3.42 Million de tonnes, soit presque 50% 
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des apports solides annuels (Figure IV.6). Les fortes précipitations sur un sol à l’état sec après 
une longue saison d’été sont à l’origine de ce taux important. En effet, les premières pluies de 
septembre et d’octobre sont capables d’arracher de grandes quantités de matières solides, qui 
seront transportées dans les cours d’eau. 

En hiver, durant les mois de janvier, février et mars on enregistre de forts apports liquides 
110.09 Hm3 avec un apport solide de 0.83 Million de tonnes de matières en suspension ce qui 
représente seulement 12% par rapport au total annuel. L’hiver présente des caractéristiques 
érosives plutôt moins prononcées. Le sol ayant été saturé après les pluies d’automne résiste plus 
à l’érosion malgré des précipitations relativement importantes. La dégradation spécifique durant 
cette saison reste relativement faible 27.96 t/ha/an. 

Au printemps, la charge en suspension n’augmente pas proportionnellement aux apports 
pluviométriques, sur les 32% de l’apport annuel en eau on remarque que seulement 17% de 
matières solides annuelles seront arrachées essentiellement du lit et des berges de l’oued (Figure 
IV.6). La présence de la végétation et l’humidification des sols réduisent considérablement la 
migration des particules par l’érosion hydrique durant cette saison. 

Mise à part la saison d’automne, l’été se distingue des autres saisons par des valeurs très fortes 
de transport solide, environ 14% des apports liquides totales, amène environ de 21% d’apport 
solide annuel avec une dégradation spécifique de l’ordre de 47.96 t/ha/saison. Durant cette 
période les précipitations sont très agressives, et la susceptibilité du sol est haute, ce qui 


augmente considérablement l’érosion (figure IV.6 (a), (b)). 
IV.5.3 Bilan des apports solides annuels 


Les variations interannuelles des apports solides montrent une grande irrégularité, liée à la 


variation et fréquence des pluies et donc des apports liquides (figure IV.4). 
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Figure IV.7 Apports liquides et solides moyens annuels à l’oued K’sob 
(Période de 1973/1974 à 1991/1992). 


L’analyse de cette figure, montre que l’apport solide annuel calculé durant l’année 1975/1976 
est estimé à 1.05 million de tonnes. C’est un apport solide relativement important, malgré la 
faiblesse des apports liquides 2.73 Hm3. Les crues de l’automne (09/10-sept.-1975) et du 
printemps (27/28-juin-1976) ont apporté respectivement 598267 tonnes et 123944 tonnes de 
matières en suspension, ce qui représente 56 et 11% de l’apport moyen annuel. 

On observe aussi que pour deux années recevant presque la même quantité de pluie la charge en 
suspension peut être différente : durant l’année 1986/1987, on a estimé 0.49 million de tonnes 
pour une précipitation de 311 mm alors qu’en 1981/1982, pour presque la même quantité de 
précipitations 312 mm il n’a été enregistré que 0.11 million de tonnes de sédiments, soit un 


rapport de 1 à 4. 


Il faut aussi rappeler qu’on peut déterminer la dégradation spécifique (en t/km?/an) à partir de la 


relation suivante : 
A, 


Ou : S est la surface du bassin versant ; A, est l’apport solide annuel (t) . 
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On trouve alors que la dégradation spécifique est très variable d’une année à une autre. Elle 
varie de 18 t/km2/an pour l’année 1983/1984 à 718 t/km2/an pour l’année 1975/1976, soit un 
rapport de 1 à 13. 

Sur les 20 ans de mesure, L’apport moyen annuel des sédiments enregistré à l’exutoire du 
bassin versant de l’oued K’sob est évalué à 340444 tonnes, ce qui correspond à une dégradation 


spécifique de 230 t/km2/an. 


A titre de comparaison, cette estimation de dégradation de sol reste très modérée par rapport à 
celles publiées pour d’autres bassins de l’Algérie et du Maghreb (Tableau IV.1). On peut citer 
par exemple des estimations très proches, 236 et 247 t/km?/an évaluées à Oued Taria et oued 
Kbir respectivement (Cherif et al. 2012 ; Khanchoul et Jansson. 2007). Hallouz et al (2017) ont 
estimé 211 t/km?/an à Oued Mina. Dans l’oued Oued Boumahdane, Louamri et al (2013) a 


trouvé un taux de dégradation de 257 t/km?/an 


En outre, la valeur estimée à oued K’sob est relativement grande par rapport à des bassins 
versants à climats similaires (Tableau IV.1), on peut mentionner, l’oued Sikkak et Isser (Tafna) 
avec respectivement 169 et 180 t/Km2/an (Bouanani et al, 2013), l’oued Leham qui fait partie 
du même grand bassin Hodna, 104 t/Km2/an (Bourouba, 1998). 


Sinon, des travaux réalisés au Maghreb ont estimé que les dégradations spécifiques varient entre 
1000 et 5000 t/Km2/an. D’autres auteurs ont proposé des valeurs plus importantes. (Probst et 
al.1992) ont avancé des valeurs supérieures à 5000 t /Km2/an, ces mêmes auteurs ont estimé 
que la dégradation spécifique moyenne varie de 400 à 610 t/Km2/an pour les rivières qui se 


jettent dans les mers et les océans (Fox, 1993). 
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Chapitre V L'érosivité des pluies 


Chapitre V: L’érosivité des pluies 


V.1 L’érosion & l’érosivité « Historique » : 

L'érosion des sols a été reconnue comme un problème grave aux États-Unis après le 
développement et la mise en culture de nouvelles terres agricoles au cours des années 1800 
et au début des années 1900. Les premières parcelles de mesures ont été établies à 
l’Université du Missouri en 1915 par le chercheur McClure, en 1916 par Vifquain et en 


1917 par Duley-Miller (Figure V.1). 


À la fin des années 1920, une campagne de sensibilisation à la société menée par le 
département de l'Agriculture des États-Unis (USDA), dirigée par Hugh Hammond 
Bennett, a contribué à attirer l'attention du public sur l'érosion des sols en tant que « 


menace nationale ». 


ee ie 
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Figure V.1 Photo prise en 1917, montrant les parcelles d'érosion de Duley-Miller 
(Gantzer et al, 2018). 


A la suite de cet effort, le Congrès des Etats-Unis a affecté des fonds en 1930 à la 
création de stations d’érosion expérimentales dans la moitié est du pays (Bennett, 1939). 
Bennett (1939) a supervisé la mise en place de 10 stations expérimentales utilisant les 
méthodes développées par Miller. Ce nombre a finalement augmenté pour atteindre 49 
stations dotées de parcelles de ruissellement naturel a partir desquelles des données de 


perte de sol ont été collectées. Ces observations et mesures, ainsi que des mesures 
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antérieures, ont fourni des données précieuses pour développer des modèles de prévision 


de l’érosion des sols. 


Cook (1937) a identifié trois catégories d'entités physiques impliquées dans le 
processus d'érosion du sol par l’eau : le sol, l'eau et les plantes, et à partir de cette 
définition, trois variables indépendantes contrôlant le processus d'érosion, à savoir 
"l'érodibilité du sol", "l'érosivité potentielle" et “ la protection de la couverture“. La 
discussion de Cook sur l'érosivité potentielle ne coïncide pas avec les définitions et 
l'utilisation actuelles du terme « érosivité », mais elle prédit le processus de réflexion qui a 
précédé le développement du concept. Cook a défini « l’érosivité potentielle » comme la 
capacité d'une pluie naturelle et de l'écoulement de ruissellement à provoquer une érosion à 
partir d'une zone « standard ». En bref, l'érosivité potentielle de Cook était une mesure de 
la capacité de toute combinaison pluie-ruissellement à causer l'érosion à partir d'une bande 
de terrain unitaire qui monte et descend la pente. Cette idée de zone « standard » était 
analogue ou constituait la base du concept ultérieur de « parcelle unitaire » de Wischmeier 


et Smith (1965), qui sera présenté ci-dessous. 


Cook a également identifié sept facteurs qui contrôlent en grande partie l’érosivité 
potentielle d’une bande de terrain : 1) précipitations totales ; 2) taux de précipitation ; 3) 
vitesses des gouttes de pluie; 4) infiltration caractéristique du sol; 5) capacités de 
stockage de la surface (y compris l'interception) ; 6) les pentes et 7) la longueur de la pente. 


(Nearing et al, 2017) 


Les premières tentatives de création de modèles sont apparues en 1940, Zingg (1940) a 
développé l'une des plus anciennes équations de prévision de l'érosion du sol basées sur 
des données quantitatives. Dans cette équation, l'érosion du sol était liée à la pente et à 
longueur de la pente sur la base de données provenant de cinq sites expérimentaux aux 


États-Unis. Cette équation n'incluait pas de facteur d'érosivité des précipitations. 


Smith (1941) et Browning et al. (1947) ont ajouté à l'équation les facteurs d'érodibilité du 
sol, de culture et de pratique auxiliaire. Musgrave (1947) a introduit une série d'équations 
d'estimation de l'érosion comprenant la longueur de la pente, l'inclinaison de la pente, 
l'érodibilité du sol, un facteur végétal et la quantité maximale de précipitation tombant sur 
une période de 30 minutes au cours d'une tempête. L'équation de Musgrave était la 


première à inclure un facteur d'érosivité des précipitations dans le processus d'estimation 
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de l'érosion du sol. Le facteur d'érosivité des précipitations était défini comme une fonction 
de la loi de puissance, avec les précipitations sur 2 ans et 30 minutes comme nombre de 


base et 1,75 comme exposant. 


V.2 Découverte de l’indice de l’érosivité R (Indice de Wischmeier) : 

Dans l'équation de l'érosion proposée par Van Doren et Bartelli (1956), le facteur de 
précipitations était l'intensité et la fréquence des précipitations de 30 minutes. 
Ces paramètres de précipitations ont fonctionné de manière satisfaisante à certains 


endroits, mais n'étaient pas suffisants pour être utilisés dans l'ensemble des États-Unis. 


Wischmeier et al. (1958) ont corrélé plus de 40 facteurs climatiques ou météorologiques 
(hauteur de précipitation, différentes intensités caractéristiques, les conditions antécédentes 
de précipitations ou d’humidité...etc) et leur relation avec l’érosion sur des parcelles 
standards nues et cultivées dans le sens de la pente. La plus grande corrélation a été 
obtenue avec le produit de l’énergie de la pluie (E) et de l’intensité maximale de la pluie 
pour une période de 30 minutes (130). Le taux de corrélation était d’environ 90% et 
l’introduction d’autres variables accroissait très peu -la précision. De plus, cette nouvelle 
variable (Elzo) s’avérait indépendante des facteurs sols, topographiques, de culture et de 
conservation. Ces valeurs Elzo reflétaient l'effet potentiel combiné des éclaboussures de 


gouttes de pluie et de l'affouillement des eaux de ruissellement sur l'érosion du sol. 


Le facteur d'érosivité des précipitations (facteur R) est défini comme la somme des valeurs 
d'événement EI 30 sur une année. Wischmeier (1959) a analysé environ 8 000 parcelles de 
ruissellement de base, de perte de sol et de précipitations associées, ainsi que des données 
connexes dans 21 États de l’est des États-Unis, recueillies par le Centre national de 
données sur les eaux de ruissellement et les pertes de sol (NRS) à partir des stations 
d’érosion fonctionnant à cette époque. Ils ont confirmé que le facteur R convenait à ces 
endroits, non seulement pour les parcelles en jachère, mais également pour les parcelles de 
cultures en rangées continues, et non seulement pour les variations entre les tempêtes, mais 


aussi pour les variations saisonnières et annuelles (Yin et al, 2017). 


Le facteur R est devenu l'un des six facteurs de l'équation universelle de perte de sol USLE 
(Wischmeier et Smith, 1965). Les cinq autres facteurs sont l'érodibilité du sol, la longueur 
de la pente, l'inclinaison de la pente, les pratiques de gestion de la couverture et les 


pratiques de conservation. Wischmeier et Smith (1965) ont remarqué : "Une différence 
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majeure entre l'équation universelle de perte de sol (USLE) et ses prédécesseurs réside 
dans la manière et la précision avec lesquelles les différences de localisation des 


précipitations sont intégrées dans les calculs de perte de sol." 


V.3 L’évolution du facteur R dans les différentes versions de l'USLE 

L’équation universelle de la perte en sol USLE (Wischmeier et Smith, 1965, 1978) et ses 
versions développés, MUSLE (Williams, 1975), RUSLE (Renard et al, 1997) et RUSLE2 
(USDA-ARS, 2013), ont été largement testés et mises en œuvre dans la pratique de la 


conservation des sols et de l'eau dans une grande partie du monde. 


Il existe certaines différences dans les érosivités pluviométriques entre les différentes 
versions de USLE et les applications de PUSLE et du RUSLE dans d’autres pays et 
régions, notamment la définition d’un événement individuel et le critère de l’événement 
érosif, l’estimation de l’énergie cinétique unitaire à partir de l'intensité, l'estimation de 


l'érosivité due à la fonte et à la fonte des neiges (Tableau V.1). 
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Tableau V.1 : Développement du facteur d'érosivité R dans les différentes versions de l'USLE 


Paramètre 


USLE (Wischmeier et 
Smith, 1965) 


MUSLE 
(Williams, 1975) 


USLE (Wischmeier et 
Smith, 1978) 


RUSLE 
(Renard et al., 1997) 


RUSLE2 
(USDA-ARS, 2013) 


Temps minimum 
entre les 


événements 


Pas de précipitations 


pendant 6 h 


Tempétes incluses 


Toutes les tempétes 


L’érosivité dans 
cette version 
modifiée de USLE 
se base 
uniquement sur le 
débit de pointe du 


ruissellement 


Idem que Wischmeier et 


Smith (1965) 


P <1.27 mm pendant 6h 


Idem que Wischmeier et 


Smith (1965) 


Cumul de pluie >12,7mm ; 


ou avec [\5s> 25,4 mm h`! 


A l'est des montagnes 
Rocheuses (USA) : comme 
Wischmeier et Smith 
(1978) ; À l'ouest des 
montagnes Rocheuses : 


toutes les tempêtes 


Cumul totale de pluies > 
12,7 mm et période de 


retour < 50 ans 


Equation 
(Énergie cinétique- 


intensité) 


e:=0.119+0.0873 logio(ir), 


R= 11.8(QXqy)"™ 
Q : le volume de la 
crue ; qp: le débit 


de pointe 


e:=0.119+0.0873 logio(ir) 
si i < 76 mm h”! ; 


e= 0.283 si i: >76 mm h! 


e=0.29[ 1-0.72exp(0.05i,)] 


e=0.29[1-0.72exp(0.082i,)] 


Erosivité due au 
dégel et a la fonte 


des neiges 


Sous-facteur R, aux 
valeurs de Elso pour la 
région Nord-Ouest : 


R;=1.5P 


Facteur équivalent R requis 
dans les zones de culture du 
NWRR* : 

Reg= -129.0+ 12.61 P 

Reg= 1.602 exp(0.2418 P), 
7.5< P <15.0 


Dans les zones de terres 
cultivées du NWRR* 
basées sur l'équation : 


Reg= -50+ 7.86 P 


NWRR* : Région Nord-ouest du champ du blé (Northwestern Wheat and Range Region) 
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V.4 L’érosivité dans le monde : 

A l’échelle mondiale, les études menées sur le rôle des précipitations dans l’érosion 
montrent que l’indice d’érosivité R subit une grande variabilité spatio-temporelle 
(Renschler et al, 1999 ; Xin et al, 2011 ; Vrieling et al, 2014 ; Naipal et al, 2015 ; Borrelli 
et al, 2016 ; Meddi et al, 2016 ; Borrelli et al, 2017). 


Panagos et al (2017) ont présenté une étude gigantesque sur la cartographie de l’érosivité 
des précipitations dans le monde (Figure V.2). Au niveau continental, l'Amérique du Sud 
connaît l’indice R moyen le plus élevé avec 5874 MJ mm ha’! h? an'!, suivie par l'Afrique 
(3053 MJ mm ha! h?! ant), l'Asie et le Moyen-Orient (1487 MJ mm ha‘! ht! an!). En 


Océanie, l’érosivité moyenne a été estimée de 1675 MJ mm ha"! h! an'!. 


o-100 [MM 1150-1700 
100-200 [ 1700 - 3100 
200-400 [M 3100-5200 
400-700 ME 5200-7400 


Figure V.2 Carte mondiale de l'érosivité des précipitations avec une résolution spatiale de 


30 secondes d'arc (Panagos et al, 2017) 


L'Afrique présente les estimations d'érosivité les plus élevées au niveau national ; Maurice 
et les Comores ont les valeurs d'érosivité annuelles moyennes les plus élevées au monde 
avec une érosivité proche de 20000 MJ mm ha! h! an!. En Afrique occidentale (Liberia, 
Sierra Leone et Guinée équatoriale), en Afrique centrale (R.D. du Congo, République du 
Congo et Cameroun) et à Madagascar, le facteur R annuel moyen est supérieur à 7 000 MJ 
mm ha! h! an!. Ces tendances concordent avec celles d'autres évaluations à l'échelle 
continentale (Vrieling et al, 2014 ; Zhu et Yu, 2015), qui ont indiqué les valeurs d'érosivité 
les plus élevées (> 10000 MJ mm ha-1 h-1 an-1) le long de la côte guinéenne de l'Afrique 
occidentale et centrale, du bassin du Congo et de Madagascar. L'Éthiopie et l'Afrique du 
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Sud ont des valeurs moyennes du facteur R proches de 2500 MJ mm ha"! h! ant, mais les 
configurations spatiales sont très variables, les hautes terres éthiopiennes ayant une 
érosivité extrêmement élevée (> 7000 MJ mm ha"! h“! an`!) tandis que les basses terres ont 
des valeurs 3 à 4 fois plus faibles. Le facteur R moyen le plus faible, avec des valeurs 
inférieures à 115 MJ mm ha’! h! an'!, a été estimé pour le Sahara occidental, la Libye et 


l'Égypte. 


En Asie, le Moyen-Orient présente les valeurs d'érosivité les plus faibles, avec un facteur R 
annuel moyen inférieur à 220 MJ mm ha’! h! an! en Jordanie, en Arabie Saoudite, au 
Koweït, en Syrie, en Iran et en Irak. La Chine a une valeur moyenne de 1 600 MJ mm ha"! 
h! ant, mais présente une forte variabilité avec une érosivité nulle dans les zones arides du 
nord-ouest (désert de Taklamakan), et une érosivité extrême (> 15000) dans les zones 
côtières du sud-est. En plus, la partie sibérienne de la Fédération de Russie et l'ancienne 
Union des républiques socialistes soviétiques (Kazakhstan, Turkménistan et Ouzbékistan) 
ont des valeurs d'érosivité moyennes très faibles (<250 MJ mm ha’! h! an!) en raison de 
leur climat continental. Au contraire, l'Asie du Sud-Est se situe presque entièrement dans la 


classe d'érosivité la plus élevée (>7 400 MJ mm ha! h? an’). 


En Amérique du Sud, le Chili a le facteur R le plus faible avec une valeur annuelle 
moyenne de 1320 MJ mm ha‘! h? an", suivi par l'Argentine (2232 MJ mm ha! h! an’). 
Les autres pays d'Amérique du Sud ont des valeurs d'érosivité moyennes élevées (>3700 
MJ mm ha‘! h? an!), les plus élevées étant celles du Brésil, de la Colombie et de 


l'Équateur (>7 000 MJ mm ha! h! an‘). 


En Amérique du Nord et dans les Caraïbes, le facteur R moyen est de 1 409 MJ mm ha?! lv 
l'an! avec des valeurs très faibles au Canada et dans la partie nord des Etats-Unis, et des 
valeurs extrêmement élevées (>8000 MJ mm ha h! an!) le long des pays du Golfe du 
Mexique et des Caraïbes. La carte de l'érosivité pour les États-Unis montre également des 
valeurs élevées le long du Golfe du Mexique et du sud de la Floride (>8 500 MJ mm ha"! hv 
1 


an‘!), tandis que des valeurs globalement faibles sont observées dans la région du 


Midwest (<690 MJ mm ha"! h! an!) (Renard et al, 1997). 


En Australie, le facteur R moyen est de 1535 MJ mm ha"! h! an! près de 1 767 MJ mm ha- 
1 h-1 an-1 estimés par Teng et al (2016). En termes de modèles spatiaux, Teng et al (2016) 


ont également trouvé des valeurs d'érosivité maximales le long des zones côtières du nord 
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et de l'est (8000 MJ mm ha’! h! an’), qui diminuent vers la région sud-centrale (<300 MJ 
mm ha’! h? an !). En Nouvelle-Zélande, les valeurs d'érosivité élevées (>4000 MJ mm ha‘! 


h“! an`!) se produisent sur la côte ouest de l'île du Sud et diminuent vers l'est. 


La valeur moyenne de l'érosivité pour l'Europe était de 488 MJ mm ha’! h! ant, ce qui est 
très inférieur à celle estimée par Panagos et al (2015) pour l'Union européenne (722 MJ 
mm ha’! h? an‘'). Cela est dû à l'inclusion de la Russie européenne, de l'Ukraine (422 MJ 
mm ha’! h! an!) et du Belarus (365 MJ mm ha! h! an‘!), qui ont tous des valeurs faibles 


par rapport aux autres pays européens. 


V.5 L’ altitude et l’érosivité des précipitations 


Lorsqu'il s'agit de l'érosivité des précipitations, il faut tenir compte de la formation des 
précipitations, de l'érosivité qui leur est associée et des caractéristiques topographiques 
liées aux précipitations d'une région. Le forçage de l'air dû aux obstacles topographiques 
est responsable des précipitations que de nombreuses régions reçoivent ainsi que de la 
création d'un effet d'ombre pluviométrique (Hewitt, 2004 ; Holden, 2005). 

Les relations entre l’altitude et l’érosivité des précipitations ont déjà été largement étudiées 
(Mikhailova et al, 1997 ; Goovaerts, 1999 ; Mello et al, 2013 ; Sanchez-Moreno et al, 
2014; Nearing et al, 2015). Diodato et Bellocchi (2007) ont développé, dans la région 
Méditerranéenne, des modèles statistiques, qui donnent une approximation de l'érosivité, 
basés sur la moyenne mensuelle des précipitations, l'altitude et la latitude. Ils ont conclu 
que le modèle qui utilisait les termes d'élévation fournissait un meilleur ajustement aux 
données d'érosivité estimées que le modèle qui ne tient pas compte le facteur d'élévation. 
En Algérie septentrionale, Meddi et al (2016) ont introduit une étude de cartographie sur la 
variabilité spatio-temporelle de l’érosivité des précipitations, les modèles utilisés pour 
estimer l’érosivité ont basé sur l’altitude, la longitude et d’autres variables. Dans le bassin 
versant de la Macta, Bouderbala et al (2019) ont trouvé que la variabilité de l’érosivité 


dépend essentiellement de plusieurs paramètres, parmi ceux-ci, on trouve I’ latitude. 


V.6 Estimation de l’érosivité dans le bassin versant K’sob 

V.6.1 Les données utilisées 
L’examination de l’érosivité des précipitations dans le bassin K’sob, est basée sur les 
données pluviométriques et hydrométriques fournies par l’Agence Nationale des 


Ressources Hydrauliques ANRH. 


107 


Chapitre V L'érosivité des pluies 


Legend 


Ca Gauging station 


® rainfall station 


drainage network 


HE Ksob Dam 


0 5 10 20 Kilometers 


CS se a a es N 


Figure V.3 Répartition des stations pluviométriques dans le bassin d’ étude 
Le nom, les coordonnées, l’altitude et la période et pas de mesure des stations sont 
présentés dans le tableau (Tableau V.3). Les hauteurs de pluies sont mesurées par le 


pluviomètre ou le pluviographe. 


Tableau V.2 Récapitulatif des données utilisées et des stations associées 


N° Nom de la station de Latitude Longitude Altitude Période Résolution 
St mesure 
a St1 Medjez 35°53'25.78"  4°0379.27" 636m 1973-1994 24h 
J E St2 Bordj Ghdir 35°54'31.36"  4°53'50.45" 1054 m 1973-1994 24h 
2 3 St3 Sidi Embarek 36°06'5.153"  4°54'31.52" 1021 m 1973-1994 24h 
as E St4 Medjana 36°07'44.51"  4°40'26.17" 1042 m 1973-1994 24h 
a 


St5 Bordj Bouarreridj 35°59'24.23"  4°45'35.26" 922m 2007-2015 15min 


1973-1992 Mesures 


8 

8 ; 
8.7 Sto Medjez 35°53'25.78" 4°0379.27" 636 m pence 
g 2 1973-1994 Débits 
A © moyens 

E4 journaliers 

= 


En ce qui concerne les lacunes, les événements de pluies sont rares, et ce n’est pas tous les 
événements de pluies qui sont érosifs, en plus il n'y a beaucoup de lacunes, donc nous 
avons décidé de ne pas prendre en considération ces périodes lacunaires et que ça ne porte 
pas préjudice à l'étude 

V.6.2 Estimation de l’indice d’érosivité R dans la station de référence St5 

Comme on l’a mentionné à l’introduction de ce chapitre, l’estimation de l’indice de 


l’érosivité nécessite des données pluviométriques d’une résolution temporelle fine, au pas 
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d’une demi-heure ou moins. Cependant et de par le monde, de telles séries données sont 
souvent de courte durée, difficile à acquérir et dont le traitement est fastidieux puisqu'il 
nécessite le dépouillement des enregistrements à partir de pluviographes (Figure V.5) 
(Bertoni and Lombardi Neto, 1985). Devant ces difficultés, des solutions alternatives sont 
proposées et où l’indice R est expliqué par des données pluviométriques facilement 
accessibles et couvrantes de longues périodes. En ce sens, l’indice de Fournier, et l’indice 
de Fournier modifié ont été mis en relation avec l’indice de l’érosivité R. Autres relations 
sont basées sur l’indice R en fonction du cumul annuel, mensuels ou journaliers de la 
pluviométrie ou encore sur les cumuls mensuels de pluie pour calculer l’indice de 


l’érosivité. 


Figure V.4 Station automatisée de mesure des données pluviométriques au pas de 15 min 
(Photo : Guesri, 2017) 

En se basant sur les données pluviométriques (au pas de 15 min) de la station automatisé 

St5 (figure V.4), nous proposons des modèles statistiques alternatifs pour estimer 

l’érosivité dans le bassin à partir les données journalières, mensuelles et annuelles de 

pluies. Les meilleurs modèles seront appliqués sur -les pluies disponibles pour calculer 

l’érosivité dans le bassin étudié. La figure V.6 montre les étapes générales suivies pour 


estimer l’érosivité dans le bassin versant K’sob. 
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Houtour de le bague au-dessus du sol : PLUVIOGRAPHE A AUGETS BASCULEURS - MODÈLE M.N. 


{4/88 [41/8] he 1///2///28/ 1/4 d 


EEEE EEE 
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E TER EE EEEa EEEH 
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Figure V.5 Bande d’enregistrement du pluviographe de l’averse du 25.02.2015 enregistrée à la station de Beni IImane 05.05.05 (M’sila) 
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Données pluviométriques 
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f (FT) 
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Estimation de l’érosivité 


Figure V.6 Organigramme de l’estimation de l’érosivité dans le bassin K’sob 
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V.6.2.1 Seuil de la pluie efficace érosive : 

La sélection d’un évènement de pluie érosif dépend du seuil minimal de pluie permettant 
une production de sédiments (Smith et Wischmeier 1978). En se basant sur des essais 
menés sur des parcelles de terrain à petite échelle, Wischmeier et Smith (1965) ont 
considéré une pluie érosive, toute une pluie pouvant produire des sédiments. 

Wischmeier et Smith (1978), trouvent qu’un évènement de pluie est érosif, 1.e., que la pluie 
est efficace en termes de production de sédiment si : (1) la pluie cumulée dépasse le seuil 
T=12,7mm, sans que le cumul ne soit inférieur à 6,35 mm pendant 06 heures ; ou bien si 
(ii) la pluie cumule plus de T/2=6,35 mm en 15 min. Les versions révisées de l’équation 
universelle de la perte en sol (RUSLE, RUSLE2), ont présenté des nouvelles conditions 


(Tableau V.1) 


Les travaux dédiés à la dégradation des sols, montrent que le flux des sédiments drainés à 
l’exutoire de bassins versants est étroitement lié à l’érosivité des pluies (e.g. Hicks et al. 
2000 ; Lal 2001) et que le seuil 12.7 mm présenté par Wischmeier et Smith (1978) est 
souvent excessif (e.g. Bollinne 1978 ; Yin et al. 2017). Denis et al (2013) ont adopté un 
seuil de précipitations érosives de 2,5 mm pour une région semi-aride en Inde. Mannaerts 
et Gabriels (2000) retiennent un seuil de 9 mm pour caractériser les pluies érosives au Cap- 
Vert et Xie et al (2016) estiment le seuil à 9.7mm pour la chine. Selon McGregor et al 
(1995) et Xie and Liu (2002), la baisse du seuil de pluies érosives génère un accroissement 


insignifiant dans la quantification de la production des sédiments. 


De ce fait, un seuil minimal est établi pour chaque saison. Par ailleurs, et vu le caractère 
local des pluies qui caractérisent la région d’étude, le seuil choisi, correspond à la pluie 
journalière minimale mesurées aux stations (St1, St2, St3 et St4) et qui précède une crue 
productive de sédiments de moins de deux ou trois jours. La hauteur minimale de pluies 
érosives d’une saison est identifiée comme le seuil de pluies érosives, T, de la saison 


considérée 
Suivant les saisons, quatre seuils de la pluie érosive efficace ont été déterminés : 


T =3mm, T, 


autumne 


=6 mm, T 


printemps 


=2,5 mm et 7,,=2 mm respectifs aux saisons : 


iver 
Automne, hiver, printemps et été. Les seuils identifiés sont comparés à l’intensité moyenne 


saisonnière des pluies calculée par l’équation : 
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I= P, (V.1) 


SSES 


n 
i=1 
Où : n est le nombre total d’événements pluvieux par saison et P; les hauteurs de pluie aux 


pas de temps de 15min. 


L’évolution suivant les saisons montre une tendance négative (Figure V.7). La 
susceptibilité des sols à l’érosion est minime en hiver, puisque l’intensité moyenne des 
pluies est faible et le détachement des particules du sol nécessite un seuil élevé de pluies. 
Contrairement à cela, en été et à moindre importance l’automne, les sols sont les plus 


susceptibles à l’érosion et les pluies sont les plus intenses. 


Hiver 


Automne 


e Eté 


Printemps 


oOo rFN WwW HRN DW NI 


Seuil de pluies érosives (mm) 


2 2.5 3 3.5 
Intensité moyenne de pluie (mm h`!) 


Figure V.7 Evolution saisonnière du seuil de pluies érosives en fonction de l’intensité 


moyenne des pluies 


L’analyse des intensités de pluies en fonction des crues productives de sédiments 
transportés par l’Oued K’sob a permis de distinguer différents seuils de pluies érosives 
pour les 4 saisons de l’année. Les seuils identifiés sont en concordance avec la nature des 
sols de la région d’étude et du climat semi-aride à influences continentales qui sévit sur les 
hauts plateaux du nord de l’Algérie. Durant l’été et l’automne, un faible seuil suffit car la 


susceptibilité des sols à l’érosion est élevée et les pluies sont intenses. 


Durant la période hivernale ou printanière, un seuil plus élevé est nécessaire car les sols 


sont moins favorables à l’érosion et les pluies sont moins intenses. 


Les seuils saisonniers identifiés pour l’Oued K’sob varient entre 2 et 6 mm. Ces valeurs 


sont faibles en comparaison aux seuils annuels proposés dans la littérature (e.g. Bollinne 
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1978 ; McGregor et al. 1995 ; Mannaerts and Gabriels 2000 ; Xie et al. 2016; Yin et al. 
2017) ou encore au seuil de 12.7 mm établis par Wischmeier et Smith (1978). 


V.6.2.2 L’érosivité d’un événement de pluie Elo 
D’après les critères de Wischmeier et Smith (1978), l’indice de l’érosivité R (en MJ mm 
ha! D, D est donné par la formule suivante : 


n 


=-) ( Ro). (V.2) 


j=1 


ale 


Où : n est le nombre d’années. 
L’érosivité Rp correspondante à une période particulière, Dp, (jour, mois, saison ou année) 
Dp p 


est donnée par : 


Rp, = X Eho) (V.3) 


Où : mj est le nombre d’événements érosifs durant la période Dp. 


L’érosivité d’un événement de pluie qui dure un temps r, est estimé par, E Izo, le produit de 
l’énergie de la pluie, E = e,v,, que multiplie l’intensité maximale pendant une durée de 30 


minutes, /39 en mm/h. 
Où : e, est l’énergie unitaire des pluies (MJ ha ! mm !) et v „le volume de pluie érosive 


(mm). 


L’ énergie unitaire e, : 

Ce paramètre mesure l'énergie nécessaire pour briser les agrégats du sol (l’effet Splash) et 

les transporter ensuite avec par le ruissellement. L'énergie cinétique des précipitations peut 

être calculée si l'on connaît la taille des gouttes de pluie et leur vitesse terminale 

(Wischmeier et Smith, 1958). La mesure de ce paramètre nécessite un disdromètre, qui 

c’est un appareil permettant de mesurer la distribution granulométrique des gouttes de 

pluie par l'analyse de leur énergie cinétique lors de leur impact. Cet appareil est utile en 

relation avec l’imagerie radar (Ellis et al, 2004). 

Les équations empiriques existantes entre l'énergie cinétique et l'intensité des précipitations 
= f(i) ont montré des paramètres différents selon la localisation géographique et la 


technique de mesure (Tableau V.3). 
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Tableau V.3 Relations entre l’énergie cinétique e,. (J m? mm'!) et l’intensité des 


précipitations I (mm h`!) selon la région 


Equation de l’énergie cinétique Région Référence 
(e, en J m? mm, J en mm h’) 
e, = 11.87 + 8.73 log! Amérique du Nord | (Wischmeier et Smith, 1978) 
e, = 8.95 + 8.44 log I Europe (Morgan et al, 1998) 
e, = 9.81 + 11.25 log! Italie centrale (Zanchi et Torri, 1980) 
e, = 9.81 + 10.60 log I Japon (Okinawa) | (Morgan, 2001) 
e, = 29.8 — (127/1) Zimbabwe (Hudson, 1965) 
er = 35.9[1 — 0.56exp (—0.034 J)] | Portugal (Coutinho and Tomas, 1995) 
er = 29.0[1 — 0.6 exp (—0.04/)] | Australie (Rosewell, 1986) 
e, = 36.8[1 — 0.69exp (—0.038 /)| | Hong Kong (Jayawardena and Rezaur, 2000) 
er = 28.3[1 — 0.52exp (—0.042 D] | / (Van Dijk et al, 2002) 
e, = 38.4[1 — 0.54exp (—0.029 I)] | Espagne (Cerro et al, 1998) 
e, = 23.41 — 18 Espagne (Usón and Ramos, 2001) 
e, = 341 — 190 France (Torres et al, 1992) 
e, = 29.0[1 — 0.72 exp (—0.05 /)] | Amérique du Sud | (Brown et Foster, 1987) 
e, = 29.0[1 — 0.72 exp (—0.082 1)] | Etats-Unis (McGregor et al, 1995) 
er = 28.8[1 — 0.45 exp (—0.033 J)| | Panama (Mclsaac, 1990) 
e, = 1(8.95 + 8.44 logio!) Etats-Unis (Brandt, 1990) 
e, = 14.18(i)°172 Japon (Meshesha et al, 2016) 


Vu la non disponibilité du disdromètre, on se limite par la formule standard de Brown et 
Foster (1987) qui est appliquée largement dans l’équation universel de perte en sol révisée, 


RUSLE (Renard et al., 1997) : 


e, = 0.29[1 — 0.72 exp(—0.05i,)] (V.4) 


Ou: i, est l’intensité moyenne des pluies pendant un intervalle de temps (mm h`!) ; 


Entre février 2007 et Mars 2015, environ de 281315 enregistrements de pluie (au pas de 15 
min) ont été effectués au niveau de la station automatisée St5. En tenant compte les critéres 
de Wischmeier et Smith (1978) cités précédemment, et des seuils saisonniers définis, la 
sélection des événements érosifs, est accomplie à l’aide d’un programme informatique 
(voir figure V.8), développé sur Excel VBA (Visual Basic pour Application). Ce dernier, 
est un langage de programmation intégré puissant et sophistiqué qui permet d'écrire des 


fonctions ou commandes dans une feuille de calcul Excel. 


115 


Chapitre V L'érosivité des pluies 


A Microsoft Visual Basic pour Applications - Classeur1 


i Fichier Affichage Insertion Format Débogage Exécution Outils Compléments Fenêtre ? 
@ 1:25, Colt 


(Général) progsource 


B-E Microsoft Excel Obje À 
#8) Feuil1 (Feuil1) 
Feuil2 (Feuil2) 

Feuil3 (Feuil3) v 


Option Explicit 

Dim i, j, k As Integer 

Sub progsource() 

For j = 0 To 2930 

For i= 0 To 95 

Propriétés - Feuil1 Sheets ("feuil2") .Cells(3, 3 + i) = Sheets("Feuill") .Cells(l + i, 2) 
Sheets ("feuil2") .Cells(4 + j, 3 + i) = Sheets("Feuill") .Cells(1 + i 
Feuili Worksheet Next i 

Alphabétique Par catégorie Next j 

End Sub 

Option Explicit 

Dim i, j, k As Integer 

Sub progsource() 

For k 0 To 70 

For j 0 To 11 

For i 0 To 30 


< 


EnableSelection|0 - xINoRestric |, £ 


Figure V.8 Programme de classification des évènements érosifs selon le langage Visual 


Basic (VBA) 


Ces fonctions ou commandes personnalisées peuvent aider à faciliter les tâches et ainsi, en 
utilisant Excel VBA, on peut manipuler tous les éléments d’un tableau, d’une feuille ou 
d’un classeur Excel. Cette extension est largement utilisée en hydrologie pour traiter et 


organiser les données (Hong-jiang, 2002 ; Yang et Niu, 2008 ; Posavec et al, 2017). 
Le tableau suivant présente les valeurs annuelles de l’érosivité, estimées à la station St5 
entre 2007 et 2015, l’indice de l’érosivité R qui est la moyenne de ces valeurs est de 232 


MJ mm ha?! h! an"! 


Tableau V.4 Estimation de |’ érosivité à la station St5 (2007/2015) 


Année 2007- |2008- |2009- |2010- |2011- |2012- |2013- |2014- R 
2008 |2009 |2010 |2011 |2012 |2013 |2014 | 2015 


(Ra) 
MJ mm ha | 300.00 | 317.82 | 280.74 | 189.22 | 112.52 | 230.00 | 224.90 | 203.96 | 232 
lh! an! 
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V.6.3 Les méthodes alternatives du calcul de l’indice R 
V.6.3.1 Indice de fournier 

Fournier (1960) a effectué une analyse de régression entre la production des sédiments 
dans les rivières et plusieurs variables dépendant de la pluie. Il a trouvé une forte 
corrélation entre la perte en sol annuelle et le coefficient de distribution des précipitations 
ou, plus communément appelé indice Fournier : 

Pn 
-P 


Ou: F est l'indice de Fournier, Pp est la pluie du mois le plus humide et P est la pluie 


F 


annuelle totale. L’indice de fournier est devenu un outil pour décrire les degrés du risque 
érosif (Tableau V.5). 
Tableau V.5 Échelle conceptuelle d'évaluation de l'indice de Fournier (Oduro-Afriyie, 


1996) 


Classe Perte en sol (t/ha/an) | Indice de fournier | Risque d’érosion 

l <5 <20 Très faible 

2 5-12 20 - 40 Faible 

3 12 — 50 41 - 60 Modéré 

4 50 — 100 61 - 80 Sévère 

5 100 — 200 81 - 100 Très grave 

6 > 200 > 200 Extrêmement grave 


Dans le bassin versant d’étude, Oued K’sob, la relation qui relie l’indice de Wischmeier à 
l’indice de Fournier est assez faible, que ce soit pour le modèle de régression linéaire 


(R2=0.09) ou pour le modèle puissance (R*=0.05) (Figure V.9). 


350 - a) Modèle linéaire 350 b) Modèle Puissance 


R, = -3.7825 FI + 250.87 


4 R, = 334.81 FI®258 
50 4 R? = 0.0944 


50 9 R? = 0.0531 


N 
S 
© 
/ 
© 
R annuelle ( MJ mm ha-1 an-1) 
> 
© 
Qj 
/ | 


R annuelle ( MJ mm ha-1 an-1) 
a 
© 


© 
+ 


0 10 20 30 0 10 20 30 
Indice de fournier (FI) | f 
Indice de fournier (FI) 
Figure V.9 Modèles R en fonction de l’indice de fournier dans le bassin versant K’sob 
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Cela confirme les résultats trouvés par (Arnoldus, 1977 ; Arnoldus, 1980 ; Meddi, 2013). 
Tandis que d’autres chercheurs ont pu obtenir des corrélations significatives et adopter 
l’indice de fournier comme une solution alternative dans l’estimation de l’érosivité 


(Tableau V.6) 


Tableau V.6 Modèles R en fonction de l’indice de fournier R = f(F) 


Equation de l’érosivité R Région Référence 
(MJ mm ha! h! Dy») 
R = 0.264 F1°° Maroc (Arnoldus, 1977) 
Rm = 2.51 F Espagne (Hernando et Romana, 2015) 
Rm = 0.07397 F+84%7 (Si F<55mm) / (Renard et Freimund, 1994) 
Rm = 95.77 — 6.08 F + 0.477F? 
(Si F > 55mm) 
Rm = 3.76 F + 42.77 Brésil (de Oliveira Junior et Medina, 
1990) 
Rm = 36.849 F10852 Brésil (Morais et al, 1991) 
Rm = 42.31 F + 69.76 Brésil (Silva, 2001) 
Rm = 12.592 F0:603 Brésil (Val et al, 1986) 
Rpm = 68.73 Fo Brésil (Neto et Moldenhauer, 1992) 


V.6.3.2 Indice modifié de fournier (Indice de Arnoldus) 

En raison de la simplicité du calcul de l'indice de Fournier, Arnoldus (1977, 1980) a tenté 
de corréler cet indice avec les valeurs connues du facteur d'érosivité des précipitations pour 
164 stations aux États-Unis et 14 stations en Afrique de l'Ouest. Le chercheur a montré que 
la corrélation entre l'indice de Fournier F et R est faible, donc les résultats n'ont pas été 
satisfaisants en raison du comportement différent de chaque indice. Alors que le facteur R 
inclut toutes les tempêtes érosives, l'indice Fournier ne tient compte que des tempêtes du 
mois où les pluies sont les plus élevées dans le dénominateur (Hernando et Romana, 2015). 
Ainsi, Arnoldus (1980) a proposé une modification de l'indice de Fournier dans laquelle les 
tempêtes qui se produisent en dehors du mois de précipitations maximales augmentent la 
valeur globale de l'indice, obtenant des coefficients de détermination sensiblement plus 


élevés (r? > 0.80) : 


i=1 
Ou: MFI est l'indice modifié de Fournier, P; est la pluie mensuelle et p est la pluie 


annuelle totale. Colotti (2004) a indiqué que l'Organisation des Nations Unies pour 
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l'alimentation et l'agriculture (FAO) a utilisé l'indice Fournier modifié comme estimation 
de l'érosion selon l'équation générale suivante : 

R=a.MFI+b 
Où R est le facteur d'érosivité des précipitations, MFI est l'indice modifié de Fournier et a 


et b sont deux paramètres régionaux d'ajustement. 


Ferro et al (1991) ont déclaré que l'indice Fp, qui représente la valeur moyenne de l'indice 
de Fournier modifié pour un intervalle de N années, est mieux corrélé avec le facteur R que 


les pluies annuelles moyennes et l'indice Fournier modifié (Hernando et Romana, 2015) : 


N 12 2 N 
AOL EX 
F NL| R| Ne j 
j=1 Li=1 j=1 


Où P; j est la pluie (mm) du mois į de l'année j, P; est la pluie annuelle (mm) de l'année j, et 


N est le nombre d'années considérées. Ferro et al (1999) ont également utilisé avec succès 
l'indice d'efficacité énergétique pour estimer le facteur R dans le sud de l'Italie et le sud-est 
de l'Australie. Le tableau suivant présente l’évaluation des risques érosifs en fonction de 


Pindice modifié de fournier. 


Tableau V.7 Classes de l’indice modifié de fournier MFI (CEC, 1992) 


Classe MFI Description 
l <60 Très faible 
2 60 — 90 Faible 

3 90 — 120 Modéré 

4 120 — 160 Fort 

5 > 160 Très fort 


Les tentatives ciblant l’estimation de l’indice de l’érosivité R en fonction de l’indice 
modifié de fournier, sont nombreuses, la tableau V.8 présente quelques exemples des 


modèles R = f(MFI) 


Tableau V.8 Modèles de régression pour l'estimation du facteur d'érosivité R = f (MFI) 


Equation de l’érosivité R Région Référence 
(MJ mm ha! ht D, 1) 
R = 4.46 MF177°30 Algérie (Meddi et al, 2016) 


R = 0.81 MFI} X014 
R=047 MEN Ker ee ag 
X et Z : coordonnées de la station 
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R = 0.297 MFI*-?3 Iraq (Hussein, 1986) 

Rm = 1.05 MFI Italie (Hernando et Romana, 2015) 
Rm = 0.264 MFI? Maroc (Arnoldus, 1977) 

R = 168.42 + 3.27 MFI Belgique (Roose, 1977) 

R = 14.485+0.602 MFI Allemagne (DVWK, 1990) 

R = 42.307 MFI + 69.763 Brésil (Est) (Da Silva, 2004) 

R = —36.46 + 1.84MFI — 0.158P Italie (Ferro et al, 1991) 

Ra = 0.3589MFI*9*62 Chine (Zhang and Fu, 2003) 


La Figure V.10 montre la relation entre l’indice modifié de fournier MFI et l’indice de 
l’érosivité R dans le bassin versant K’sob. Selon le coefficient de détermination (R?=0.34 
pour la régression linéaire, R°=0.27 pour le modèle puissance), la relation est relativement 
significative par rapport à la relation précédente R = f(F), mais elle reste insuffisante 
pour remplacer l’indice de l’érosivité, donc on ne peut pas adopter, les indices FI et MFI 


comme une solution alternative dans le bassin Oued K’ sob. 


350 a) Modèle linéaire 350 a) Modèle puissance 
© o 

— 300 — 300 
2T © 2T © 
2 & 250 2 & 250 
8 À 200 , E ° 
: i i 2 = 200 x 
g = 150 R, = 6.3106 MFI-0.8753 3 = 150 R, = 3.0771 MFI! I83 
Se ioo fe) R? = 0.3455 AC ioo © R? = 0.2768 


Indice modifié de Fournier (MFI) Indice modifié de Fournier (MFI) 


Figure V.10 Relations entre R et l’indice modifié de fournier MFI 


V.6.3.3 Modèles R en fonction des pluies journalières 
De nombreuses corrélations ont été menées pour explique R en fonction de paramètres 
disponibles et plus simple à calculer. Richardson et al (1983) est le 1% qui a suggéré un 


modèle exponentiel liant l’érosivité d’un évènement à la pluie journalière : 
— pp 
Roan = a P°? + €, 


Où a et b sont des paramètres empiriques et € l'erreur normale aléatoire distribuée. 
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Ce modèle est repris par plusieurs auteurs (Avila et Avila 2015 ; Capolongo et al. 2008 ; 
Elangovan et Seetharaman 2011; Karami et al. 2012; Mannaerts et Gabriels 2000 ; 
Nigussie et al. 2014 ; Sanchez-Moreno et al. 2014). D’autres chercheurs ont proposé des 
modèles linéaires (Ali 2015 ; Lee et Lin. 2015 ; Pauwelyn et al. 1988). Yu and Rosewell 
(1996) ont modifié le modèle de Richardson et al (1983) en introduisant des fonctions 
périodiques. PetkovSek et Mikoë (2004) ont appliqué ce modèle pour estimer l’érosivité 
dans le sud-ouest de la Slovénie. Angulo-Martinez et Begueria (2009) et Xie et al (2016) 
ont comparé les différents types des modèles qui se basent sur les pluies journalières, en 


démontrant les avantages et les inconvénients de chaque méthode. 


Au Maghreb, Heusch (1971) propose un modèle en fonction de la pluie journalière 
maximale : R = 143 log(Pan moy: Pan: 107°) + 89.7 


où pest la pluie annuelle. Meddi et al (2016) a introduit pour l’Algérie des modèles 
similaires en se basant sur les pluies journalières maximales. Choukri et al (2016) ont 
proposé des modèles journaliers établis pour chaque mois. Quelques autres exemples des 


modèles, R = f(P;,) sont montrés dans le tableau suivant : 


Tableau V.9 : Modèles de régression pour l'estimation du facteur d'érosivité R = f (P24) 


Equation de l’érosivité R Région Référence 
(MJ mm ha! h! Dp!) 
Roan = a PP +E Etats-Unis (Richardson et al, 1983) 
R = 143 log( Pan. moy: Pan: 1076) + 89.7 | Maroc (Heusch, 1970) 


P2,4n : pluie journalière maximales pour une 
période de retour de 20 ans (mm) 


Ran = 0.87 — 0.037P + 0.016P? Nord du Maroc | (Pepin et Hammouda, 
2012) 
R>4n = 0:03 Pe"? Algérie du (Touaibia et al, 2000) 
nord 
R = 0322, 1% po Algérie (Meddi et al, 2016) 
ee acre ee (2x ( | 3) Sud de l'Italie | (Capolongo et al, 2008) 
Roan = 0.533P + 3.1 Inde (Babu et al, 1969) 
Roan = 0.569P1702 + € Inde (Ali, 2015) 
Ru = 023? Colombie (Avila et Avila, 2015) 
Ros = 0.193P1895 Inde (Elangovan et 
Seetharaman, 2011) 
Reg, = CATR? Tran (Karami et al, 2012) 
Roan = 0.0723 P18 Cap-Vert (Mannaerts et Gabriels, 
2000) 
Roan = 11.64 P — 144.87 Ethiopie (Nigussie et al, 2014) 
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Roan = 2.08 P**” si P<10mm ; Zambie (Pauwelyn et al, 1988) 
Roan = 2.36 P*?" si P>10mm ; 
R4 026 P131 Cap-Vert (Sanchez-Moreno et 
Mannaerts, 2014) 
Roan = à P17263 Chine (Xie et al, 2016) 


a —0,3937 (la saison sèche) ; 
a =0,3101 (la saison humide) 
R34 = 19.14 P Sud de Taïwan | (Lee et Lin, 2015) 

Rog = 0.50 P166 


Dans le bassin versant K’sob, la mesure haute-fréquence de la pluie à la station St5 a 
permis de constituer une base de données journalières associant l’érosivité un évènement 
(R;-El39) avec une pluie journalière P24, sur l’ensemble de la chronique de 7ans (2007- 
2015). Cette base a été utilisée pour caler un modèle de prédiction de l’érosivité journalière 
sous la forme des relations linéaires et puissance. L’évolution de ces deux paramètres est 


présentée pour les différentes échelles temporelles (Figure V.11,12,13). 


a) Fev 2007-Mars 2015 1000 b) Fev 2007-Mars 2015 


5 F150 Qu) eas. : 
Be 8 R2=0,5482 "y 
5 £100 5 = 10 
2 E =E 
BF so OR = 2.3567P-9.6896 SE | 
à À R? = 0.6361 à À 

0 0.1 

0 20 40 60 80 0.1 1 10 100 
Pluie-24h (mm) Pluie-24h (mm) 


Figure V.11 : Evolution de l’érosivité des pluies en fonction de la pluviométrie journalière 


à Péchelle interannuelle : a) modèle linéaire b) modèle puissance 


Au regard des coefficients de détermination, la relation est significative aux échelles 
saisonnières. En effet, pour le modèle puissance, entre 68 et 78% de la variance de 
l’érosivité des pluies est expliquée par la pluviométrie journalière donnant lieu à 
l’événement pluvieux érosif. Pendant qu’on trouve des pourcentages élevés pour le modèle 


linéaire (Entre 79% et 93%) 
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Figure V.12 : L’érosivité des pluies en fonction de la pluviométrie journalière à l’échelle 


saisonnière : a) modèle linéaire b) modèle puissance 
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Figure V.13 : Modèles de régression R = f(P;4n) aux périodes, humides et sèches : a) 


modèle linéaire b) modèle puissance 


V.6.3.4 Modèles R en fonction des pluies mensuelles 
Le développent des modèles qui servent à estimer l’érosivité à partir des cumuls mensuels 
de pluie, était l’objet de nombreuses études. Le tableaux V.10 ci-dessous présente quelques 


modèles de la fonction R = f (Pn) dans plusieurs régions de monde. 


Tableau V.10 : Modèles de régression pour l'estimation du facteur d'érosivité R = f (Pn) 


Equation de l’érosivité R Région Référence 
(MJ mm ha! ht D, 1) 
12 /p2 Maroc (Rango et 
Log R = 1.74.Log 2 (=) + 1.29 Arnoldus,1987) 
i=1 


P; : Pluie mensuelle pour le mois i ; 
P : Pluie annuelle 


LS Portugal (deSantos Loureio et 
R= N > (7.05 Pm10 — 88.92 n10) deAzevedo Coutinho, 
met 2001) 


P19: pluie mensuelle >10 mm ; 219: nombre 
des jours de pluie P>10mm pour chaque mois 
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R = 10 (—-1.48 + 1.48 P,) Allemagne (Rogler et Schwertmann, 
P, : Pluie moyenne en été (Mai-octobre) 1981) 
Rm = 0.33Pt?8 Iran (Karami et al, 2012) 
a Chine (Huang et al, 1992) 
Ra = > 0.0199 (P > 20mm)1°°82 
i=1 
Rm = 0.66 Pnoy + 8.88 Brésil (Oliveira, 1988) 
Rii— Ole E Brésil (Leprun, 1981) 
Rm = 19.55 + 4.2 Proy Brésil (Rufino et al, 1993) 
Rm = 17.97 P Sud de Taiwan | (Lee et Lin, 2015) 
Rm = 0:60:P 27" 
Rm = 0.125 P1629 Chine (Wu, 1994) 
Rm = 3.046P — 26.398 Chine (Zhou et al, 1995) 


Il convient de noter que l’érosivité mensuelle Rm d’un mois est égale à la somme des 
évenements erosifs R; pendant ce mois. Entre Février 2007 et Mars 2015, l’équation de 
régression entre les pluies mensuelles et les valeurs mensuelles de l’érosivité montre un 
coefficient de détermination égal à 0,31 pour le modèle linéaire et 0.45 pour le modèle 


puissance (Figure V.14). 
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Figure V.14 : Modèles de régression R = f (Pn) pour la période Février 2007- Mars2015 


La qualité des modèles (linéaire et puissance) de la relation R = f (Pn) s’est améliorée à 
l’échelle saisonnièrre et à l’échelle des périodes humides et sèches. 
Pour le modèle liénaire: entre 38% et 71% de la variance de l’érosivité mensuelle est 


expliquée par les pluies mensuelle (Figure V.15,16). 
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Figure V.15 : Modèles de régression R = f (Pn) à l’échelle saisonnière 
Le taux de la variance expliqué est plus élevé pour le modèle puissance que pour le modèle 
linéiare (entre 56% et 89%). Un des inconvéients des modèles mensuelles c’est qu’on ne 


peut pas savoir si un mois contient ou pas des pluies érosives. 
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Figure V.16 : Modèles de régression R = f (Pm) en fonction des périodes sèches et 


humides 


V.6.3.5 Modèles R en fonction des pluies annuelles 

Les données annuelles de pluies qui sont largement accessibles, ont été utilisées comme 
substitut de l’indice R. Les modèles qui relient la pluie annuelle à l’érosivité se 
caractérisent pas la simplicité de calcul, c’est-à-dire qu’on peut estimer directement 


l’érosivité annuelle sans évaluer l’érosivité journalière ou mensuelle (Tableau V.11). 


Tableau V.11 : Modèles de régression en fonction des pluies annuelles R = f(P,) 


Equation de l’érosivité R 
(MJ mm hat h! Dy!) 


Région Référence 


R = O.72P PSE 

Ra = 0.001P172 

Ra = 0.35 P(+0.05) pour une région 
montagneuse 

Ra = 0.6 P(+0.05) sur la frange 
côtière 

Ra = 0.5 P(+0.05) pour les autres 
zones 

R = 0210 "PR ENGPE 

q : l’altitude de la station ; P : pluie 


Algérie 
Algérie 
Afrique de l'Ouest 


(Meddi et al, 2016) 
(Touaibia et al, 2014) 
(Roose, 1977) 


Italie (D'Asaro et Santoro, 1983) 
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annuelle moyenne ; NGP: Nombre 

annuelle moyenne des jours de pluie 
Ra = 0.0483 P1:610 

Ra = 587.8 — 1.219P + 0.004105P? 

P<850 mm 


P>850 mm 
Ra =0.5P 
Ra = 0.15 P 
R = —12.98 + 0.0783 P 
Ra = 14.50 P 
Ra = 2.74 Pt? 


R = 2.28P — 8.838 
R = (0.27 P75)/100 
R = 9.17 p920 
R = 23.61 e (00048 P) 
R = 0.0438 P161 
R = 29 + 0.363P 
R = 79 + 0.363P 
R = 0.03 P19 
Ra = 0.028 P1534 
R = 0.028 Pnoy ™”* 
R = 3.48 P + 38.46 
Ra = 0.0534 P16548 
R = 0.0668 P16266 


États-Unis 


Afrique centrale 
Espagne 


Allemagne 
Sud de Taïwan 


La Malaisie 
Etats Unis 

/ 

La Jordanie 
l'Australie 

Le Jharkhand 
Inde 

/ 

Le Chili 


Hawaii 
Chine 


(Renard et Freimund, 1994) 


(Roose, 1975) 
(Hernando et Romana, 
2015) 

(Deumlich et al, 2006) 
(Lee et Lin, 2015) 


(Morgan, 1974) 

(Foster et al, 1981) 
(Cooper, 2011) 

(Eltaif et al, 2010) 

(Yu and Roswell, 1996) 
(Parveen and Kumar, 2012) 
(Singh et al, 1981) 
(Arnoldus, 1980) 

(Bonilla et vidal, 2011) 


(Lo et al, 1985) 
(Zhang et Fu, 2003) 


La figure V.17 montre les résultats de la régression entre les valeurs annuelles de pluie et 


de l’érosivité. Le modèle linéaire a permis d’obtenir l’équation suivante : 


Ra = 1.1836 P — 62.787 


L’équation obtenue pour le modèle puissance est la suivante : 


Ra = 0.0078P1:86 
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Figure V.17 : Modèles de régression en fonction des pluies annuelles R = f(P,) 
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V.6.4 Variabilité spatio-temporelle des pluies et érosivités annuelles 

Dans le but d’estimer l’érosivité dans les stations Stl, St2, St3 et St4, nous avons basé sur 
les modèles journaliers qui sont les plus robustes (en se référant au coefficient de 
détermination) (Figure V.12). On préfère suivre le modèle puissance établi car le modèle 


linéaire est généralement guidé par les valeurs extrêmes (Saporta, 2006) 


Les précipitations et les valeurs interannuelles de l’érosivité estimées aux stations 
pluviométriques (Medjaz (St1), Bordj Ghdir (St2), Sidi Embarek (St3) et Medjana (St4)) 
présentent une évolution temporelle irrégulière mais restent globalement similaires (Figure 


V.18). 
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Figure V.18 Evolution interannuelle de l’érosivité et des précipitations aux 4 stations 


pluviométriques (1973-1994) 
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A l’échelle spatiale, la moyenne interannuelle des indices d’érosivité (228, 386, 342 et 255 
MJ mm ha! h! y`’) est similaire aux moyennes interannuelles des pluies (224, 374, 305 et 
284 mm) estimées aux stations pluviométriques respectives. Cependant, aux stations, 
Medjaz St1 et Sidi Mebarek St3, la pluviométrie et l’érosivité présentent une tendance à la 
hausse. L’accroissement respectif moyen de la pluviométrie est de 1.6 et 3.1 mm/an, soit 
en valeur relative 0.71-1.41%. Ces accroissements restent modérés en comparaison de la 
hausse annuelle moyenne de l’érosivité estimée respectivement à 6.3 et 10.2 MJ mm ha’! h` 
l yl, soit un accroissement relatif de 2.8 et 3%. Pour les stations Bordj Ghdir St2 et 


Medjana St4, aucune tendance significative n’est observée à l’exception d’une légère de 


l’ordre de 1.2% pour l’érosivité à la St2. 


A l'échelle interannuelle, l'érosivité et la pluviométrie présentent des fluctuations 
similaires. Des différences sont mises en évidence à l’échelle saisonnière au regard de la 
densité de l’érosivité estimée par le quotient de l’érosivité par la pluie (Panagos et al. 
2016; Ballabio et al. 2017). Sur les 4 stations, Stl, St2, St3 et St4, les densités les plus 
élevées sont enregistrées durant les saisons d’été et de printemps, où elles oscillent entre 
89-98% (figure V.19). En automne et en hiver, elles fluctuent respectivement entre 87 et 


94% et 66-74%. 


Stl oe 
10_ Automne 28 Hiver Printemps 8 Été 
— 
90 C. 92 92 
% % 
St2 ee Ghdir i 
12 Automne 26 Hiver Printemps 7 Eté 
s _ 
88 — 92 93 
% % 
St3 Sidi Embarek j 
8 Automne 28 Hiver 5 Printemps Eté 2 
% 9 % % 
92 95 98 
% % % 
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Figure V.19 Répartition spatiale et saisonnière de la densité d’érosivité des pluies (%). 


Globalement la densité de l’érosivité reste élevée et semble être étroitement liée à 
l’intensité des pluies annuelles. La variabilité spatiale, l’indice de l’érosivité est 
principalement corrélé à l’abondance des pluies (Figure V.20) dont la densité de l’érosivité 


reste spatialement quasiment similaire (Figure V.19). 
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Figure V.20 Variation spatiale de l’indice d’érosivité en fonction de la pluie annuelle 


Dans le bassin de l’Oued K’sob, la moyenne spatiale de l’érosivité annuelle varie entre 228 
et 386 MJ mm ha! h! y'!. L’indice de l’érosivité est de 302 MJ mm ha! h! y`. Sur le 
même bassin le modèle de Yu and Rosewell (1996), appliqué par Benkadja et al (2015), 
donne des estimations interannuelles 24 à 45% inférieur à celles établies dans cette étude 
(Table 4). 
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En comparaison des estimations de l’érosivité interannuelle du nord du Maghreb, les 
estimations effectuées dans cette étude sont dans la fourchette proposée pour le nord de 
l’Algérie par Meddi et al. (2016). Elles restent modérées par rapport à celles faites pour 
Oued Mina en Algérie (Benchettouh et al. 2017) et Boulabbouz en Tunisie (Kefi et 
Yoshino, 2010). Sinon, les travaux menés sur d’autres régions du Maghreb publient des 


estimations relativement faibles de l’érosivité des pluies (Tableau 4). 


Tableau V.12 Exemple d'estimation de l'indice d'érosivité des précipitations dans le 


Maghreb 


Pays Bassin d’étude R Auteur 
(MJ mm ha’! h! y!) 
Nord d’ Algérie 200-1000 Meddi et al. 2016 
Oued Mina 370-773 Benchettouh et al. 2017 
Bv Bouhamdane 160-240 Bouguerra et al, 2017 
Montures de Sebaa 62-87 Meghraoui et al, 2017 
Chioukh 
Algérie Bv de centre Aurés 47-192 Hassen et al, 2017 
Bv K’sob 125-292 Benkadja et al, 2015 
Oued El Maleh 72-109 Benselama et al, 2018 
Oued El-Ham 45-69 Djoukbala et al, 2018a 
Oued Gazouana 89-130 Djoukbala et al, 2018b 
Oued Sahouat 324-536 Toubal et al, 2018 
Soummam 44-140 Sahli et al, 2019 
Bv Kalaya 95-100 Issa et al, 2016 
Maroc Bv Ourika 55-100 Meliho et al, 2016 
BV de Tahaddart 33-46 Tahiri et al, 2017 
BV Oued El Sourrag 33-38 Medhioub et al, 2017 
Tunisie Bv Lebna 73-95 Gaubi et al, 2017 
Bv Boulabbouz 396 Kefi and Yoshino 2010 


V.6.5 Validation du Modèle d’ estimation de l’érosivité 


Durant la période allant de février 2007 à mars 2015, les pluies aux pas de 15 minutes sont 
utilisées pour calculer l’érosivité annuelle R-RUSLE suivant l’approche établie dans 
l’équation universelle révisée de perte en sol ‘RUSLE’ (Renard et al. 1997). Les valeurs 
trouvées sont comparées à l’érosivité annuelle, R-Modèle, estimées par la présente étude et 


par d’autres modèles issus de la littérature et couramment utilisés. 


Cette comparaison est introduite comme un critère pour apprécier la qualité des modèles à 
estimer l’érosivité annuelle des pluies. A cet effet, une régression simple linéaire est établie 


entre les points R-RUSLE et R-modèle (Figure V.21). 
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Figure V.21 Comparaison des valeurs estimées du facteur R-RUSLE à partir des modèles 


basés sur les pluies journalières 


Un ‘bon’ modèle est celui qui donne une droite ayant une pente proche de 1 et un affixe 
proche de 0. Plus le coefficient de détermination est proche de 1 et plus est le nombre de 
points où l’estimation par le R-Modèle est en bonne adéquation avec l’estimation par le R- 
RUSLE. 


Lorsque le coefficient de détermination est proche de 1, la pente de la droite de régression 
indique si le modèle sous-estime ou surestime l’érosivité R-RUSLE. A partir de ces 
critères, il apparaît clairement que le modèle proposé dans cette étude est celui qui montre 
le plus de performance à s’approcher des estimations de l’érosivité R-RUSLE avec une 
sous-estimation moyenne de l’ordre de 6%. Alors que les modèles établis par Avila and 
Avila (2015) et Elangovan and Seetharaman (2011) sous-estiment respectivement de 33 à 
78% l’érosivité R-RUSLE. Les modèles proposés par Sanchez-Moreno and Mannaerts 
(2014), Mannaerts and Gabriels (2000) et Lee and Lin (2015) surestiment entre 49 et 138% 
l’érosivité R-RUSLE. 
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V.6.6 Erosivité moyenne dans le bassin versant K’sob 

Afin d’estimer l’érosivité moyenne dans le bassin d’étude, nous avons utilisé la méthode 
des polygones de Thiessen (Figure V.22) (Avila et Avila, 2015 ; de Santos Loureio and de 
Azevedo Coutinho, 2001). 
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Figure V.22 Polygones de Thiessen et les stations pluviométriques associées 


Où : 
_ LAR, 
may = 
Avec : 


Rmoy : Erosivité moyenne sur le bassin ; 
R; : Erosivité estimée à la station i ; 
À; : Superficie du polygone associée à la station i ; 


A : Aire totale du bassin. 


La moyenne de surface obtenue par la méthode de Thiessen donne des résultats 
raisonnables pour un réseau peu dense et irrégulièrement espacé. Alors que d'autres 
méthodes telles que la moyenne arithmétique ou la méthode isohyète nécessitent un réseau 


pluviométrique relativement dense et uniformément espacé (Bruce et Clark, 1966). 
V.6.7 Relation de l’érosivité des précipitations avec la production des sédiments : 


La liaison entre l’érosivité et la production des sédiments à l’exutoire du bassin versant est 
examiné dans l’intervalle (1973-1994) à l’échelle saisonnière. L’estimation des apports 


solides saisonniers est basée sur les données instantanées des débits liquides et de la 
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concentration mesurées au niveau de la station hydrométrique Medjez St0 qui se situe en 


amont du barrage K’sob (Voir chapitre IV) 


A une échelle pluriannuelle, l’érosivité des événements pluvieux montre une tendance en 
puissance en fonction de la production des sédiments avec une dispersion assez 
remarquable autour de la courbe de régression (Figure V.23). Seulement 40% de la 
variation de la production des sédiments est expliquée par l’érosivité des pluies. La 


représentativité des modèles en puissance s’améliore nettement à l’échelle saisonnière. 


En automne, 68% de la variance de production des sédiments est expliquée par l’érosivité 


des pluies, contre 58% en été, 53% en hiver et seulement 42% au printemps. 
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Figure V.23 Relation de l’érosivité des précipitations avec la production des sédiments 
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Ces pourcentages sont jugés très appréciables puisque la production des sédiments est 
conditionnée par de nombreux paramètres, autres que l’érosivité des pluies, comme la 
durée et la fréquence des pluies, la variation de la température, l’occupation du sol, les 
pratique culturale et les actions anthropiques (e.g. Walling and Fang 2003 ; Syvitski et al. 
2005 ; Zhang and Nearing 2005 ; Walling 2006 ; Megnounif et al. 2013). 


Les pourcentages les plus élevés correspondent aux saisons d’été et d’automne, où 
l’érosivité des pluies est importante sur les sols secs avec une faible couverture végétale 
(Scott et al. 2006). Au printemps, seulement 42% de la variance de la production des 
sédiments est expliquée par l’érosivité des pluies. En effet, durant cette période de l’année, 
le sol est souvent humide et couvert par la végétation, ce qui les protègent contre l’érosivité 
des pluies (Langbein and Schumm 1958; Seeger et al. 2004). D'autre part, ces 
pourcentages attestent que dans cette région, l’érosvité des pluies est le facteur déterminant 
dans la production des sédiments comme rapporté dans plusieurs études menées sur la* 
production des sédiments sous climat semiaride (e.g. Langbein and Schumm 1958 ; Scott 


2006). 


Conclusion 

Ce chapitre présente l’étude de l’érosivité des précipitations dans le bassin de l’Oued 
K’sob (1480 km?). Dans cette région sévit un climat semi-aride à influence continentale. 
Des modèles sont élaborés pour estimer l’indice d’érosivité à partir de pluies journalières, 
mensuelles et annuelles. L’introduction dans le modèle de seuils de pluies érosives 
distincts pour chaque saison a permis de prendre en considération la variabilité intra- 
annuelle de la pluviométrie et les changements de l’état du sol. Les seuils de pluies 


érosives sont établis en fonction de la production des sédiments. 


Durant les saisons d’été et d’automne, les pluies sont souvent intenses et la susceptibilité 
des particules du sol au détachement est élevée. Durant cette période de l’année, le seuil de 
pluies érosive est faible. Par contre les pluies hivernales sont en moyenne moins intense et 


nécessitent un seuil plus élevé pour provoquer le détachement des particules du sol. 


Un modèle d’estimation de l’indice d’érosivité en fonction des pluies journalières est 
élaboré. L’étude révèle une similarité spatio-temporelle entre l’indice d’érosivité la 
pluviométrie. La relation est étroite durant les saisons d’été et d’automne où les pluies sont 


intenses. De même, l’étude montre que l’érosivité des pluies est le facteur déterminant 
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dans la production des sédiments avec un degré de liaison différent au cours de l’année. En 
automne, 68% de la variance de la production des sédiments est expliquée par l’érosivité 
de pluie, contre seulement 42% au printemps dû au changement de l’état du sol, 


notamment la présence d’une couverture végétale qui protège le sol contre l’érosivité des 


pluies. 
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CONCLUSION GENERALE 


Le processus de production et de transfert des matières solides en suspension dans les 
bassins versants est complexe. Il est étroitement lié à l’intensité des précipitations, à la 
configuration du bassin versant et aux caractéristiques hydrauliques du cours d’eau. Parmi 
les facteurs contrôlant l'érosion du sol et la production des sédiments, l’érosivité des 
précipitations, appréciée par l’indice (R), a une importance primordiale. Elle est utilisée 
largement dans les différentes versions de l'équation universelle de la perte en sol pour 
évaluer l’impact énergétique d’une averse sur le sol. Cependant, l’estimation de ce facteur 
nécessite des données pluviométriques avec une haute résolution (Intensité maximale 
pendant 30min). Or, les données de bonne qualité et sur une longue période sont rares et 


difficiles à acquérir. 


Cette thèse est une contribution dans la mesure de l’agressivité des pluies et ses 
conséquences sur la dégradation des sols et la genèse des sédiments susceptibles d’être 
transportés en suspension par les écoulements de surface en milieu semi-aride. L’étude est 
menée sur Oued K’sob (1480 km?), sous bassin de la Hodna situé au Nord-Est de |’ Algérie 
entre l'Atlas saharien au sud et de l'Atlas Tellien au Nord. 

Le bassin de Oued K’sob est de forme allongée avec un coefficient de compacité KG=1.35. 
Le relief est modéré avec un indice de pente global de 1.03 %. L’exposant climatique, 


établi à partir des courbe IDF, est de 0.75. 


L’aboutissement de la thèse a permis de conclure les points suivants : 


e L'analyse des pluies agressives, c’est-à-dire génératrices de crues transportant une 
charge sédimentaire, a montré une variabilité spatiale du seuil d'érosivité. 

e En comparaison à 12.7 mm seuil d’érosivité établi par Wischmeir et utilisé dans les 
différentes versions de PUSLE, le seuils d’érosivité dans l’Oued Ksob sont faibles, 
ce qui met en évidence un processus d'érosivité typiquement élevé dans les régions 
semi-arides. 

e En été les pluies au-delà de 2 mm sont généralement érosives, alors qu’elles ne le 


sont qu’à partir de 6 mm en hiver. 
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La disponibilité de la mesure des précipitations au pas de 15 min recouvrant une 
période de 08 ans a permis la construction des modèles reliant l’énergie cinétique 
des pluies aux pluies journalières. 

Les modèles établis fournissent un moyen fiable de prédire l'érosivité à long terme 
en utilisant les mesures de pluies facilement accessibles. 

L’érosivité des pluies est le facteur déterminant dans la production des sédiments 


dans le bassin d’étude. 


On recommande de 


Confirmer l’étude sur d’autres bassins limitrophes : 

Généraliser l’étude à un modèle régional permettant d’estimer l’érosivité des pluies 
aux pluies journalières couvrant la région d’étude. 

Élaborer d’autres modèles reliant l’érosivité à d’autres paramètres de pluies 
facilement accessibles et ou mesurables. 

Élaborer des outils d'aménagement des cours d’eau et bassins versants par 
l’ouverture de nouvelles voies de recherche dans l’estimation et la prévision des 
taux d’érosivité pour des périodes de retour données. Cela passera par exemple 


par : 


° Exploitation de nouveaux outils tel que les courbes IDF dans l’estimation et la 
prévision du taux d’érosivité des précipitations pour des temps d’intervalles de 


récurrences donnés. 


° Par l’intégration de nouvelles approches (autre que la régression) dans la 
modélisation, notamment les techniques de l'intelligence artificielle utilisées 


dans les modèles de prévision. 
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